06-361 page | Spring, 2001 


Chapter 1. Course Introduction 


WHAT ARE "UNIT OPERATIONS"? 


A typical process which a chemical engineer might work with is the production of gasoline 
from crude oil 


process = sequence of “unit operations” (physical changes) + chemical reactors 


Any process can be subdivided into a number of steps which are performed in sequence to go 
from some initial starting material (crude oil, in this case) to some final material (gasoline). For 
example, we might start by heating the crude oil to lower its viscosity, then pump the oil to the 
distillation column, where we then separate various components of the crude. A unit operation is 
one of the steps of this sequence. Usually the term refers to steps intended primarily to perform 
some physical transformation (as opposed to chemical transformation) of the input stream. 


Examples of unit operations: 
e heat exchange (change temperature of a stream) 
e fluid flow (transportation) 


e distillation (separation of mixture into multiple streams which are richer in some components 
that original) 


e evaporation (remove water from liquid) 

e humidification (increase water content of gas) 

e gas absorption (remove one component of a gas mixture) 

e sedimentation (separate solid from liquid) 

e classification (divide mixture of particles into different "classes" on the basis of size) 

In this course, we will apply the principles learned in the engineering science courses you have 


had the last three semester to the design of equipment for physical transformations. We might 
divide up the courses in the curriculum as follows: 
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Thermo 1) conservation of energy molecular 
2) criteria for phase/chemical equilibria 

fluids, heat/mass rates of physical processes microscopic 

transfer 


kinetics 1) rate of chemical processes molecular 
2) design equipment for chemical macroscopic 
transformations 


process design divide process into sequence of steps supramacro 
(synthesize flowsheet) 
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Chapter 2. Heat Exchange Equipment 


DOUBLE-PIPE HXER 


One particularly simple  heat-exchange Fluid A inlet 
device is the double-pipe heat exchanger, 
shown in the figure at right. It is called this 
because the basic component is constructed 
from two pipes, one inside the other. One fluid 
(Fluid B) flows inside the inner pipe while the 
second fluid (Fluid A) flows in the annular 
space between the two pipes. Of course, one 
fluid is hotter than the other and heat flows | g^ 
through the wall of the inner pipe from the hot 
fluid to the cold fluid. To obtain larger heat 
exchange area, several pipes are arranged side-by-side and fittings are attached to allow the fluids 
to contact the pipes in series so the overall flow is the same as if only one very long pipe were 
used rather than several shorter ones. 


SHELL-AND-TUBE HXER 


If very large heat exchange areas are required, the double-pipe heat exchanger design 
becomes impractical. Instead, the shell-and-tube heat exchanger design is used. 


Shown above is a bundle of small-diameter tubes which are arranged parallel to each other and 
reside inside a much larger-diameter tube called the "shell", much like strands of uncooked 
spagetti come in a tube-shaped container. The tubes are all manifolded together at either end so 
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that the “tube fluid" enters the left side and is distributed equally among all the tubes. At the 
right side, the fluid exits from each tube, is mixed together in a second manifold, then leaves as a 
single stream. The second fluid, called the “shell fluid", flows in the space in between the 
outside of tubes. Baffle plates inside the shell force the shell fluid to flow across the tubes 
repeatedly as the fluid moves along the length of the shell. 


In the design above (called a 1-1 hxer), the flow is from left to right inside the tubes and from 
right to left for the shell-side fluid. Another common design is shown below (called a 1-2 hxer). 


orn 


FIGURE 15.3 
1-2 parallel-counterflow exchanger. 


Half of the tubes have flow from left to right while the other half have flow in the opposite 
direction. 


FIGURE 15.4 
2-4 exchanger. 
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Chapter 3. Fundamentals of Heat Transfer: Review 


LoG-MEAN DRIVING FORCE 
The design equation used for sizing heat exchangers generally takes the form of 
q =U AAT (1) 


where q is the heat duty of the exchanger (rate of energy transfer) U is some heat-transfer 
coefficient, A is the area of the heat exchange surface, and AT is the driving force for heat 
transfer. 


In any real heat exchanger, the local driving 
force AT will vary along the length of the 
exchanger: say it varies between AT! at one end 
of the exchanger to A7» at the other end. Then 
the question arises: which value of AT do you 
use in [P The answer turns out to be: use the 
log mean of AT; and AT»: 


This turns out to be correct in nearly every 
situation. We will prove it now for one particular 
case and leave some other cases to be proved in 
the homework (HWK #1, Probl. 1) 


Proof: When the driving force varies with position along , one slice d i 
the exchanger, the usual approach to math modelling is : 
to divide the length of the hxer (L) into a large number 
of thin slices, such that the thickness of each slice is thin 
enough so that the temperature of either stream is 
essentially constant over the slice, leaving the driving 
force constant over the slice: 


AT = T(x) - T(x) 


Depending on which x we choose, we get a different driving force. The rate of heat transfer from 
the hot fluid to the cold fluid in any slice is given by a local heat transfer coefficient U(x): 


dq = U(x) AT(x) dA (2) 
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where q(x) is the rate of heat transferred from the hot fluid to the cold fluid in the length of hxer 
from x=0 and x=x; and dq is 


q(x*dx) = q(x) + dq 


and dA is the heat exchange area in this slice of length dx. If the inner pipe (of a double-pipe 
hxer) is a cylinder of diameter D, then the heat exchange area of this slice is 


dA — xD dx 


pipe wall 


Integrating as a function of x is difficult since 
we'd have to predict A7(x). There is a easier way 
(but a “trick” is required). Consider steady-state 
counter-current flow through a double-pipe heat 
exchanger (shown at right). A steady-state energy 
balance on the hot fluid gives 


A A A 
rate in = rate out 4 
m'H' -m'Hi +q (3) x=0 X=xX x=L 
A similar balance on the cold fluid gives 
rate in = rate out 
(4) 
mH, +q = mH 
Solving (3yland laylror q:m |H (x)- H; | =q(x)= m'| H' (x)- Hi | (5) 
eR —wW eR 
rate of energy rate of energy 
gained by lost by 
cold stream hot stream 


where m is the mass flowrate of the stream (units are Kg/sec), H is the enthalpy per unit mass 
(units are J/Kg), the subscripts and primes have the following meaning: 


e unprimed refers to the cold stream, 

e prime refers to the hot stream, 

e “1” refers to the left end of hxer, and 

e "2" refers to the right end of hxer. 

Now, let's try to relate each enthalpy to the appropriate temperature. In general, enthalpy 


changes can have both /atent heat (caused by change in phase) and sensible heat (caused by 
change in temperature) contributions. 
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dH = cy dT + Xdy 
—— 


: lat 
sensible latent 
heat heat 


where 

e dH is the change in enthalpy (energy/mass) 

* Cp is the heat capacity (energy/mass/degree) of the fluid, also known as "specific heat" 
e dT is the change in temperature (degree) 

e i is the heat of vaporization (energy/mass) 

e dy is the mass fraction of fluid vaporized 


When the fluid consists of two phases (liquid and vapor), which have different heat capacities 
CpL and Cp y, we use weighed average of the two heat capacities 


Cp = (l-y)cpL + Y €py 


where y is the mass fraction of vapor present and 1-y is the mass fraction of liquid present. For 
the time being, let's assume: 


Assumption #1: no phase changes occur in either fluid 

Now we only have to worry about contributions from changes in the temperature (sensible heat): 
dH = c, dT 

Integrating over the temperature change allows us to evaluate the change in enthalpy 


using using 


assumption assumption 
#1 Tj #2 
AH = JepaT = — e,(%-T%) 
T, 


In general, the specific heat (c,) depends on temperature. But if the temperature change is not too 
large — say 10-20? — then the specific heat doesn't change too much. Let's assume we can 
neglect these small changes: 


Assumption #2: constant specific heat for both streams 
Then we can factor c, out of the integral and the integration becomes trivial: 


AH =c, AT 
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In particular, H-H = cy(T-T|) 

and H'- H* - cy (T' - T1) 

where c, and c’, are the specific heats of the cold and hot streams. Thus (5)]can be rewritten as: 
at x=x: mcy(TG) - T) = mc, (TG) - T4) = qo) (6) 
For example, 

at x=0: mcy(T, - T) » m'es(T'4 - T1) = «(0 «0 


In particular, we are interested in the total rate of heat transfer over the entire length L of the 
exchanger: 


at x=L: mc p(T -n)em'e;(T* -1\)=q(L) =4r 


Now let me show you the "trick": using ko) I can show that the temperature of either stream [i.e. 
T(x) or T (x)] is a linear function of the amount of 
heat transfered, q(x): 


T(x) =F, | : je 


MC p 


AR 
T'(x) = +| —— |g) 
m 


Cp 


AT 
Since both temperatures are linear functions of q, 

their difference is also a linear function of q. A 
straight line has the very nice property that its slope 

is the same at every point: 


d(AT) onst- 22-4 (7) 
dq qr —9 

I can easily evaluate the driving force (AT) at the two ends of the heat exchanger, and I know the 

total heat duty (q7), so I can calculate this slope. The dg and d(AT) which appear in K7)|represent 

the incremental rate of heat transfer and the change in driving force AT which occurred over a 

different length dx of the heat exchanger(from x to x+dx). They are related by 


dq = U dA AT (8) 


where dA is the heat exchange area in that slice of the heat exchanger from x to x*dx. Treating 
the heat exchange surface as a cylinder of radius R: 
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dA =2nR dx 


lotus n d(NT) _ AD -AT 


U(AT(dA) — ar 


d(AT) AT -AT, 
AT qr 


or 


(44) 


Next, we integrate over the length of the heat exchanger from x=0 (where AT = AT!) to x=L 
(where AT = AT»): 


AT, x= L 

T = 
J ae » B9 5d JU(x):2xR dx 
aw OE qT x dA 


The integral on the right-hand-side becomes trivial to evaluate if we can make one additional 
assumption: 


Assumption #3: U is constant 


Then we can factor U outside the integral and integrate the remainder: 


x=L L 
J2nRdx = 2nR [dx =2nRL = Ar 
x=0 0 


which is just the total heat exchange area. Thus we have 


AT, AT, - AT, 
1 2.27752 1 UA 
ATi qT 
AT, — AT, 
Solving for the total heat duty qr = UAL : 1 
AT, 
AT, 
— — AT, — AT, 
or qr -UAT AT, where AT, = : L 
l AT, 
n——- 
AT 


is called the Jog-mean temperature difference (LMTD). It represents the average driving force 
for heat transfer. This result provides us with the main design equation for heat exchangers. 
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Summary of assumptions: 


e Uand c,'s are constant 


Spring, 2001 


e AH of neither stream has contributions from both sensible and latent heats (i.e. AH 


is either all sensible heat or all latent heat) 


Actually, we assumed in the analysis above that AH was all sensible heat for both streams. In the 
homework, you will show that AH can also be all latent heat and the LMTD is still the 
appropriate average driving force. For example, if you use condensing steam as your hot fluid, as 


long as the steam remains saturated: 


for saturated steam: T (x) = const = T», 


We also assumed in the analysis above, 
that the flow was counter-current. You 
will also show in the first homework 
problem that the same design equation 
applies to co-current flow. 


RESISTANCES IN SERIES 


Consider heat transfer through the 
inner pipe of a  double-pipe heat 
exchanger. Let's take a slice out of the 
heat exchanger and examine the 
temperature profile along a radius to the 
pipe at one particular axial position. 


aL 


oblique view 
of slice 
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One can identify three regions in which the temperature varies sharply. These correspond to 
the three resistances to heat transfer: the pipe wall itself, the fluid film inside the pipe and the 
fluid film outside the pipe. These resistances act in series. The total resistance is related to the 
resistance of each film by: 


1 1 x, D 
U, ho k„ D, | 
total resistance of resistance of end VISW 
resistance outer film pipe wall of slice 

1 Do 

h; D; 

yoa 

resistance of 
inner film 


That resistances add in this way is a very 
important concept, so I'm going to take 
the time to prove it to you again. This 
will also give us an opportunity to 
review the definitions of the quantities. 


Consider the temperature profile 
along some line drawn perpendicular to 
the pipe wall through the center of the 
pipe. These temperature variations drive 
heat transfer down the temperature hill. 
The rate dq of heat transfered through 
some small element of the inner surface 
of the pipe having area dA; is 


resistance #1 


dq 2h; dA; (Ty - Twn) (9) 
a) 
nD; dL 


where h; is a single-phase heat-transfer coefficient for the inner film of 
fluid next to the pipe and dL is the length of the pipe segment. At steady-state, this same rate of 
heat transfer must occur through the liquid film outside the pipe: 


dq=h, dA, [Ps 
—— 
xD, dL 


z) (10) 


At steady-state, this same rate of heat transfer must also occur through the pipe wall itself: 


k. — 
dq=—“dA, (Ton = Tye) (11) 
Xw 
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where kọ is the thermal conductivity of the pipe wall and xy, is the wall thickness. Since the 
inner and outer areas are somewhat different, we use an average of them in 


"E cud. D, — D; dL 
where dA; = dham dd oa - i) 


L= = 
dA D 
inf o 4 4] nf 7» 


is the log-mean of dA,, and dAj. Now we define the overall heat-transfer coefficient Ug as 


- xD, dL 


dq =U, dA, (Tj, — To) (12) 


This is called “overall” because we are using the overall driving force in the definition. Solving 
(12)|for the overall temperature difference, then writing this overall as the sum of the individual 
temperature differences 


dq 
= T, To =T; — Thy + Thy — Tew + Tew — Te 
Uo dA, — ———— 
dq dq dq 
Xw 
Substituting and 
dq__ @ , dq  , dq 
Uo dA, h; dA; kw Gay ho dA, 
Xw 
Dividing through by dq/dA, leaves: 
1 " ] D, " M Do n 1 
U, h; Di ky Dy ho 


after substituting the expressions relating the dA’s to D’s and dL. 


Heat transfer coefficients (either U or A) are like conductances in electricity: the heat flux is 
proportional to them (like electrical current through a resistor is proportional to its conductance). 
The reciprocal of the heat transfer coefficients (1/U or 1/A) is like an electrical resistance. The 
formula above (relating 1/U to 1/h's) is like that relating the overall resistance of a circuit 
composed of three resistors in series: the total resistance is the sum of each of the three 
individual resistances. 
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FOULING FACTORS 


The usefulness of analogies is that we can extend our knowledge in one area using 
knowledge in other. For example, if there are additional layers through which heat must be 
conducted, then we can just add their resistances to get the total. A situation involving additional 
layers is: 


fouling — deposition of dirt, scale or any solid deposit on one or both sides of the tube 


inner fluid film 


inner fouling layer 
pipe wall 


outer fouling layer 


outer fluid film 


Fouling generally reduces heat transfer by imposing an additional resistance. It's effect is 
accounted for by use of 


fouling factors: hg = -+ 


There may be a fouling factor for both the inside and the outside of the pipe. Adding these 
resistances to the others: 


lao. gi. aae Dy e Jobs 
U, hy hdo ky Dr; hii Dj h; D; 
—3 —— 
fouling fouling 
outside pipe inside pipe 


By means of the electrical analogy, we are able to guess the form of the correction without 
having to go through the analysis. 


CORRELATIONS TO ESTIMATE A’S 


In all the problems we have worked so far, we have been given the values of the 1-phase heat 
transfer coefficients (the h; and h,). Now we are going to turn our attention to how these can be 
estimated if they are not otherwise known. The empirical correlations appropriate for a given 
situation usually depend on three questions: 


1. What is the geometry? (e.g flow inside pipe, flow across bank of tubes, flow around sphere) 
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2. Isthe flow turbulent or laminar? 

3. Does a phase change occur? 

Our book presents a number of common correlations which cover the most common situations: 


Sect. 4.5, 4.6 - without phase change 
laminar flow inside tubes — p238 (not covered in class) 
turbulent flow inside tubes — 239-2444 


Sect. 4.8 - with phase change (not covered in class) 


For flow inside long tubes, the criteria for turbulent flow is simple: just calculate the 
Reynolds number: 


Di GD 
Nee = p(v) bo l 
Hu u 
where p = fluid density 


u = fluid viscosity 
(v) is the cross-sectional average fluid velocity 


The latter quantity is defined as (v) = 
i 
where Q = volumetric flowrate through tube (units of volume per time) 


S; = nDj?/4 = cross-section area of tube 


Another term which is sometimes used to calculate Reynolds number is the mass velocity: 


G=p (v) = eo Pa (units of mass per area per time) 
i i 


The mass velocity is often used in problems involving gas flow where both density and velocity 
depend on pressure, but their product is independent of pressure (since it represents mass 
flowrate and at steady state, a mass balance requires the mass flowrate to be the same at all axial 
positions). However Reynolds number is calculated, the criteria for turbulent flow is as follows: 


NRe > 6000 = flow is turbulent 
NRe < 2100 = flow is laminar 
2100 « NRe < 6000 = flow is in transition 


* You should make a special effort to read these pages before working the HWK problems, since I don’t have time 
in class to present all the details. 
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hj for Turbulent Flow In Long Pipes 


For turbulent heat transfer in long pipes, the individual film coefficient inside a pipe can be 
estimated using the Sieder-Tate equation [eq. (4.5-8) in G]:* 


0.8 1/3 0.14 
2 = 002 028 ER | T 
u m 
N Nu N N k 
Re Pr 
where k — thermal conductivity of fluid inside tube 


cp = heat capacity (specific heat) of fluid inside tube 
u = viscosity of fluid inside tube 


Unless specifically indicated, all fluid properties are evaluated at the mixing-cup temperature of 
the fluid: 


u, k, Cp evaluated at T' 


€. 


The subscript “w” indicates the property is to evaluated at the wall temperature: 


uy evaluated at T, 


Experimentally, one finds that heat transfer coefficients can be different for heating and cooling 
over the same temperature range. The reason, is that the viscosity at the wall is very important in 
determining transport and py, can be very different from u. This is accounted for in the wu 
factor. 


Which T do you use for property evaluations? 


When the inlet and outlet 7"s differ, what T do you use to evaluate the properties at? Use 
arithmetic mean values for T and T: 


T, +T = Toa +T, 
a b and T, = = wb 


T= 
where the subscript “a” refers to the inlet and “b” refers to the outlet. Since the outlet T and the 


outlet Tọ depend on the heat transfer coefficient, this makes the determination of / a trial-and- 
error (see Example 4.5-2). 


* The coefficient of 0.023 in this equation is the value which appears in McCabe, Smith & Harriott. Geankoplis 
gives a value of 0.027 instead. In homeworks and exam, you can use either value. The Notes use 0.023 because 
most of the homework problems using this equation were taken from MS&H. 
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Chapter 4. Shell-and-Tube Heat Exchangers 


The value of A, for shell-and-tube heat exchangers depends on arrangement of tubes and 
baffles inside the shell. We will first describe typical arrangements and then give a common 
correlation for ho. 


TUBE PITCH 


The bundle of tubes in a shell-and-tube heat exchanger can be stacked in one of two different 
Ways 


OO WOO eo O 
eo O OO O 
eeoeoo 000000 


triangular pitch square pitch 
(advantage: closer packing) (advantage: easier cleaning) 


The use of triangular pitch allows the tubes to be more tightly packed -- more tubes and therefore 
more area per unit volume of shell. This makes the shell cheaper. 


On the other hand, the square pitch has the advantage that it is easier to clean. You can get a 
long brush in between two rows or between two columns. Of course, the shell has to be pulled 
off of the tube bundle before you can get at it in either 


case. 
Regardless of which geometry is used, the center-to- r 
center distances between adjacent tubes is called the P N p [| p 
tube pitch. 

P p 

5 i 
BAFFLES pez); p= Do + 1 inch 


To some extent, it is advantageous to reduce the 
cross-sectional area available for flow. For the same 
volumetric flowrate, this will increase the mass flux or the mass velocity of the fluid flow, which 
will increase h,. Of course, this also increases the pressure drop and the pumping cost. 
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A common technique for thus increasing o is to 
install baffle plates, which partially block the ^ baffle window 
cross-sectional area on the shell side. A baffle f 
plate is a disk whose diameter is the inside 
diameter of the shell: 


D, = diameter of shell 


with holes drilled through it for the tubes to pass 
through. Part of the disk is cut off to form what is 
called the baffle window, which is the cross 
section available to the shell-side of the fluid. 
When the height of the window is 1/4 of the shell 
diameter, we call the plate a 25% baffle. 


A number of these plates is placed along the length of the exchanger with the location of the 
baffle window alternating between the top and bottom of the plate. 


4 baffles 
tube 
side 
flow 


tube sheets ^ 


shell-side flow 


A second advantage of these plates is that they prevent the formation of large stagnant regions 
(having low heat transfer rates) which usually accompany “channelling” of the fluid in a direction 
path connecting inlet to outlet. The spacing between the baffle plates is called the baffle pitch: 


P — baffle pitch 
Typical baffle pitch is a fraction of the shell diameter: 
0.2 < P/D, « 1 
Besides reducing the area available for flow on the shell side, the baffles also force the fluid 
across all of the tubes, which normally will enhance the heat transfer coefficient. 
ho FOR SHELL-AND-TUBE EXCHANGER 


Reference: McCabe, Smith & Harriott, 5th ed., p432-436 
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An approximate but generally useful correlation for estimating the shell-side coefficient in a 
shell-and-tube heat exchanger is the Donohue equation: 


0.6 0.33 0.14 
C 
h, D, =02{ 222e) | e | u ) a3) 
k u k Hw 


where u, k and Cp are properties of the shell-side 
fluid. The form of this equation is very similar _ Shell flow 
to the Sieder-Tate equation, except the / Ln ube 
exponents and coefficient are different. Another 
difference is that the mass velocity Ge is a 
weighted average 


Go =./G,G, (geometric mean) 


shell flow L P d 


of the mass velocity through the baffle window tangent to tubes 


G, =— 
b F 


and the mass velocity for crossflow perpendicular to the tubes 


where Sj is the area available for flow of the shell-side fluid through the baffle window and S, is 
the intersticial area available for crossflow perpendicular to the bank of tubes at the widest point 
in the shell. 


2 2 
S; - af 2) mE) (14) 


4 4 
area of 
blue region 


^ baffle window 


where — D, — diameter of entire shell 
Np = number of tubes in baffle window 
fp = area fraction of baffle plate that is window 


For a 25% baffle plate, we have 
fp = 0.1955 (15) 


The second area is calculated from 
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D 
s= rofi) 
P 


p = tube pitch 
P = baffle pitch 
Appendix: Derivation for Sp and Sç 


In case you are curious as to how these formulas we obtained, I have put their derivation in 
this appendix. 


widest 
row of 
tubes 
P 
< > 
A 
ba 
v 
side view area for end view 
cross-flow 
(shaded) 


Sc is the area open to fluid flow across the widest row of tubes in the shell and between two 
adjacent baffle plates. First the total area (1gnoring tubes) is 


total area = PD, 


where P is the baffle pitch. To calculate the fraction of PDg which is open. let's subdivide Dg 
into repeating unit cells: 


length of unit cell = p (tube pitch) 
length which is open = p-Dg 
fraction open = (p-Do)/p = 1 - Do/p 
Multiplying this fraction by the total area gives the area open: 


S, = Pp 1-22) (16) 
p 
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Of course, we have ignored some end effects which could be significant if the number of tubes is 
not large. 


Likewise the open area of a baffle window is the total area of the window minus the area blocked 
off by tubes: 


2 2 
TD TD 
s-^| 3 zl 2) 


where fp = 0.1955 (for a 25% baffle) ~- sector 
"4 (heavy outline) 
and Np = no. of tubes thru window // segment 


dark shading) 
To get this value of fp, we have to do a little 


geometry. First, recall that the area of a circle of 
radius r is x72. 
= nr? 


A 


circle 


A “sector” is a portion of a circle (see figure at 


right): 
1,2 
Asector = af (20) 
total angle 
subtended 
by sector 


where 20 is the total angle subtended by the arc of the sector. Note that if we substitute 2x for 
the 20, we get the Asector = Acircle. The area of the triangle whose sides are x, y, and r equals 


—1,,_1 
triangle = 5 bh — 7Y 


Finally, the area of the segment is calculated from 


Asegment = Asector ~ 2A triangle 


For a 25% baffle, the height of the segment, r-x, will be 25% of the diameter -- or 50% of the 
radius: 


l(rY 43 3» 2 
Atriangle = | j j- 3 r^ =0.217r 
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2 T 2 
Asec tor 7 r^ = a 
Asegment = 7 - td i E = zp = 0.614187? 


Now fp is the fraction of the total area of the circle which is represented by this segment: 


E 232 
A segment 3 4 


h= - =— = 0.1955 
A rele nr? 3.1416 


MULTIPASS CONSTRUCTION 


Just as we can improve the heat transfer coefficient on the shell side by reducing the baffle 
pitch, we can improve the heat transfer coefficient on the tube side by restricting the flow to half 
of the tubes. Both have the effect of reducing the area available for flow, thereby increasing the 
fluid velocity and the heat 
transfer coefficient. 
tube 
side 
flow 


Suppose we partition the 
manifold in such a way 
that the inlet fluid can only 
enter half of the tubes. 
The other half of the tubes 
are used for return flow. 
We say that such an 
exchanger has fwo tube 
passes. 


tube sheets 


shell-side flow 


shell-side flow 


Although two tube passes 
improves hj, one of the passes 
will be co-current A 2-4 
exchanger improves both hi 
and fg, and in addition, the 
flow is more nearly 
countercurrent. We could also 
partition the shell side with a 
horizontal baffle in addition to 
the vertical baffles, producing 
two Shell passes. 


tube 
side 
flow 


The most commonly used 
designs are: tube sheets 
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e 1-2 exchanger 
e 2-4 exchanger 


where the first number is the number of shell passes and the second number is the number of tube 
passes. 


An even number of tube passes is desireable because then input and output of the tube-side fluid 
occur on the same end of the exchanger. This allows the other end of the exchanger to "float" in 
response to thermal expansion of the metal, which can otherwise be a serious problem when the 
shell changes temperature during startup or shutdown. 


"TRUE" MEAN TEMPERATURE DIFFERENCE 


For multipass construction, the flow is neither counter-current nor co-current. As a consequence, 
the average driving force is no longer the LMTD. As it turns out, what is called the true mean 
temperature difference (MTD) is somewhat less than the log-mean: 


ATm = FTATL (17) 


This is the AT, which when multiplied by UA gives the total rate of heat transfer for a shell-and- 
tube exchanger: 


gT = UoAoATin 
F ris a correction factor which is unity for pure co-current or pure counter-current flow: 
Fr=1 


For pure crossflow to a bank of tubes, or for a mixture of crossflow and counter-current flow, the 
factor is less than unity: 


Fr= Fr (Y, Z, type of flow) < 1 


where Y and Z are ratios of temperature differences: 


Y= Teb E Tra and Z= Tha B Thb 
Tha = Teb Teb z Tea 


where the subscript “c” refers to the cold fluid while the subscript “h” refers to the hot fluid and 
the subscript “a” refers to the inlet while “b” refers to the outlet. 
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Fr is usually plotted as a function of Y for different 
values of Z. Such a plot is sketched at right. More 
accurate plots can be found in Geankoplis's book: 


1-2 HX > Fig. 4.9-4a 
2-4 HX Fig. 4.9-4b 


Note that as Z0, Fr — 1 over the entire range of 
Y«]. This special case is applicable for condensing 
vapors. As long as the vapor is condensing, its 0 Y 
temperature stays at the boiling point: 


Tha = Thb = Thp 
so that Z=0Oand Fr= 1 


for condensing steam as the hot fluid. Similarly if we have boiling water as the “cold” fluid, its 
temperature doesn’t change (i.e. Teq = Teb = Thp) and then we have 


Y=Oand Fr=1 


Although the flow in the exchanger is complex, we calculate AT z in lazylas if the flow were 
counter-current. The deviation caused by not being counter-current is included in the values of 
the empirical correction factor FT. 


Example: A steel mill is considering whether a spare heat-exchanger will be able to cool of 3000 
gal/hr of a dilute acid from 250?F to 180?F using 5000 gal/hr of 80?F river water. You are asked 
to estimate: a) the true-mean driving force, AT, and b) the shell-side coefficient, ho. 


Description of available heat exchanger: 


type: 1-2 

tubes: 158 1-inch, 14BWG stainless tubes, 16 ft long, stacked in 11⁄4 inch square 
pitch 

baffles: 25% with 6 inch pitch 

shell: 20 inch inside diameter, mild steel 


Solution: First we have to determine the outlet temperature of the cold stream. The rate of 
energy lost by the hot stream is the rate of energy gained by the cold stream 


gy TC ph (Tha = Thap) = TeC pe (Top = fa) 
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Assuming the physical properties of dilute acid and river water are the same as pure water, we 
can take cph = cpc so that heat capacity can be cancelled out of the second equation above, 
leaving 


Top -Tea _™hCph my _ 3000 gal/hr _ 
lh = Thb McC pe Me 5000 gal/hr 


(18) 


From the problem statement we have 
Tha =250°F Tpp=180°F Teq= 80°F 
Substituting these into [18)]and solving for Tep, we obtain 
Tep = 122°F* 


Next, we calculate the two quantities needed to evaluate the true-mean driving force: 


T,, —T, T.» — T, 
gaa Tib. 1 .167 and yo tee le 29247 
Teb = Tea 0.6 Tha = Toa 


Since this is a 1-2 exchanger, we read FT from Fig. Fig. 4.9-4a: 
Fr= 0.97 


Although the flow in the exchanger is complex, we calculate AT z as if the flow were counter- 
current: 


AT Tggq-Tep = 128°F and AT = Tpp-Teg = 100?F 


— AÑ AT 
AT, =— 1 2 = 113°F 
AT, 
In —. 
AT» 
Finally AT, = Fr AT 1 =110°F 


The shell-side coefficient (họ) can be estimated using Donohue’s correlation Since acid is 
corrosive, we would be well-advised to place it inside the stainless tubes, leaving river water on 
the shell side (which is mild steel). With this choice, the mass flowrate is 


* Such a high temperature in the river could be harmful to fish and other wildlife. Engineers need to be mindful of 
the impact of their designs on the environment as well as on the health and safety of people in the vicinity. Faced 
with this result, the engineer should suggest altering the design (e.g. increase the flowrate of cooling water to obtain a 
lower increase in temperatue). For the purpose of this academic example, we will ignore this objection and proceed 
through the analysis. 
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DES (s000 2) 8349. | 417x104 Œ 
hr gal hr 


For crossflow, the maximum available area for flow can be calculated from 


S.= ep (1-2) = 0167 ft? 
p 


using P = 6 inch, D, = 20 inch, Dg = 1 inch (see Appendix A.5-2), and p = 1.25 inch. The mass 
velocity for cross-flow is 


n aisi = 
Se ft^ - hr 


Let's assume that the fraction of tubes which pass through the baffle window equals the fraction 
of the shell cross-section represented by the baffle window. In other words: 


Np = fpx158 = 0.1955x158 = 30.8 


where 158 was given as the total number of tubes in the shell and fp = 0.1955 applies for 25% 
baffle windows (see (15)land Appendix on page bo]. Now everything in (14) lis known: 


2 2 
D D 
Sp - if 22 | Ji - 0258 ft? 


The mass velocity through the baffle window is 


Gy = 2162 x 10? es 
Sp ft^ - hr 
and the effective mass velocity is 
G, = JGG, = 2.01 x105 = 
ft^ - hr 


Before we can estimate ho from we need some fluid properties. Since we have no idea 
what the wall temperature is, let's just get a rough estimate of hy by setting wuy = 1. The other 
properties (besides wọ) are evaluated at the average temperature of the shell-side fluid [T = 
0.5(80°F + 122°F) = 101°F (38.3°C)]. Using Appendix A.2: 


= 0.681 cP, &- 0.363 BTU/ft-hr-°F, p= 62.0 lb/ft3, cp = 0.999 BTU/Ib-?F 


(13)]yields ho = 363 BTU/ft2-hr-°F 
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Once we have ho and h; we can estimate Tọ and iterate until the area required no longer changes. 
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Chapter 5. Evaporation 


evaporation -- concentrating a solution containing a nonvolatile solute by boiling away the 
solvent (usually water) 


Examples: vapor 


e production of orange juice concentrate 
e production of concentrated H5SO, feed 


in heat 
e production of distilled water 
Usually the product is a concentrated solution (called the 
“liquor”) and the vapor is condensed and discarded. 
liquor 


liquor - the concentrated solution 


An important exception is the production of distilled 
water from tap water. Here the condensed vapor is the product and the more concentrated 
solution of minerals or (in the case of desalination) salts is discarded. 


Evaporation is not the only unit operation that involves boiling. Some others that involve boiling 
include: 


e distillation - involves two or more volatile components (evaporation has just one) 

e drying - product is a solid (in evaporation, product is usually a more concentrated liquid 
solution) 

e crystallization - product is a slurry of crystals precipitated from solution (not a simple 
solution) 


Usually these differences are sufficient to require different equipment to perform the operation. 


EQUIPMENT FOR EVAPORATION 
There are several types of evaporators, but they all have the following parts: 


Part #1: HX (for adding the latent heat of vaporization to the liquid feed). Usually HX is a 
shell-and-tube with condensing steam as the hot fluid on the shell side (because steam 
is free of minerals which form scale). 


Part #2: vapor space or vapor head (larger chamber in which liquid, entrained in the vapor as 
droplets or foam, can be separated from the vapor, usually by impinging stream onto a 
plate.) 


You also need something to move the liquid through the heat exchanger. Sometimes a pump is 
used and sometimes it's just left to gravity, either in the form of a falling film on the outside of 
the HX tubes or the entrainment of liquid by rising bubbles. 
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Part #3: liquid mover 


e pump 
e gravity - falling film 


vapor 


- Vent 
steam 
(shell side) 


condensate 


feed 


liquor 


e gravity - entrainment by rising bubbles 


Long-Tube Vertical Evaporator 


When the heat exchanger tubes are oriented vertically, the bubbles produced by boiling rise 
and carry liquid with them. Tubes are typically 1 or 2 inches in diameter and 12 to 32 feet in 
length. The liquid entrained in the vapor coming out of the HX is richer in solute than that 
entering the HX. Vapor and fine liquid droplets are separated in the vapor head by placing a 
baffle plate in the path of the two-phase fluid coming out of the HXer. 


Notice the vent at the top of the shell-side of the HX. Usually a very small fraction of the steam 
flow is allowed to bleed out of the shell. This is done to allow air which is often mixed with the 
steam to escape. Otherwise non-condensable gases will accumulate in the shell, causing the 
pressure to rise and decreasing the steam flowrate. 


Long-tube vertical evaporators are especially effective in concentrating liquids that tend to foam. 
Foam is broken when the two-phase mixture hits the baffle plate. 
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Forced-Circulation Evaporator 


Many concentrated solutions are highly 
viscous. For example, a concentrated sugar 
solution can have a viscosity that is 106 
times larger than water. If the concentrated 
solution is highly viscous, the velocity of 
rising bubbles will be very small and the 
heat transfer coefficient will also suffer. 


vapor 


For very viscous liquids, a pump is often 
used to improve the circulation. Of course a 
pump also raises the pressure significantly 
which can prevent the liquid from boiling. 
In a forced-circulation evaporator, boiling 
usually occurs near the end of the HX where the pressure is lower or in the pipe leading from the 
HX to the vapor space. 


condensate liquor 


feed 


HX tubes are now oriented horizontally so they are easier to clean. Again a baffle plate is used to 
separate vapor from entrained liquid. 

EVAPORATOR PERFORMANCE 

Regardless of which type of evaporator we use, there are two standard measures of performance: 
capacity — pounds water evaporated per hour 


economy — lb water evaporated per lb steam used (usually somewhat less than unity) 


mr -m 

EXAMPLE: SIMPLE T,,H, 
EVAPORATOR PROBLEM 4inHg 

Suppose we are trying to gi | ODE ms, Iş, Hs 
concentrate a solution of orange juice e steam@ I5psig 
from 10wt% solids to 50wt% in a Hy 
long-tube vertical evaporator. Our He 
past experience has shown that a m 
reasonable heat transfer coef. is n 


Given: U= 500 Btu/h-ft2-°F 


To keep the boiling T low enough, 
we keep the vapor space under vacuum: 


p=4 inHg 
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If the feed is 55,000Ib/hr at 70°F, what area of heat exchanger is required and what flowrate of 
15 psig steam? 


Find: ms, A 

First we do a mass balance to determine the flowrates. A mass balance on the solids yields: 
10%(55,000 Ib/hr) = 50%(m) 

so m = 11,000 Ib/hr (liquor) 


is the flowrate of the liquor. The corresponding flowrate of the vapor must be (from an total 
mass balance): 


m p-m = 44,000 Ib/hr (vapor) 


The heat duty of the HXer is just the heat gained by the feed. In terms of enthalpies of the 
streams: 


q = vapor + liquor - feed 
q—(my-m)H,* MH- mH, (19) 


For high-molecular-weight substances like orange juice and other “organic colloids"* we can 
neglect boiling-point elevation and heat of dilution. Then the enthalpies of the feed and the 
liquor are essentially the properties of pure water under those conditions. The enthalpy of liquid 
water and water vapor can be read from the Steam Tables found in Appendices A.2-9 and A.2-10 
(pages 857-861).* 


Y “Organic colloids” are solutions of high molecular-weight polymers or dispersions of fine particles. Then the 
molar concentration of these solutions is negligibly small and their properties are close to that of pure water: see 
example on page 35. An example of an organic colloid is orange juice. 


* For consistency with later examples in which heat of dilution effects are important, we read all the enthapies from 
the steam tables, which have a reference state (where H is defined to be zero btu/Ib) corresponding to saturated liquid 
water at its freezing point. In Example 8.4-1 on p498, Geankoplis works a very similar problem in which all 
enthalpies are calculated using the conditions of the liquor as the reference state. Then (by definition) H in (19) is 
zero; the enthalpy of the feed 77, differs from that of the liquor just in sensible heat, which can be estimated as (8.4- 
9); and the enthalpy of the vapor H, differs from that of the liquor just in latent heat (i.e. H, is just the heat of 
vaporization at the T,p of the vapor space). See the discussion on “steam economy" on page 33. Try re-working 
Example 8.4-1 evaluating the enthalpies appearing in (8.4-7) using the steam tables (as above): you should get very 
nearly the same answer. 
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at 70 ?F:^ Hy = 38.09 btu/lbm 


But what is the temperature of the liquor? Well, since the liquor is in equilibrium with the vapor, 
their temperature and pressure should be the same. Since the vapor is just saturated steam at 4 
inHg, we look up its properties in the steam table: 


Py = 4 inHg = 1.97 psia 
from App. A.2-9 (p858): T= 125.1 °F 
Interpolating the tables,* we obtain 
at 125.1 °F: A, = 1115.7 btu/lb,, and H-93.1 btu/lb,, 


Substituting back into 


g= [44000 Al 1157 pe) + í 1,000 LAC pw) 2 (s 5,000 EAE 8.09 2) 
hr ib hr ib hr ib 


= 491x107 4010x107 —021x 107 = 480 x 107 = 


r 
Now, using the design equation, we can calculate the required area: 
q = UAAT,, (20) 


In a typical evaporator problem, steam is condensing on the shell side while water is evaporating 
on the tube side. When a phase change occurs with a single component (water in this case), the 
temperature stays at the boiling point. Thus the driving force AT is practically constant along the 
length of the heat exchanger. 


^ Strictly speaking, the enthalpies tabulated in A.2-9 are for saturated liquid, meaning that the pressure on the liquid 
equals the vapor pressure. Here the pressure on the feed is probably 1 atm (or 14.7 psia) rather than 0.3622 psia (the 
vapor pressure). Pressure has a very weak effect on the enthalpy of liquid water and the effect of pressure usually 
neglected. For example, the effect of a change on pressure equal to latm on the enthapy of liquid water can be 
estimated as 


S 
AH, - Ap V «(1 es [i58 = 0.044 — 
g 


which is equivalent to a change in temperature of only 0.044 °F. As long as we are dealing with much larger sensible 
temperature changes, we can neglect the effect of pressure on the enthalpy of liquid water. 


* In working the homework assignments, you will have frequent need to perform this interpolation. A Mathcad 
document has been prepared to assist you: http://www.andrew.cmu.edu/course/06-202/Steamtab.mcd This file reads 
the steam table from http://www.andrew.cmu.edu/course/06-202/STEAMTAB.PRN. Links to both files can be 
found at our homepage. 
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However, some changes might still occur at either end 

of the heat exchanger. For example, steam might enter Steam enters 
superheated, meaning that it's temperature is above the superheated 4; 
boiling point. Until the temperature drops to the “3 A 
boiling point, condensation does not occur. Similarly, “\ Steam leaves 
so much heat might be removed from the steam, that it subcooled 
completely condenses and is then subcooled below the 
boiling point. Finally, the feed solution has to be 
heated up to its boiling point before evaporation 
occurs." The temperature profile at right illustrates the 
most general situation. 


a 


N 


Feed solution 
enters subcooled 


Temperature 


Position in Hxer 


For a preliminary design (which is all we will do in 
this course), it is customary to approximate the “true mean” MDT as: 


AT, 


m IIT 
where T, — T of sat'd steam at p in steam chest 
T= Tp, of conc'd liquor at p in vapor space 
Knowing the pressure of the steam, we can find its temperature in the steam tables: 


at 15 psig — 29.7 psia: T, = 250°F 


The heat of vaporization à, is the difference between the enthalpy of saturated vapor H, and the 
enthalpy of the condensate H, leaving (assumed here to be saturated liquid): 


ùs = H,— H,- 1164.2 btu/Ib - 218.6 btu/Ib = 945.6 btu/Ib 


The driving force for heat transfer is the difference in temperature of saturated steam and 
saturated vapor: 


AT yy © T-T = 250-125 = 125?F 


Now everything in (20)lis known except for A: 


477x107 BU = [so 
h hee = 


Y The heat required to preheat the feed from 70°F (the feed conditions) to 125°F (the boiling point at the pressure in 
the vapor space) is 


55000? ag 38,09 bt" -53g 10! BRE 
hr Ib Ib hr 


This is only 6.3% of the heat duty g. We are willing to tolerate this size of error in our preliminary design. 
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or A=768 ft? 


The amount of steam required is just the heat duty divided by the latent heat of the steam: 


480% 107 "n A 
ri, = A = Bat = 50,8007 
so 956 n 


Recall that our two main performance indicators were "capacity" and “economy”. Capacity is 
just the rate of evaporation, which will be the flowrate of the vapor. For this example, 


capacity = my —m = 44,000 
r 


The economy is the mass of vapor produced per mass of steam consumed: 


T 44,000” 

economy = —— = hr = 0.87 
$ 50,800 — 
hr 


Steam Economy 

Why is the economy less than one? The answer becomes clear, if (instead of looking up the 
enthalpies of water and its vapor in a table) we estimate them as sensible and latent heat. First, 
let's select the liquor as the reference state for the calculations of the enthalpies. Then 


choose liquor as reference: H=0 


The vapor is at the same temperature as the liquor so there are not sensible heat changes, but the 
vapor has a much higher enthalpy owing to its heat of vaporization: 


H,7X 


Since the feed is also a liquid (like the liquor), its enthalpy differs from the liquor only in sensible 
heat: 


A= cy(Ty- T) 


Substituting these enthalpies in 


(mp m) — pe, (Ty - T) = rhs 
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mr-m m 
or f A -L e, (T T,) 
ms 2, Ms ` 
<l 


Generally the heat of vaporization of water decreases with increases in temperature. Since the 
steam must be hotter than the water being boiled, this means that A, < à. So the first term in the 
equation for economy is already less than one. The second term decreases it further; this term 
represents the energy required to preheat the feed up to the temperature of the liquor where 
boiling occurs. So a pound of steam boils less than one pound of water because some of the 
energy of the steam is required to preheat the feed. However, even if the feed were at the same 
temperature as the liquor, a pound of steam would boil less than one pound of water because the 
latent heat of the hot steam is less than the latent heat of the water being boiled. 


OVERVIEW OF EVAPORATOR DESIGN 


The design equation for evaporators can be summarized as 


heat duty: (mp -m)H, -mp Hy + mH =q = (H; — He) (21) 
————— > liquor : 
vapor feed steam 
HXer area: q=UAAT,, 


For preliminary designs, we can approximate 
AT, = const. = T, - T 


where T is the boiling point of the liquor at the pressure in the vapor space. This is approximate 
when there is significant BPE because 7 will vary along the length of the heat exchanger tube as 
the concentration of solute also changes. 


e Evaluation of enthalpies with negligible BPE and negligible heat of dilution (example on 
page p9 Jand also Example 8.4-1 on p498). Assuming the solvent is water, we can read the 
enthalpies of all three streams from the steam tables. 


read T, T,, H, Hy and H, from steam tables 


e Evaluation of enthalpies with complications (Example 8.4-3 on p501): 


read H, Hy from special table or chart 
read T from Dühring chart; T,, H, from steam tables* 


* Since enthalpies are being taken from more than one source, we should take care that the reference state is the 
same for all tables and charts used in a given problem. The steam tables in Geankoplis use saturated liquid water at 
32°F and 1 atm as the reference state. The enthalpy chart for NaOH solutions in Geankoplis also use 0 wt% NaOH 
at 32?F and 1 atm as the reference state. 
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COMPLICATIONS 


Boiling-Point Elevation 


Some complications arise in other evaporator problems because the enthalpy and partial 
pressure of a liquid more generally depends on its concentration as well as its temperature. 
Whenever a solute is added to water, the partial pressure of the water is reduced. Recall Raoult's 
law from thermodynamics. 


Roault's law: px. T ) =x; pP(T) 
where P; = partial pressure of component i 
x; = mole fraction component i 
pf = vapor pressure of pure i at some T 
Note that adding any solute with water causes x,,<1 so that 
Pw < Py? 


In general, this implies that I must heat the solution up to a higher temperature (i.e. higher than 
the b.p. of pure solvent) before it will boil. This increase in boiling point due to addition of 
solute is called 


boiling-point elevation (BPE) - increase in b.p. due to addition of solute 


This effect can be quite significant for concentrated solutions of inorganic salts. For example, for 
at a pressure of 1 atm, the boiling point of pure water is 


1 atm, 0% NaOH in water: Typ = 212°F 


If the pressure is other than 1 atm, the temperature can be read from the steam tables. Given this 
boiling point for pure water, we can look up the boiling point of solutions of NaOH from a 
Dühring chart (see p500 of Geonkoplis): 


50wt% NaOH in water: Typ = 283°F at latm 


Note that the boiling point of the solution is significantly greater than the boiling point of pure 
solvent. BPE is a colligative property, which means 


colligative property -- depends on molar (not mass) concentration 


Example: Use Raoult’s law to estimate the partial pressure of water at 100°C of a) a 50wt% 
NaOH solution and b) a 50wt% suspension of “organic colloids” (use MWayerage = 1000 g/mol). 
Estimate the BPE for each. 
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Solution: To calculate the mole fraction of water in the solution, take 1 gram of solution as a 
basis. This 1 gram mixture consists of 0.5 g of water (MW = 18) and 0.5 g of NaOH (MW = 40). 
We calculate the corresponding number of moles of each and the mole fraction: 


NaOH (MW-40) Organic Colloids 
(MW-1000) 


0.5g 05g 
Jn Hn 
Xw 0.5g " 0.5g RUARO 0.5g " 0.5g uis 
18-5. 49-8. 18-5. 1000-8. 
mol mol mol mol 


At 100°C (the b.p. of pure water), the vapor pressure of pure water is 1 atm (= p,,?). Using 
Raoult's law, the partial pressure is estimated as 


Pw 0.690(1 atm) = 0.690 atm | 0.982(1 atm) = 0.982 atm | 


To cause the solution to boil, the partial pressure must be raised to 1 atm. For the same mole 
fraction of NaOH, this means raising the vapor pressure to 


1 atm 
0.982 


] atm 
Pw? 0.690 


= 1.450 atm | = 1.018 atm | 


Using the steam tables (Appendix A.2-9), we find that the temperature has to be 


T | 110.0°C 100.8°C | 
Thus using Raoult’s law, the BPE is estimated to be 


BPE from 
Raoult’s law 


110.0°C - 100°C = 10.0°C | 100.8?C - 100°C = 0.8°C 


Of course, this is only an estimate. Raoult’s law is only asymptotically correct as x,,—1. The 
actual BPE is 39°C. 


BPE 283°F - 212°F = 71°F = = 
determined 39°C 
from Fig. 
8.4-2 


Although BPE is significant for a low-molecular-weight solute like NaOH, the same 50wt9 6 
concentration of an "organic colloid" (e.g. orange juice) may have a negligible BPE because the 
molecular weight is much higher. Thus the same wt% gives a much smaller molar concentration 
of a high-molecular-weight solute and therefore a much smaller BPE. 
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Heat of Dilution 


If you dilute concentrated sulfuric acid by adding water, large quantities of heat can be 
released. If the mixing is adiabatic and both fluids are initially at room temperature, the final 
temperature after mixing can by high enough to causing boiling. You were probably warned to 
use great caution when diluting sulfuric acid in chemistry lab. The problem is that the enthalpy 
of a solution depends not only on its temperature, but also on its concentration: 


H= K(T,c) 


which means you might need more than just a heat capacity to calculate sensible heat changes. 
Fig. 8.4-3 gives the enthalpy-concentration diagram for NaOH solutions. 


Example: Suppose we dilute 1-lb of 50wt% NaOH solution at 70°F with 1-lb of pure water, also 
at 70°F. What is the final temperature if the mixing is adiabatic? 


Solution: read the enthapies of the two solutions from 8.4-3: 


H(5096, 70°F) = 123 BTU/Ib 
H(0%, 70°F) = 40 BTU/Ib 


Multiplying each by the corresponding weight and adding, we obtain the final enthalpy of the 
mixture (“adiabatic” just means that no energy is lost during mixing): 


(123 BTU/Ib)(1 Ib) + (40 BTU/Ib)(1 Ib) = 163 BTU 
Dividing by the mass of the final mixture, we obtain the enthalpy per pound: 


163 BTU = sí 5 BIU = H(2595,T - ?) 
2 lb Ib 


This represents the specific enthalpy of the final mixture, which we also know is 25wt%. Again 
using Fig. 8.4-3, to get this H, we need to choose the temperature as 


T= 123°F 
Thus the final solution is warmer than the two initial solutions, due to the “heat of dilution.” 


Question: when diluting sulfuric acid (which also exhibits a significant heat of dilution), why is 
is wiser to add (concentrated) acid to water, rather than water to (concentrated) acid? 


We will now look at two methods to improve steam economy of an evaporator: 


1. multiple “effects” 
2. vapor recompression 
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MULTIPLE ÉFFECTS 


In the simple evaporators we have dealt with so far, one 
of the product streams is water vapor — also known as | N 
steam. Why can't this vapor be used in place of steam? The | rv | 
answer is simple: if the vapor replaced the steam, there V 
would be no driving force for heat transfer because the 
"steam" and the vapor would have the same temperature. To 
get a driving force in the right direction, we have to use 
steam with a higher pressure. One method of making use of 
the latent heat of the vapor is the compress it. We will look 
at that option below, but first let's examine another option. 


One way to improve economy is to divide the heat duty 
among several evaporators which are piped together in such a way that the vapor produced by 
one can be used to boil the liquid in another. This reduces the steam usage and improves 
economy. 


Ty Ds 


Effect #1 Effect #2 


You can have any number of evaporators connnected together in this way. Each evaporator in 
this sequence is called an effect, and the entire sequence is called a multiple-effect evaporator. 
In the figure above, we show a double-effect evaporator. In order to maintain a driving force for 
heat transfer in each effect, we must have 


Dpener 
This can be accomplished by having? — p; < py «p, 


This requires a vacuum pump or compressor: 


* Assuming we have saturated liquid and vapor in each effect, the pressure in any effect is just the vapor pressure at 
the temperature for that effect. Recall that vapor pressure increases with temperature. 
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po*py:* AI g= Ty -1n>0 

q2 = UpA5AT 5 (22) 
Of course, we still have to supply steam to the first effect 
AT ,=T,-T,>090 

qı = U\A\AT | (23) 


For simple designs where BPE and AH, are negligible and where the sensible heat needed to 
preheat the feed up to its boiling point is negligible, the same latent heat added to the vapor in the 
first effect is recovered by condensing in the second effect: 


ELT gg (24) 


Usually the pressure in the last effect and the pressure of the steam are known. Using the steam 
tables, we can convert these pressures into temperatures. The unknowns are the temperatures 
and pressure in the other effects. 


Example: given T, and 75, find Tj. Solution: (22), and [23)linto ields 
Uj Aj (T; — T1) = UjAxX(1 - T)) 


This represents one equation in one unknown (71), assuming that the U’s and A's are known. 


If instead, we have N effects, the analog of and [24)] yield N-1 equations in N-1 
unknowns. These equations are strongly coupled. A simpler approach is to try and solve for the 


AT;'s instead of the 77's. We will illustrate this by re-solving the last example. From|(22] and 
we have 


qı 
AT, = 21. 25 
1 UA (25) 
and AT, = —2 (26) 
Uy Ay 
Me — dau 
Addin and : AT; = —— + —— 27 
g 2. UA Ud Q7) 


Dividing [25)|by 


* The expression for AT which follows assumes that the vapor produced in the first effect condenses in the steam 
chest of the second effect at the same temperature at which it boiled. If BPE is significant, the vapor will condense 
at a lower temperature than it boiled (even 1f the pressure is the same). Why? 
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qi 
AT, _ Va M = iA (28) 
AT q2 ^X Ak 
2 A LM Uj4 /U54) 


Since XAT can be easily calculated from the info given in the problem statement: 


S AAT-ANAD-(0n-n)«(n-05)-21,-T, (29) 


we can easily calculate AT, from 


AT, = Lm (T, 7) 
JU, * pa 


The analog of e (29) for N effects are 


Les AT, =i for i-12,..N (30) 


qi 
m 
rcs F +f, T" + ANZ, 
U14 U24 Uy Ay 


Applying the analog of (24) ka; =q =... = qy = q), we can cancel out the q;, leaving: 
Y Naat xe + Fa Ay 


Now the sum of the driving forces can also be expressed as 


Y, AT = AT, + AD +... +AT y 


= (7, - 1) + (Gj -T )+..+Ty-2 - Ty) + (Ty - Ty) (32) 
= 1; -Ty 
(32)linto [31)] 
AT; = (T, - Ty) JA i-12,..N 


LU LO i DEA 
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In the simple case in which all the A's are equal and all the U 's are equal, this last result reduces 
to 


AT; = (r, -Ty) for i-12,.N 


In short, the total driving force is divided equally among each of the N effects. 


Performance 


Each effect results in the evaporation of water at a rate equal to q/A;. Summing the rate of 
evaporation in each effect, we obtain the total rate of evaporation: 


Nn 
ry ow 
ip 


Steam is used only in the first effect. The rate of steam consumption in the first effect is 


. qı 
m, = — 
S As 
The steam economy can then be calculated as 
di 
tiny Hj n= os =4N neve N 
m, q i3 |i 
À 


If we neglect changes in the latent heat with temperature and pressure, then all the A's are equal 
and the sum reduces to N. If changes in the latent heat with temperature are considered, the sum 
will be less than N. 


Summary: For the same capacity, steam requirements are significantly less for multiple-effect 
evaporators, than for single-effect evaporators: generally, the economy increases with the number 
of effects. The price that is paid for this increased steam economy is greater equipment costs: 
generally the cost of the equipment is proportional to the number of effects. 


For a single effect: economy < 1 
For N effects: economy < N 


The reason that the economy increased is that we need to supply steam only to the first effect. 
The source of heat for 2”4, 3rd, ... Nth effects is supplied by condensing the vapor produced in the 
previous effect. 
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VAPOR RECOMPRESSION 


Multiple effects recover some of the latent heat of 
the overhead vapor by using it in place of steam to boil 
the liquor in downstream effects. Is there any way we 
can use the vapor in place of steam in a single-effect 
evaporator? The problem with using the vapor in place feed 
of steam is that the temperature of the vapor is the 
same as the temperature of the liquor: the liquor and 
vapor leaving the same effect are already in thermal 
equilibrium with each other. In other words, there is no 


compressor 


condensate 


driving force for heat transfer. liquor 

without compression: AT=T,-T=0 

One way to increase the temperature of the 

vapor is to compress it. Fig. 1 shows how 

adiabatic compression Ware ind vapor 

(treated as an ideal gas) increases its qc 

temperature.* Note that the temperature Ei joo 

rises faster with increases in pressure than £ 

the boiling point increases. Thus adiabatic E 

compression of a saturated vapor nearly E 

always produces a superheated vapor. "= 

When heat is removed from this 

compressed vapor, it will eventually 

condense at its boiling point for that P 1 308 

pressure. Pressure (psia) 

Increasing the pressure of the vapor also Fig. 1 

increases its boiling point, so that as the Colored curves show how temperature 

vapor condenses, it stays at a high increases with adiabatic compression, 

temperature T, so we always have some starting with saturated vapor at 1 psia 

driving force. (green), 14.7 psia (red) or 400 psia 
(yellow). For comparison, the boiling 

With compression: AT- T, -T»0 point is shown by the black curve. 


Unfortunately, raising the temperature and 
pressure of the vapor actually lowers its latent heat per pound 


As <A (because T, > T) 


* Temperatures calculated from (7.30) in Smith & van Ness (4" ed., p237). 
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so the vapor no longer has enough latent heat to boil the same mass of liquid. Thus make-up 
steam is required, although generally much less steam is needed than if the vapor’s latent heat is 
not recycled. 


The reduction in steam is not without cost. We still have to pay the utilities (electricity or 
steam turbine) to operate the compressor. Generally, compressors are expensive to buy, 
maintain, and operate; so you have to do some economic calculations to see if using a 
compressor is economical. At the usual conditions of an evaporator, the vapor can be treated as 
an ideal gas. Then the work required for adiabatic compression can be estimated from * 


i 
PinY Pou d Btu 
pr | - 1| [=] — (33) 


(Y —1)Pin |\ Pin 


where “out” and “in” refer to the outlet and inlet of the compressor and where y is the heat 
capacity ratio, which is 


— 
for ideal gas 
for water vapor using 


c = 0.45 Btu/Ib-^F, 
R = ideal gas constant = 1.986 Btu/Ibmol-R = 1.986 cal/mol-K 
M = mol. wet. = 18 Ib/Ibmol = 18 g/mol 


Now gives the work needed to compress a unit mass of gas, under ideal conditions (no 
irreversible losses of energy). To get the power requirements, we must multiply by the mass flow 
rate and divide by the mechanical efficiency (n) of the compressor: 


ui i = [=] kW 
r 


P= 
n 


By compressing the vapor to raise its temperature, we are lowering its latent heat of vaporization. 
This means that, even if we completely condense all the vapor, we still won’t have enough latent 
heat to boil the liquor at its lower temperature. The difference must come from make-up steam. 
To estimate the make-up steam requirements, let’s return to the design equation (19)]for the total 
heat duty: 


* See Smith & van Ness, 4" ed., footnote on page 238. 
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q-(my —m)H, — Hy + mH 
— liquor 
vapor feed 


In the simple case where BPE and A77; are negligible (i.e. for inorganic colloids), Hy and H 
differ from the enthalpy of saturated liquid only in sensible heat. To the extent that sensible heat 
is negligible,* we can substitute the enthalpy of saturated liquid for Hyand H: 


heat to boil vapor: qm (m p= iJ (34) 


This heat duty is partially offset by the latent heat recovered by condensing the compressed 
vapor. The remainder comes from make up steam: 


q= (my-m)], + mj, (35) 
o ( —— 
compressed vapor make-up steam 


where we are assuming we are compressing the vapor to the same pressure as the make-up steam. 
Eliminating q between BA land we can solve for the steam requirements: 


m (iy (t 


S 


Notice that any compression allows us to recover all the latent heat of the vapor, thus 
significantly reducing steam consumption without compression. However greater compression 
actually increases steam consumption (because A, decreases with pressure). Rearranging this 
expression, we can calculate the economy: 


* Sensible heat is often negligible compared to latent heat. To see why, note that the latent heat of vaporization of 
water at latm, 212?F 1s 970.3 btu/Ib, where 1 btu is the amount of sensible heat required to raise the temperature of 
one Ib of water 1?F. This means that, for sensible heat to equal latent heat, I need to raise the temperature of water 
970?F. A temperature rise of only 30°F, say, represents a sensible heat change of 30/970 or 3% of the latent heat of 
water. This is an estimate of the relative error made in this estimate of m, . 
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which is now very much greater than unity. However, this way of computing economy is unfair 
because the reduction in steam consumption is at the cost of operating the compressor. In 
computing the “economy” for an evaporator with vapor recompression, it is customary to add an 
equivalent steam rate for the compressor to the make-up steam rate before dividing into the 
capacity. 


e If the plant uses a steam-powered turbine or if the plant uses steam generators to make its 
own electricity, then the apparent steam usage is 


UT " aid P [=] btu/hr 


Ib 
] 
Any  —btu/Ib 


rm 


where P is the electrical power supplied to the compressor, A, is the latent heat of the steam 


used (the prime is added because the pressure of the steam used here might be different from 
that used in the evaporator) and n, is the efficiency of electricity generation from steam 
(typically 0.35). 


e On the other hand, if electricity is purchased from a utility to drive the compressor, then it 
makes sense to convert the electricity usage into an equivalent amount of steam (on a 


= ms P X —— 


So of stm 


(ris Vp 


If P is in kW then the last term yields lb/hr of steam. 


Regardless of which method of converting the power usage is appropriate, the economy is 
computed as: 


my-m 
economy = l f 5 $300,000 - 
m 
S, app _ $250,000 
5 $200,000 —Hxer area 
. . 2 . s 
In Hwk #3, Prob. 5, we will consider the trade off  & $150,000 — electricity 
of operating costs (direct steam usage plus 4% Sandie Ed 
electricity for compressor) and capital costs (heat © 
exchanger area) as a function of the output $20,909 
pressure of the compressor. The results are LUNES: ERREUR ESSE: 
: N : 15 20 25 30 
summarized in the graph at right, which shows 
" n i Pressure of compressed vapor 
that an optimum pressure exists at which total (psia) 


costs are minimized. Notice that both power 
requirements for the compressor as well as steam 
usage increase with pressure. Only when the cost 
of the heat exchangers is considered does an optimum exit. 
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Chapter 6. Vapor-Liquid Equilibrium 
Reference: Smith & Van Ness, 4th Ed., Sect. 10.5 


Today we switch from heat-transfer operations (heat exchange, evaporation, condensation) to 
mass-transfer operations (distillation, extraction, absorption, humidification). Instead of sizing 
equipment to exchange heat between two phases we will size equipment to exchange mass 
between two phases. 


NO. OF PHASES PRESENT 
Suppose I have an apparatus like that shown in the figure 1 atm 
at right in which I can place an arbitrary mixture of 


components and control the temperature and pressure. How 
many phases will be present at equilibrium? 


One Component 


For example, suppose I put just water in the piston at 1 
atm and room temperature (say 70 °F) Just one component — AG 
no air! Is the water liquid, solid, or vapor? We know from 
our own experience that water is a liquid near room 
temperature. Thus we have just one phase present at equilibrium.* 


~ piston 


If I repeated this experiment at 250 °F, we would have the water present as all vapor — again 
I have just one phase present. Only if the temperature corresponds to the boiling point (212 °F) 
will I have two phases present. If I repeated the experiment at some other pressure there will be 
only one temperature at which two phases are present. 


The relationship between temperature and 
pressure for which two phases co-exist at equilibrium 
is called the vapor pressure curve. This diagram 
summarizes all the vapor-liquid phase behavior for a 
one-component system. 

1 atm 


Two Components 


212°F T 
Now suppose I put two components in the pot: 


let's say a mixture of n-heptane (C7) and n-octane 


* If air were present as a second component, we would also have a gas phase. 
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(Cg). which are important components in gasoline. The vapor-liquid phase behavior of two- 
component systems is a little more complicated, because I have an additional variable — 
composition — as well as T and P. Instead of a single temperature at which both liquid and 
vapor can co-exist, there is a range of temperatures. 


At a given pressure, we can summarize the phase 
behavior by means of a 7xy diagram like that shown at right. 
To understand the meaning of this diagram, consider what 
happens if we charge the piston with a mixture of heptane 
and octane having a particular composition. If the 
temperature is low enough, we will have just one phase: 
liquid. 


As we heat this mixture to the temperature corresponding 
to the lower curve, we would begin to see bubbles of vapor mole fraction C7 
formed in the liquid as more heat is added. This temperature 
at which the bubbles first form is called the 


bubble point — when heating a subcooled liquid, the T 


. tie line — 
temperature at which the first bubble forms OU 


The composition of this first bubble can also be read from 
the Txy diagram. In general, the vapor will be richer than 
the liquid in one component and leaner in the other. This 
difference in composition of the vapor and liquid can serve 
as the basis for separating two or more components of a 
mixture. By completely condensing the vapor, I produce a 
liquid which is richer in one component than the original 
liquid. Thus I have achieved a degree of separation of the 
two components. 


Thub 


For a binary system having two phases, we denote the composition of the liquid and vapor using 
x and y: 


x = mole fraction of more volatile component* in liquid phase 
y = mole fraction of more volatile component in vapor phase 


Keep in mind that liquid and vapor must have the same temperature at equilibrium. The line 
connecting the two points representing vapor and liquid is called a tie line. It's horizontal on a 
Txy diagram. 


* By “more volatile component" we mean that component which has either 1) the higher vapor pressure at a given 
temperature, or 2) the lower boiling point for a given pressure. Example: for a mixture of C} (heptane) and Cg 
(octane), heptane has the lower normal boiling point and so heptane is the more volatile component. 
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As we continue to heat the mixture, the liquid becomes less 
rich in the more volatile component since the more of the 
more volatile component is leaving for the vapor phase. The | — tie line 
compositions of the liquid and vapor move along the heavy 
solid line shown in the figure at right. The highest 
temperature at which liquid remains is called the 


dew point — when cooling a superheated vapor, the 
temperature at which the first drop of dew 
forms 


Eventually, as we continue to heat, all the mixture ends up in 
the vapor. * 


IDEAL SOLUTIONS OBEY RAOULT'S LAW 
Reference: Smith & Van Ness, Example 10.1 (p305) 
vapor Yi 


Now let’s try to predict the Txy diagram for the simplest case: an ideal 
solution in equilibrium with an ideal gas. So let’s suppose we are at some 
temperature for which both phases co-exist at equilibrium. 


Let P;' = vapor pressure of pure i 
P = total pressure 
x;y; = mole fractions in liq or vapor at equil. 
Dj * yjP = partial pressure of component i 


Raoult's law: pi ^ xP; 


IDEAL BINARY MIXTURES 


What we have said so far about ideal solutions applies to any number of components. Let's now 
restrict our attention to systems with only two components, which we will call A and B. 


The mole fractions must sum to one for each phase, so we can express the mole fraction of B in 
terms of that for A: 


XB 1-x4 


* The vertical line in left margin indicates that this material was not covered during class. I do however expect you 
to read it and you will be responsible for this material on the exams. 
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Yg = l-y4 
Applying Raoult's law to each component: 
P4A=YaP=x4P4' (36) 
Pg = (l-y4)P = (1-x4)Pp' (37) 


Adding the partial pressures from (36) hnd [37)]we obtain the total pressure: 


P=py4t pg =x4P4'(T) + (1-x4)Pg(T) 


Solving for the mole fraction x4: x4 = (P-Pg)(P4-Pp) (38) 
from [36) y47 x4P47P (39) 


Keep in might that the vapor pressures in these equations depend on temperature. So if you 
specify a T and P, I can determine the vapor pressures and calculate the composition of the vapor 
and liquid which will be at equilibrium. 


Given: TP 
Look up: P4'(T) and Pg T) 
Calculate: X4,y4 from (38)]and 


Thus I can construct a complete phase diagram for the mixture using only vapor-pressure data for 
the pure components. 


Example. n-heptane and n-octane form a nearly ideal mixture. Use Raoult's law to construct a 
phase diagram at constant pressure (1 atm). 


Solution. We first look up the normal boiling points of the two pure components. 


component Typ (CC) denoted by 


n-heptane 98.4 A 
n-octane 125.6 B 


Y Illustration 9.1 from Treybal. A similar example is 10.1 from S&VN (p305). 
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We choose the "low boiler" as T(°C) P'4 P'p x y a 
component A since it has the larger (mmHg) (mmHg) 

vapor pressure over this range of 98.4 760 333 1.000 1.000 2.28 
temperatures. Next, we choose a 105 940 417 0.656 0.811 2.25 
temperature between the two boiling 110 1050 484 0.488 0.674 2.17 
points, look up vapor pressures of the 115 1200 561 0.311 0.492 2.14 
pure components and substitute into (B8)] 120 — 1350 650 0.157 0.279 2.08 
and(39)| The results are summarized in 125.6 1540 760 0.000 0.000 2.03 


the table at right and the graph below. 


Temperature (deg C) 


0 0.2 0.4 0.6 0.8 1 


mole fraction of heptane 
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XY AND TXY DIAGRAMS 


The Txy diagram is a convenient graphical method 
for summarizing vapor-liquid equilibrium (VLE) data at 
a particular pressure. For binary mixtures, we usually 
denote the mole fraction of the more volatile component 
as x in the liquid phase as as y in the vapor phase. For 
this particular choice of x and y, the Txy diagram will 
always slope downward as shown at right. 


From the Txy diagram, one can construct an “xy” 
diagram, which is the one we will use in distillation 
calculations. The xy diagram is shown at right. For xy 
diagrams which resemble is one, the relationship 
between x and y can be sometimes be represented by a 
single number: the relative volatility, which is defined as 


Q dS 
;2 
Kj 


where the K 's are called the distribution coefficient, 
which in turn are defined as: 


pa” 
Xi 


Ti bp,B 


Spring, 2001 


mole fraction of A 


Ti bp,A 


Since this term is also useful in mixtures having more than two components, we will use 
subscripts to denote components. Generally the most volatile component will have the largest 
distribution coefficient and will have its concentration increased (in the vapor compared to the 
liquid) by the largest percentage. The ratio of the distribution coefficients of two components is 


the relative volatility. 
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NONIDEAL BEHAVIOR 


J 


Temperature, °C 
Temperature ,?C 


ONE ATM 


e 0.2 0.4 0.6 0.8 1.0 


O 0.2 0.4 0.6 08 1.0 x, y"=mole fraction acetone 


X, y*7 mole fraction CS» 


(b) (b) 
1.0 is 
S ONE ATM 3 
o o 
S OB a 
= o 
o S 
c o 
2 c 
5 2 
B8 0.4 $04 
pe E 
o o i 
E o2 202 | ONE ATM 
* " d F : 
` A 
0 O | 
O 0.2 0.4 0.6 08 1.0 O 0.2 0.4 0.6 0.8 1.0 
x = mole fraction CS2 in liquid x = mole fraction acetone in liquid 


Fig. 2. Txy and xy diagrams for nonideal systems. On left is carbon disulfide in acetone at 
one atmosphere pressure. On right is acetone in chloroform. Taken from Treybal, 2" ed., 
pages 290 and 295. 


Examples Txy and xy diagrams for binary mixtures which display significant nonideal behavior 
are illustrated in In particular notice that the curve on the xy-diagram can cross the 45? 
line. This point of crossing is called an azeotrope: 


azeotrope - vapor and liquid have the exactly the same composition at equilibrium (1.e. y — x) 


The occurance of azeotropes can pose difficulty for distillation since producing a new phase by 
condensation or by evaporation produces no separation of components. 
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Chapter 7. Flash Distillation 


Now let's apply these concepts about vapor-liquid 
equilibria to distillation. The simplest device for continuous Hp 
distillation of some liquid feed is the “flash unit" This is 
just a HXer to add enthapy of the feed and a chamber to F 
allow the vapor and liquid to separate. The schematic at XF 
right is the same as that used for evaporation, however the Hp 
vapor now contains all of the components — not just the 
"solvent". 


Given the feed composition and the amount of heat Hp 
added to the feed, we would like to determine the 
composition and flowrates of the distillate and bottoms 
products. 
Given: F, xp and q 
Find: D, B, yp and xg 
The mole and enthalpy balances are: F=D+B (40) 
xgF — ypD + xpB (41) 
HFF + q = HpD + HgB (42) 


where F,D,B are total molar flowrates of each stream and x and y are mole fractions of the more 
volatile component. The remaining equation comes from requiring the vapor and liquid products 
to be at equilibrium: 


(xg, Vp) lies on equil.curve (43) 


represent four equations in the four unknowns. 


APPROXIMATE SOLUTION 


The solution to these four equations can be obtained in a number of ways. Let's start with the 
simplest method. In many cases, the molar heat of vaporization is nearly the same for both 
components. (For this reason, we tend to solve all distillation problems using mole fraction 
rather than mass fractions, and express enthalpies per unit mole rather than enthalpy per unit 
mass.) If you can assume that the molar heat of vaporization of the mixture is completely 
independent of its composition, then solution is quite easy: 


e Assume Hp-Hp = const, say AH, [=] Btu/mol (independent of yp) 
e neglect sensible heat compared to latent heat 
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Let's choose saturated liquid having the composition of the feed as our reference state for 
enthalpies. Then 


Hp = Hg-0 
and Hp = AH, 


Even if the liquid is not saturated liquid, its enthapy differs from that of saturated liquid only 
through sensible heat, which is neglected here. The first equality is the result of neglecting 
sensible heats 


Then (42) |reduces to: q =x AH, D (44) 
Now we can define and solve for the fraction of the feed that is vaporized: 


D 4 
F AH,F 


f= 


which can be considered known. We will now express the molar flowrates in terms of this f: 


Dividing m and rearranging: £ =1 2 =1-f (45) 

Dividing [41)]by F: xp 7 ypf * xg(1- f) (46) 
; -1 

Solving for yp: yp= E Xp ue (47) 


which is a straight line on the xy diagram. To 
plot this line, we need to know the slope and one 
point which lies on the line. 


| 7 
Z A 
4 | . : 
4 operating line 


AU 


To obtain one point, we can substitute an | 
a 


arbitrary value in for xp, say we substitute 


xp-xp. Then becomes 


P pin = 
f 


f 


Thus one point on the line given by (47) lis 


| phase equilibrium 


XB Xp x 


YD = XF 


one point: (xB, YD) = Gt, XF) 
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Knowing one point and the slope, we can draw the line. In effect, this line represents the 
relationship between yp and xp which is imposed by the mass and energy balances for a given q 
(or a given f) and for the given assumptions on calculating enthalpy. We call this line the 


operating line: relationship between yp and xg which is imposed by the mass and 
energy balances 


To find which point on this line is the solution for a given q, we need also to satisfy which 
represents phase equilibrium. This relationship is a curve on the xy diagram, known as the 


equilibrium curve: relationship between yp and xg which is imposed by vapor- 
liquid equilibria 


The intersection of operating line and the equilibrium curve satisfies all four of the original 
equations and yields the answer. 


Inverse Lever-Arm Rule on Txy Diagram 


There is also a simple graphical procedure for estimating the T 
fraction vaporized on a Txy diagram. Suppose you are given a 
feed composition and (instead of specifying the fraction evaporat- 
ed) the final temperature of the mixture at equilibrium. Then you 
can locate a point (xp, T) on the Txy diagram. A horizontal line 
through this point intersects the vapor and liquid curves at the 
composition of those phases: (xg, T) and (yp, T). 


As it turns out, the amount of vapor produced is inversely Xs Xr Vp 
proportional to the distance between the feed point (xp, 7) and the 

vapor point (yp, T). Similarly, the amount of liquid remaining is 

inversely proportional to the distance between the feed point (xp, T) and the liquid point (xp, T). 
More precisely, the ratio of vapor D to liquid B flowrates in an continuous flash unit is given by 


D. Xr —Xg 


= 48 
B yp-Xr i 
where D, B and F are proportional to the distances 
indicated on the Txy diagram above. Because of the 
inverse nature of this relationship, this is called the 
inverse lever-arm rule. As a useful analogy, think of 
the flowrates as being like masses being balanced on a 
beam, resting on a fulcrum (see sketch at right). The 
further the mass is from the pivot, the small the mass 
can be. At (mechanical) equilibruim: 
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l 
Mıh = Moh or Sla 


Notice the inverse relationship between the mass and the distance from the pivot. Now we will 


try to prove 


Proof: 
Solving kolor f f= E (49) 
F yp-x*xg 
uni t D — — 
Substituting into m 2 =] ee eee a a (50) 
F F YD-XB YD-XB 


Dividing [49)|by (50)] we obtain 


RIGOROUS SOLUTION 


Although the approximate solution works reasonable well for some mixtures, for others it 
might not be good enough. Then you will need a lot more thermodynamic data besides latent 
heat and heat capacity. You will need an Hxy 
diagram. 

ideal case 
An Hxy diagram is at least superficially enthalpy 
similar to a 7xy diagram in that it gives you per mole sat'd vapor 


some idea of the thermal condition of the A 

mixture as a function of composition, except " Au 

that instead of T on the y-axis, we have H. 5 mcm 
sat'd liquid Y 


At right we show the Hxy diagram which 
would apply if 1) the molar heats of 
vaporization of the two components are equal, 
and 2) sensible heat is negligible compared to 
the latent heat: two parallel, horizontal lines. T a4 4 
Below it is an Hxy diagram for a real mixture. 


Unlike a 7xy diagram, the saturated liquid and 

saturated vapor curves do not connect at the ends or enthalpy 
at any other composition. The vertical distance 
between these two curves represents the latent heat 

of vaporization at a particular composition, which hs 
is usually substantial for all compositions. Unlike 
temperature, enthalpy of liquid and vapor will be 

quite different, so the fie lines are not horizontal 


.— — tie line | 
/ sae: i 
J^ oes auxiliary 
line 


XB YD 


mole fraction A 
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like on a Txy diagram. In fact, tie lines tend to be more nearly vertical. 


As we shall soon see, we will need to determine the tie line through a particular point. No matter 
how many tie lines your graph has, it would be a real stroke of luck if one passes through the 
particular point specified. So we need some way to interpolate tie lines. This is accomplished by 
means of: 


auxiliary line -- a curve on an Hxy diagram for interpolating tie lines 


If we form a right triangle out of each tie line, then the locus of the points where the right angle is 
located is the auxiliary line. Of course, the ends of the auxiliary line coincides with the ends of 
the saturated vapor line. By reversing this procedure, we can locate more tie lines. 


Now let's see how this diagram can be use do solve our flash distillation problems. First, I will 
give you the receipe and then later I will show why 
the receipe works. 


receipe: - find tie line 


e locate the feed point (xp, Hp) thru feedpoint 


H `N 
e locate a second point q/F above the feed point a E OU" 
(xp, Hp-q/F) 
ROME H F+ a ————— 
e find the tie line which passes through this 
second point B| X— | 
e ends of the tie line are (xg, Hg) and (vp, Hp). Xp! YD 
e f can be determined with the helg of (49): XF 
_ZF_*¥B 
YD- *¥*B 


Proof: now let's show that this procedure really works. 


The main idea which we need to prove is that the point (xp, Hptq/F) [call this point “F”] lies on 
the tie line. We know that (xp, Hp) [call this point *B"] and (xp, Hp) [call this point *D"] lie on 
the tie line since the two streams leaving the flash unit must be in phase equilibrium with each 
other. Our approach will be show that the slope of the line drawn from B to F is the same as the 
slope of the tie line drawn from F to D. 


Substituting (40) into xp(D* B) =ypD + xgB 
Dividing by D and solve for B/D: SU (51) 


D Xr —Xp 
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D 

Divide 2b F: Hp tte] p H 

y Ftp aD B 


Using (40) |to eliminate D: 


Solving for B/F: = (52) 


= (53) 


L IL] 


(52) + (53): = z (54) 


Equating K51)]and we could show that: 


q q 
Ap -|Hr+— Age | 
p-(we+4) (ne £)-m 


JD-*F XF—XB 


Note that both sides of this equation have the form of rise/run on an Hxy diagram; thus each side 
represents the slope of a line drawn between two points on an Hxy diagram. These two slopes 
involve three distinct points: 


(xp. Hg); (xp. Hp-*q/F); and (yp, Hp) 


The equality of these two slopes means that these three points must lie on the same line. That 
line connecting B and D must be a tie line. The conclusion is that 


(xp, Hp q/F) lies on tie line 


This is the basis for the receipe. 
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Inverse Lever-Arm Rule on Hxy Diagram 


find tie line 
thru feedpoint 


(54)|provides a basis for an inverse lever-arm rule 
which is similar n E The total difference in 
enthalpy between s ted vapor and saturated liquid 


is Hp-Hg.* This difference can be partitioned into Hp ——————— i 
(Hp*qlF) - Hg and Hp - (Hp*q/F). Thus the red tie j 
line can be divided into two parts which are Hp + 4 
proportional to this difference in enthalpies. says F : 
that ratio of the lengths of these two parts of the tie Hp Sa : 
line equal the molar flowrate ratio B/D. Hp : 
*Bi YD 
BATCH (DIFFERENTIAL) DISTILLATION id 
The oldest device for distillation consisted of a condenser 


kettle into which a liquid mixture was placed. By heating 
the kettle either by fire or steam or electrically, the liquid 
is partially vaporized. After condensing the vapor, you 
have a product which is richer in one of the components. 


If both of the components are volatile (e.g. ethanol and 
water) it wouldn't make sense to vaporize all of the 
liquid in the kettle. A mass balance can tell you that, if 
you vaporize the entire charge, the cumulative condensate 
composition will be essentially the same as the initial 
charge. Oh, perhaps you have removed some nonvolatile 
minerals which might have been in the water, but you 
won't enrich the ethanol concentration in this way. And 
enriching the composition of the more volatile component is the objective of distillation. 


condensate 
collector 


steam 


The solution is to evaporate only partially, because the first vapors which come off will be the 
richest in the more volatile component. Vapors produced later will only dilute the condensate. 
This process is called batch distillation because we process one batch at a time rather than 
continously processing as with flash distillation. 


Let's now try to predict how the composition of the liquid and vapor change as we evaporate 
more and more of the charge. 


Let n4(t), np(t) = moles of component A or component B left in liq. 


n(t) = n4(t) + npg(t) = total moles of liq. left in kettle at time t 


* Hp-Hg is not the “heat of vaporization”, which is the heat to totally vaporize one mole of feed. On the Hxy 
diagram, heat of vaporization is the length of a vertical line between the two curves. 
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x (t) = = mole fraction of A in liq. remaining (55) 


y(t) = instantaneous mole fraction of vapor produced 
Suppose we evaporate a small amount of liquid at time ¢ over a time interval dt: 
n(t*dt) - n(t) = dn « 0 
-dn — total moles evaporated in time df 
A fraction of this is component A: 
-dn4 = y(-dn) (56) 
(55)]into -d(xn) = y(-dn) 


xdn + ndx = ydn 


ndx — (y-x)dn 
dn dx 
n y-x 


Integrating from £—0, at which time x,n have their initial values: 
x(0)=x9 and n(0)=n9 
to some later time, at which time x,n have values: 


x(t) and nð 


we obtain — 
no ý X0 y VE 
n dx 
Integrating: In x = J 
is x(t) nod 
X0 d 
Solving for n(t): n (t) = ng exp| - f = 
-x 


x(t)? 
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which is called Rayleigh's equation. This serves as the design equation 
for batch distillation. To use this equation, we must recall that the vapor y 
coming out of the kettle at any time is in equilibrium with the remaining 
liquid: 


at equilibrium: y — f(x) 


For each x between x(0) and x(t), we can look up the 


corresponding y(t) and compute the integrand in 
Rayleigh's equation l/(y-x). The integral is the area »—* b | 
Xo 


under the curve between these two limits. 


Example: A 50mol% mixture of heptane and octane -Il 
is placed in a batch and heated. Plot the composition 

of the remaining liquid as a function of the fraction of — > 
the initial charge which remains. At the pressure in x(t) x0 

the still, the vapor-liquid equilibrium can be 

reasonably represented by a constant relative 

volatility of a = 1.7 (for a discussion on when this might arise, see p05). 


Solution: For binary solutions, it's customary to define x and y to be the mole fraction of the 
more volatile component (heptane in this case); this choice yields relative volatilities greater than 
unity. 


We will have to guess values of x(t). Below we go through the procedure for one particular 
guess, say x=0.2. 


Step 1: Guess x. We expect the vapor to be richer in the more volatile component 


Step 2: Calculate the corresponding y. Recalling the definition of relative volatility (see page 


B1] 
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boss Ki _ Yil Xi " y/x 
* Ky ya (1-»)-x) 


Dropping the subscripts and solving for y: 


= o 
d (a-1)x+1 = 


This function is evaluated for a=1.7 and plotted in the 
figure at right. For x=0.2, this yields, y-0.298. 


Step 3: We evaluate n/ng from Rayleigh's equation: 


0.5 
n(t) 7 =f dx 15 
no 02^ * 


The integral is evaluated by assuming x's between 0.2 « 
x «0.5. The integrand is plotted in the figure at right. 
The area under the curve between the two dotted lines 
is 


[-] 
Um I NL E E 
UA 
EN 


0.5 
f dx 245 
027 7 


This represents the fraction of the initial charge which remains in the still as liquid. The fraction 
of the initial charge which has been distilled is just unity minus this 


no 


1 = 0.914 
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Step 4: Calculate the condensate concentration 


If all the condensate is collected in a container 
and mixed together, the average concentration of 
the condensate will be 


x _ moles of A depleted 
cond — total moles depleted 


— X00 HO 0528 
no n 


The graph at right shows what typically happens 
as the initial liquid is more and more distilled. 
Notice that the vapor is initially richer than the 
liquid. Since the leaving vapor is richer, the 
remaining liquid becomes leaner. Thus both x 
and y become leaner as distillation proceeds. 


0.8 


Spring, 2001 


0 0.5 1 


x(t) fraction of initial charge distilled = 
— yt) [n(0)-n(t)]/n(0) 

xbar condensate 

en x(0) 
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Chapter 8. Multistage Operations 

In both batch distillation and flash distillation, the separation of the two components is limited by 

equilibrium. For example, suppose I have a 50-50 mixture of hexane and octane. If I vaporize 


half of this mixture, I produce a vapor having: 


xp=0.5, f-D/F-0.5 > yp70.55, xpg-0.45 


pure product? If I vaporized a smaller fraction of the feed, 0.55 
that would increase the composition of the distillate, but the 
best I can do is: 


which we read off of an xy diagram. Suppose you want a 95% — 0.58 4 


xp=0.5, fx0 > yp=0.58, xg-0.5 N 
This helps a little but it's still a long way from 0.95 and you 05 
don't generate much product with such a small fraction : 
vaporized. mole fraction hexane 


One solution is to take the distillate product of the first flash 
and partially condense it. This will produce a liquid phase which will be leaner in the more 


y1,Di7 0.5 y2,D2 = 0.25 


x1,B4— 0.5 X5,B5- 0.25 


mole fraction hexane 


volatile component, leaving the remaining vapor richer. 


This could be repeated many times until the vapor composition reaches or exceeds the desired 
value. Although we have produced a stream of the desired composition, it turns out that we are 
"throwing out" alot of the desired component in the "waste" streams -- the bottoms. 


From the first flash unit alone, we are losing 25% of component A in the initial feed [(1-f)x, ~ 
0.25]. If this component A is valuable, we don't want to throw it away. 


e success at achieving high yp, but low recovery (recovery is % of A in feed which 
ends up on final product) 
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e requires n+] HX'ers 

Another problem is that each flash unit requires a heat exchanger and then we will need one 

additional heat exchanger to condense the product from the final flash unit. 

COUNTER-CURRENT CASCADE 


Both of these drawbacks can be removed by a slight modification: we recycle the bottoms of 
one flash unit into the feed of the previous flash unit. 
»2,T) 


yp, ys. T: 


Instead of n+1 HXers, we now need only two. To condense the vapor, we need to transfer some 
of its latent heat to the liquid recycled from upstream. For this to work, the temperature of the 
vapor must be greater than that of the liquid. Let's see: generally, the vapor will get richer in the 
more volatile component as we move left to right: 


V1 $y2 S. SYN 


Assuming that the vapor and liquid streams leaving each flash 
are in equilibrium, then the corresponding inequalities for 
temperature can be deduced from the typical Txy diagram at 
right: 

Ty ete iw Ty 


So the recycled liquid will be cooler than the vapor it's mixed 
with. Thus the liquid will absorb some of the heat of the vapor. 
This condenses some of the vapor and vaporizes some of the 
liquid. Now that we have shown that the heat exchangers 
might be eliminated, let’s return to a discussion of the countercurrent cascade. 


The product recovery for a countercurrent cascade is much better than for the series of flash 
units without recycle. First, we have only one waste stream (B) -- instead of one from each flash. 
Second, the composition of this waste stream will usually turn out to be leaner in the valuable 
component -- because we are enriching the feed to the first flash by recycling from upstream. 
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Note that we also recycled some of the product liquid to condense the vapor feed to the last stage. 
This recycling of the product is called 


reflux -- that part of the product stream recycled to the cascade 
reflux ratio -- L,/D (ratio of molar flowrates) 
stage -- one unit of a cascade (e.g. one flash unit) 


ideal stage -- a stage whose exit streams are at equilibrium 
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Chapter 9. Equipment for Tray Towers 


In our schematic of a countercurrent cascade, we had the streams moving horizontally. 
Generally, a pump would then be required between each stage to overcome friction. With vapor 
and liquid, the density difference is sufficiently large, that we can use gravity to move the fluids 
if we stack the stages vertically rather than horizontally. A vertical stack of stages is called a tray 


tower. 
ie nJ* 


At right we see a schematic of a typical fray tower 
(Treybal, p129). | Each stage or flash unit is 
represented by one tray. Five trays are evident in — 4, _ 
Liquid flows down from the tray above and 

flows across the tray where in comes into contact 
with vapor flowing upward from the tray below. 
After mixing, the liquid and vapor are separated, 


os Ligand in 


Cewrapart -.2 Ji! 


with the new liquid flow down to the tray below and 
the vapor flowing to the tray above. Through the 
act of mixing, the composition of the liquid and 
vapor streams has changed. 


Fig. 4. An expanded view of a single 
“bubble cap” 


The contacting of vapor and liquid on any given tray 
is facilitated by bubble caps, one of which is shown 


in more detail in| Fig. 4. Each tray contains many 
bubble caps arranged in a array covering the central 


region of the tray (see Fig. 5]. 


Eunale cnr 


"ES Scopcrt 


rirj— — 


Tray " 
hitene: -U 


1 
d Sidgstren v 
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= 
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Fig. 3. A typical tray tower. 
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In this particular design shown ih Fig.] 4, vapor 
penetrates the tray through a circular hole 2.47 


inches in diameter, upward through a 3-inch sO Ac C No C 
length of pipe before being diverted downward by S 


i opo Doo 


ES Q e an n 


The flow pattern under both normal and abnormal | O05! Ae (3 
operation is shown in Ei | c AO PS bm d 
' : M a n e» ‘ee 
O26 ea 
= Pi Qe =OS 
28686 
JOS C 
- 5O 1 m 


wy " 


Dey 


the cap and forced through the liquid outside the 
pipe. / 


Fig. 5. Top view of arrangement of bubble caps on 
a tray. 


Pita dps Gas aki 1g 
d "nr Wreriul aie Hi 
Rack -trapp nig EM opercrior cap skit 


Fig. 6. Flow patterns under normal and abnormal operating of a bubble cap 
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An alternative tray design is the sieve-plate tray (see 

i . Sieve-plate trays are much cheaper to fabricate 
ea to bubble-cap trays since the tray consists of 
a circular plate with holes drilled through for the 
passage of vapor. These holes are generally much 
smaller in diameter than the bubble caps. 


A disadvantage of sieve plates is that they are much 
more prone to weaping and flooding (conditions 
corresponding to poor contact of liquid and vapor). 


—-|-lna3:2n 1* 
apo 


"— Perforared 
siert 


Fig. 7. Side view of a “sieve plate” tray. 
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Chapter 10. Performance of a C-C Cascade 


ANALYSIS OF A C-C CASCADE 


At right is a schematic representation of a countercurrent cascade. It 
ignores some of the details of plumbing inside the stages and, of course, in 
reality there are no pipes to carry the liquid and vapor between the 
adjacent stages. But these are unimportant details as far as mass and 
energy balances are concerned. The schematic preserves the essential 
feature -- which is that the flow of liquid and vapor occurs 
countercurrently. 


Since we have so many streams to deal with, some convention for 
numbering the stages and labelling the streams will be helpful. McCabe, 
Smith & Harriott use the following conventions: 


tray numbers: increase in direction of liquid flow (12i2.N) 
stream labels: subscripts denote stream 


e Lj V; are molar flowrates of total (A+B) liquid and vapor leaving the 
ith stage. 


e La Va= V; are the flowrates at the top of the cascade while Ly, Vp are 
the flows at the bottom of the cascade ("b" for bottom)." 


e x; y; are the mole fraction of the more volatile component in the liquid and vapor leaving the 
ith stage. 


Suppose we know everything about the two input streams: their flowrate, composition, and 
enthalpy. Suppose we also have all the thermodynamic data: Hxy diagram and xy diagram. 


Given: L4X4,Vpy,yp; Hxy and xy diagrams 


We want to determine the flowrate, composition and enthalpy of the two output streams and of 
all the intermediate streams: 


Find: Ljixj,Vjy for i=1,...,.N (= AN unknowns) 


I haven't bothered to list the enthalpies or temperatures here, because if I know the mole fraction 
and pressure, I can look up the enthalpy and temperature of saturated liquid or vapor. 


"aW 


"In heat transfer, the subscripts "a" and "b" were used for entering and leaving streams, respectively. That 
convention is abandoned here. 
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"n 
I 


As our system, let's select the first 
requires: 


stages. An overall mole balance for this system 


total: Lat Vi = Lit Vos (57) 
while a mole balance on the more volatile component requires: 

component A: Lg, + Viyia = Lixit Vaya. (58) 
Finally, we can also perform an enthalpy balance: 

enthalpy: Laly a ViaH, ja = LH, + Valya (59) 


Do I have enough information to calculate the result? How many independent equations do I 
have? 


equations: (sT)kssW59)x N ( = 3N eqns) 
x; and y; at equil. (=N eqns) 


So, at least when I am dealing with ideal stages, I have as many equations as unknowns. 


Many distillation towers have N=20 stages or more, so I could easily be talking about 80 
equations -- a formidible problem. Fortunately, there are assumptions which can make this 
problem tractable so that it can be solved easily with pencil and paper. 


CONSTANT MOLAL OVERFLOW 


In our discussion of flash distillation, we noted that many pure 
components have the same heat of vaporization per mole. If you can 
assume that the molar heat of vaporization of the mixture is completely 
independent of its composition, then solution is greatly simplified. 


Assumptions: 


1) H,-H, for all i. 
2) neglect sensible heat compared to latent heat 


Let's choose saturated liquid having the composition of the feed as our 
reference state for enthalpies. Since the only difference between the 
various liquid streams is sensible heat, we take them all to be equal -- 
equal to zero: 


Aly g =H =9 for all i 


and Fy i+) = Ay g 5A 
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With these assumptions, the enthalpy balance about the top i stages 


LaHy aq * Vit Ay ig = Lili + Vaya (60) 

simplifies to: Vi4 = Va for alli (61) 
and the total mole balance about the top i stages 

Lat Vpn =Lj+ V, (62) 


Subtracting (16)|from LaF” Li for all i 


In other words, all of the vapor streams have the same molar flowrate within the cascade 
(although their composition and mass flowrate will change) and all the liquid streams have the 
same molar flowrate: 


VaV 5n... =Vy= Vp (= V, say) 


La Ly Ly ee LyeL, (=L, say) 


which is called the state of constant molal overflow. Since the input flowrates are known, we 
know all the flowrates. Thus we have just halved the number of unknowns from 4N to 2N. We 
can then drop the subscripts on the flowrates and forget the enthalpy balance (60) ]and the total 
mole balance This leaves just the component mole balance 


Laxa T Vis itl = Lx; + VaVa- 


L L 
which can be rearranged to: Jiu = y + b. = aa (63) 


or V4] = mx; + b 

which is a straight line on an xy diagram. This line is called the: 

operating line — y;,4 vs. x;; a graphical representation of the mass balance. 
Another important curve on an xy diagram is the: 

equilibrium curve — y; vs. xj; a graphical representation of VLE 


In the next section, we will see how these two curves can be used to determine the number of 
equilibrium stages required to achieve some desired degree of separation. 
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A RECTIFYING CASCADE 


In a typical distillation tower, the liquid fed to 
the top of the enriching cascade is a portion of the 
vapor product which has been condensed and 
refluxed back into the tower, as shown at right. 
The ratio of the amount refluxed to the amount not 
refluxed is called the reflux ratio 


2 
D 


R= 


(64) 


Since all of the vapor stream leaving the top of the 
cascade is condensed to form the distillate product 
and the liquid fed to the top, the three streams must 
have the same composition: 


XD — Xa — Ya (65) Fx,=0.5 


Then the operating line (63) |can be simplified to^ 


L 
Vido zs d >*D (66) 
zm D 
V 


The ratio of the flowrates can be expressed in terms of the reflux ratio R using 


L Lt _ RD _R (67) 
V L+D RD+D R«l 
D D | D _ 1 
V L«D RD+D R«l 
The operating line (66)]can be rewritten as 
R 1 
HQ x; + x 68 
Yi tt Ra (68) 


Thus the slope of the operating line (a plot of y;,, versus x;) is R/(R+1). 


^ Alternatively, we could "start over" and perform a component mole balance on the top i stages of the cascade 
(indicated by the blue line in the sketch). 
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Now, one point on the operating line is x; = y= xp. In other 
words, if we substitute xp for both x; A. (68) is 


automatically satisfied. xp 
Example. Suppose I have a feed stream containing 5096 
alcohol in water which I want to enrich to 95% alcohol. How A 
many stages do I need if I use a reflux ratio of unity and what 
will be the flowrate of this product? 
Given: xp, xp, F and R 
xp 


Find: Nand D 


Solution: The first step is to draw the operating line on the xy- 
diagram. One point is (xp,xp) or (0.95, 0.95). And the slope is 
R/(R +1) = 1/2 for R=1. Armed with one point and the slope, we could draw the line. But it is 


usually easier to construct a line if we know two points. A second convenient point corresponds 
to the y-intercept. Substituting x; = 0 into (68)lyields 


x XD 


R-«1 


xj = 0: Yi+1 = 


which can be conveniently calculated from what’s given. Thus 

XD 
(o zn.) is a second point on our operating line. Knowing the 9,1 
+ 


slope and one point, or two points, we can draw the operating 
line.” Once the operating line is plotted on the xy diagram 
(which is the equilibrium curve), we can determine the XD 
composition of other streams as follows. 


The two streams leaving a given ideal stage are at equilibrium, 

thus the compositions of the two streams leaving a given stage must lie on the equil. curve. 
Since we already know the composition of the vapor stream leaving stage 1 (y, = xp = 0.95), we 
can determine x, by moving horizontally from one point on the operating line (xp, xp) to one 
point on the equilibrium curve and then read off the value of x): 


(x11) on equil. curve: gives x 


* The point (xp,xp) and slope of 1/2 are not drawn to scale in the sketches. For that matter, neither is the 
equilibrium curve. However, given the equilibrium curve and operating line as drawn, the method of stepping off the 
number of stages is correct. 
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On the other hand, the operating line represents a component 4 
mole balance about an arbitrary number of stages measured from 5 PU 
; . : : Y2 
the top. In particular, for just the top stage, the operating line 
(66)]eives me a relationship between x, and y»: 


(x1,y2) on oper. line: gives y» 


Knowing x, I can determine y». In general, the compositions of 

the two streams between a given pair of adjacent stages must 
lie on the operating line. Now 
I can repeat the process: 


“TF 
y3 | Gf (x2,y2) on equil. curve: gives x» 


(x2,V3) on oper. line: gives y3 


X] Xp 


and so on. Continuing to 
step off stages in this manner 
until the vapor mole fraction 

X2 equals or drops below xp = 
yp. requires four steps: 


feed line for sat'd vapor 


answer: N=4 


Comment #1: Although the feed might be a liquid, it’s the 
vapor mole fraction which must be at or below the feed 
composition, because the actual feed to be bottom of the 
cascade is a vapor (see stream whose composition is labelled 
yp in the flowsheet on page 7b). The horizontal line y = xp 
(= y,) on the xy-diagram is called the feed line. 


Comment #2: After taking 4 steps (or any integer number), we 

don't end up exactly at the feed composition (except by some extreme stroke of luck). Instead 
we end up stepping below the feed composition. If the trays really are perfectly ideal, and the 
feed composition is 0.5, the distillate composition would turn out to be slightly richer than 0.95. 
Although any real column has an integer number of trays, it is customary to report the answer as 
a fraction; for example, to go from a feed of 0.50 to a distillate of 0.95, we require 3.7 ideal 
stages (or 3 + 0.7, where 0.7 is the fraction of the last step which corresponds to the feed 
composition) in the current problem. 


Comment #3: Suppose that our feed composition had been a little lower, say yp’. No matter how 
many steps we took, we could never get a vapor composition lower than that at the intersection 
of the operating and equilibrium curves. In that case, any solution for the specified reflux ratio is 
impossible. Instead the reflux ratio would have to be increased, which will rotate the operating 
line counter-clockwise about the specified point. This will lower the composition at the 
intersection point and make the separation possible. It turns out there exists a minimum value of 
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the reflux ratio which allows a particular distillate composition xp to be produced from a given 
feed composition xp. This called the minimum reflux ratio. 


Let's finish this problem by calculating the flowrates. Since the vapor flowrate leaving each 
tray is the same, and because we are creating the vapor to the first tray by completely vaporizing 
the feed, we can say that 


T 
ES 


The ratio of liquid to vapor flowrate inside the cascade is given in terms of the known reflux ratio 


by [67)] 


L R RN ns 1 


y R+1 2 2 2 


Since the bottoms stream B is one of the liquid streams 


palel? 
2 


Finally, the distillate flowrate can be calculated from the definition of the reflux ratio: recall 


A STRIPPING CASCADE 


Use of a single cascade with 3.7 stages allows us to convert 
a feed of 50mol% alcohol into one containing 95mol%. But a FEx;=0.5 
significant amount of the alcohol is being lost in the second a {hy Va 
product stream. It’s possible to strip the alcohol from this 
stream by feeding it to a second cascade. 


In a stripping cascade, liquid is fed into the top. Liquid out 
the bottom of the cascade is partially vaporized and the vapor is 
recycled into the bottom of the cascade. The vapor stream at 
the top represents the second product stream. 


Because the bottoms product is only partially vaporized to form 
the vapor stream which is recycled to the cascade, its 
composition is not the same as the bottoms product. In other 
words 


XN *YN+1 * Xp 


So no analogy can be made to (65). However we still have a 
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operating-line relationship like Or which is imposed by mole balances. From a balance 
over the entire stripping cascad h e volatile component, we know that 


Lx, Fxp = Bxg + Vya 
——— — 


in out 
so Vya - LXq = -Bxp 


Substituting this result into the operating line for a stripping 
cascade becomes^ 


L B 
Jia = y Xi y *B (69) 4, 
Like all the operating lines we've encountered so far, this one has Z 


a slope of L/V. One point on this operating line is (xg.xg) since 
substituting y;,, = x; = xp automatically satisfies (69) provided L 


- B = V which is always true since this represents a total mole xB 

balance. Another point is (xp,y1) if we happen to know or can 

somehow calculate yj. Then we have two points and can draw XB XF 
the operating line. Otherwise, we use one point and the slope. 


Example. Suppose I have a feed stream containing 50% alcohol 

in water which I want to strip the alcohol from to produce a stream containing only 5% alcohol. 
How many stages do I need if I boil half of the liquid leaving the bottom of the cascade and 
return it as vapor? 


Given: Xp, Xp and V/L 
Find: N 


Solution: Once again, vapor and liquid streams leaving the 

same stage are assumed to at equilibrium; thus (xj,y;) lies on 

the equilibrium curve. On the other hand, the composition of 
adjacent streams (xjyj,,) lies on an operating line for this 
cascade. One point on this operating line is x; = yj;4 = xp. If 

we “boil half the liquid leaving the bottom of the cascade and yy), 
return it as vapor" we know that V — L/2, so the slope of the 
operating line is: 


LIV=2 XB XN 


Knowing the slope and one point, we can immediately draw 


^ Alternatively, we could “start over” and perform a component mole balance on the bottom N-i stages of the 
cascade (indicated by the blue line in the sketch). 
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the operating line. The liquid product being drawn off at xg is in equilibrium with the vapor 
being sent back into the cascade; thus (xg, yy+1) lies on the equibrium curve. On the other hand, 
(xy, YN44) are adjacent streams in the cascade and so their compositions must lie on the operating 
line. 


Similarly, (xy, yy) lies on the equilibrium curve, while (xy.4, feed line for a sat'd liq 
yy) lies on the operating line. 


We continue to step off stages in this manner until we step past 
a liquid composition equal to the feed xp; in other words, we 
continue stepping until we cross the blue vertical feed line. In YN 
this case, three steps are required. But one of these steps is the 


reboiler; only two additional equilibrium trays are required: YN+1 
answer: N =2 plus reboiler 
Comment #1: Notice that the “feed line” for a stripping XB XN XN-1 


cascade is vertical, whereas the feed line for a rectifying 

cascade was horizontal (see page This is because the feed 

into the bottom of a rectifying cascade is a vapor, whereas the feed to a stripping cascade is a 
liquid. 


Comment #2: Suppose that our feed composition had been a little higher, say xp’. No matter how 
many steps we took, we could never get a liquid composition higher than that at the intersection 
of the operating and equilibrium curves. In that case, any solution for the specified L/V ratio is 
impossible. Instead the L/V ratio would have to be decreased, which will rotate the operating line 
clockwise about the specified point. This will raise the liquid composition at the intersection 
point and make the separation possible. It turns out there exists maximum Z/V, or a minimum 
vapor rate V which allows a particular bottoms composition xp to be produced from a given feed 
composition xp and liquid flowrate L. This minimum V/L is called the minimum boil-up ratio. 


Let's finish the problem by determining the flowrates. Once again, all the flowrates can be 
deduced from the equimolal overflow assumption and the reboil fraction. All the liquid 
flowrates correspond to the flowrate of the liquid fed to the top: 


L-F 
Since I’m reboiling half of this liquid to return to the cascade as vapor, all the vapor flowrates are 
V=L/2 = F/2 
The bottoms flowrate is what remains after half of the liquid to the reboiler is boiled: 


B=L-V=F/2 
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3:205 
STEAM DISTILLATION Lx4 yo, 
Instead of using a heat exchanger to partially boil the liquid 
stream leaving the bottom stage to produce the vapor entering the Lx, ! | by; 
bottom stage, we could directly inject steam as the vapor. This 2 
has the advantage that it avoids the cost of a heat exchanger and L, xy i Vy, 
you are getting direct contact between the hot and cold fluids. 
Of course, use of steam as the vapor is only useful if water is LXy, i ii A 
one of the two components being separated (e.g. alcohol/water) N 
and further, it must be the less volatile component. If instead you L, X $ Vy, =0 
are trying to separate heptane from octane, injecting steam would (live stean 
be produce a second liquid phase (since oil and water are not 
miscible). 


This produces the operating line shown at right. Note that 
one point on the operating line is (xg, 0) since the “vapor” 
entering the bottom is pure steam (water). Of course, part of the 1 
operating line lies below the 45° line (x=y). But there are no 
restrictions that y must be larger than x. As with the stripping 
cascade, we continue to step off stages between the operating line 
and the equilibrium curve until our liquid composition exceeds 


XF. 


Comment: Suppose that our feed composition had been a little 
higher, say xr’. No matter how many steps we took, we could 
never get a liquid composition higher than that at the intersection 
of the operating and equilibrium curves. In that case, any solution for the specified L/V ratio is 
impossible. Instead the L/V ratio would have to be decreased, which will rotate the operating line 
clockwise about the specified point. This will raise the composition at the intersection point and 
make the separation possible. It turns out there exists maximum Z/V, or a minimum steam rate V 
which allows a particular bottoms composition xg to be produced from a given feed composition 
xp and liquid flowrate L. 
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Chapter 11. McCabe-Thiele Method 


A typical distillation tower will consist of at least two V 
cascades -- one to enrich the feed in the more volatile 
component and a second one to strip the more volatile 
component out of the bottoms stream before disposing of 


condenser 


it. enriching 
"TP cascade 
rectification -- enrichment of vapor by contact with liquid 


reflux 


With enough trays, we can usually make D as pure as we 

like (an exception occurs for an azeotropic mixture). 

However, the liquid stream at the bottom of the cascade H, 

might still contain a significant fraction of the more 

volatile component. This is where a second cascade is 

useful: stripping 
cascade 


feed tray 


stripping -- removal of the more volatile component from 
a liquid by contact with reboled bottoms (or direct 
injection of steam) 


e 


A typical problem in designing a distillation column is the 
following: 


Given: zp Xp,Xg.L/D, Hp 

Find: Nand N 

where N and N are the number of stages in enriching and stripping sections. Having two 
cascades is a little tricker than one, so let me first tell you how to solve the problem by giving 


you the receipe, and later we will rationalize some of the steps in this receipe. 


This receipe is called the McCabe-Thiele Method and it works when you have only two 
components and you can assume equal molal overflow. 


Step 1) Plot equilibrium data on an xy diagram. 
Step 2) Locate xp,zp, and xp on 45? line 
We use "zp" to denote the mole fraction ignoring which phase the feed is in: 


. moles of A in liquid + vapor 


FA = 
F — total moles in liquid + vapor 
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Step 3) Calculate q (a dimensionless variable representing the 
thermal condition of the feed): 


n Bue 
TUH y -Hr 
where Hp = enthalpy of feed 


Hy = H of sat'd vapor with y=zp 


Hy =H of sat'd liq. with x=zp 


Step 4) Plot q-line, which is a line given by the following 
equation (Q.L.): 


XD 
Of course, the name "g-line" is given because this line p41 
depends on g — the thermal condition of the feed. For 
example, suppose the feed is ... 


Xj Zi Xp 
superheated vapor: mole fraction hexane 
Hp>Hy > q0 > 0< i «1 L&V satdL 
sat'd vapor: 
Hpg-Hy > q-0 > E = 0 (horizontal line) 
2-phases mixture: Hr<Hp<Hy — O<q<l > aA <0 


In this case, q represents the fraction of the feed which is liquid. 


sat'd liquid: Hpg-Hyj > q=l > 2 i = oo (vertical line) 


1«—- «re 


q-— 


subcooled liquid: Hr<H, > el > 


As it turns out, the two operating lines for the two cascades intersect somewhere along this line. 
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one point: (zp.zr) 
slope: qK(q-1) 


Stop 5) Plot the operating line for the rectification section (R.O.L.): 


R XD 
y=—— x+ 
R+1 R+1 
one point: (Xp.xp) 
another point: (0,xp/(R+1)) 


Step 6) Plot the operating line for the stripping section (S.O.L): 


L Bx B 
LA SSS 
A V V 
one point: (Xg,Xg) 
another point: intersection of R.O.L. and Q.L. 


Step 7) Step off stages 

You could start at either end (x=y=xp or x=y=xp), although it is customary to start at the top. In 
stepping off the stages, you eventually have a choice as to which operating line to drop down to. 
LOCATION OF FEED PLATE 


At some point during the stepping process, you will have a choice of either operating line to 
step down to. Either operating line can be used, but once you switch operating lines, you cannot 
go back to the original line: switching operating lines corresponds to switching from one cascade 
to the other. The stage at which we switch operating lines is the 


feed plate: stage at which you switch operating lines 
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Either of the two stepping schemes above is OK, but one scheme leads to fewer total trays than 
the other. 


optimum feed plate: to get the fewest number of total stages, you switch operating lines when 
that allows you to take larger steps 


In other words, if you have a choice of which operating line to use, use that one which takes you 
furthest from the equilibrium curve. This choice gives the biggest step and fewest total stages. 


CONVERGENCE OF OPERATING LINE AND EQUILIBRIUM CURVE 


In some of the homework problems, the distance between the operating line and equilibrium 
curve becomes very small near the point of their intersection. When this happens, a very large 
number of steps might be required to make any progress on separation. 


Rule: an infinite number of steps is required to get to an intersection of the operating line and 
equilibrium curve. 


Proof: when the size of your steps are no longer large compared to the width of the pen or pencil 
you are using to draw them, you should "blow up" the drawing to reduce drawing errors. 


To demonstrate this graphically, let me present the results of stepping off stages in a particular 
example: 


Example: given the xy-diagram for ethanol and water at 1 atm pressure (has an azeotrope at x = y 
= 0.8), step off stages for the following conditions: 


Given: Zp = 0.25, xp = 0.7999, xp = 0.01, q = 0.3 


Find: Nand N for the optimum feed location 
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Enter values: 


0.8 
.d ici 
nel Murphree plate efficiency 


xprime 


x p=0.7999 Mol frac in distillate 
x p=0.25 Mol frac in feed 


x p=0.01 Mol frac in bottoms 


q=0.3 Thermal conditions of 
feed (liquid fraction) 


R=8.618 Reflux ratio (to R nin? 4-483 
obtain total reflux, 
use R20) 


After N steps, this is actual xB 


-3 
xB calo = 5-57-10 


M? 0.2 0.4 0.6 0.8 Nrect - 60 Feed tray 
—— equilibrium curve 
Qr y=x line N =67 Total number of 
—— Rectifying section oper. line stages 
—— Qline 
— — Stripping section oper. line 
— — stepping off stages 


Solution: This is a particularly difficult separation because the distillation composition xp = 
0.7999 is very close to the azeotrope. We anticipate a large number of steps. Available online is 
a Mathcad document [McCabe-Thiele (Pro).med | which steps off stages automatically using the 
optimum feed location. 


The results are shown in the figure at right. The program reports (not shown) the total number of 
steps and the feed location. In this figure, there are a total of 67 steps with the feed introduced on 
the 60" from the top. In the upper right corner, the steps are becoming so small as to be 
invisible. 


The below on the left is a blow-up of the region from 0.7 to 0.8. Even on this blow-up the steps 
are becoming invisible in the upper right corner. So the second graph below on the right is 
another blow-up: this from 0.79 to 0.80. Still the steps are becoming vanishingly small. 


The red (equilibrium) curve does not show up on the second graph below on the right. 
However, the equilibrium curve is just the locus of points from the upper-left corner of each step. 
Similarly, the lower right corner of each step is the operating line. Notice that in these blow-ups, 
both the equilibrium curve and the operating line very linear. 
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a 
+ 0.98 ES 
| 0.96 
0.794 
0.792 
| | l | nee | | l | 
0.72 0.74 0.76 0.78 0.79 0.792 0.794 0.796 0.798 


Over a narrow range, any continuous curve becomes virtually linear. When both the operating 
line and equilibrium curve are straight lines, an analytic expression exists for calculating the 
number of steps: 


which is called Kremser's equation. The subscripts a and b 
denote streams at the top and bottom of the cascade, respectively; the y* are vapor compositions 
in equilibrium with the liquid at the top and bottom of the cascade. 


Proof: the linear equilibrium curve can be expressed as 
y; = mx; + b (71) 


A mole balance about the top 7 stages of the cascade yields 
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L L 
JAS Sat Ja ~ y Xa 
Using (71) |to eliminate x;: 


L y;-b L 
i+ = a Ya "y" 


~ 


Now we define a new parameter, which represents the ratio of the slope 
of the operating line (L/V) to the slope of the equilibrium curve (m): 


In terms of A, our mole balance becomes 


[ ] Yi+1 = A(yj - b) + ya — Amx, = Ay; + ya — A(mx, +b) 
— 


Ya* 


From (71), we see that mx,+b is just the vapor concentration in equilibrium with liquid at 
concentration x,, which we denote at y,*. After minor rearrangement, the above equation 
becomes 


Yi+1 = AV; - Aya * * ya (72) 


This serves as a recursion formula, which we can applied repeatedly to obtain the concentration 
on successive stages, starting at the top. For i=1, recalling y, = yg, [72) yields 


i-l: Y2 = Ayy- Aya * ya = Ya (1+ A)— Aya * 
Ya 


y3 = Aya — Aya * ya = A| ya (1+ 4) - Aya * ] - Ava * y; 


= ys (1e 4* P)». t (a 4’) 


where the second equality above results from substituting y) from the previous step. Now we 
begin to see a pattern developing: 


m Yast 7 Ya (E dt 4 soe a") vut (n eoe a") (73) 


The sums represent partial sums of geometric series, whose sum is known analytically: 
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_ ,4"041 
bP ALA Sot A" S uc dig 
1-A 
Substituting this result into 
]- 4”+! ]- 4" 
Yn+1 — Ya I_A Ya*4A T 


For the bottom stage of the cascade, n=N and yy+1 = yp: 


= HEN CUN RN MEO eA 
Yb = Ya {=a Ya ld 


Multiplying through by 1-4 and collecting terms multiplied by A or by AN*1: 
N+1 
Jp 7 Ya t A(¥q *-y5) 2 4" (Ya *-Ya) (74) 


This can be simplified. Recall that A 1s defined as the ratio of the slopes of the operating line to 
the equilibrium curve. This ratio can be calculated from the concentration defined on the xy- 


diagram above: 
[95 B i/i 
"T slope of OL _ scs) s Ju Yb (75) 


 slopeofEC (ya *-yp p Ya * yp * 
(xa = xp) 


We can use this result to obtain y, - yp = A(y,*-yp*), which will now be substituted into 


A( yp *— Ya *)+ A( Ya * —yp) = ave (Ya * —ya) 
eh SS 
A( yp *- Ya * Ya *-Yb )- A(Yo * y») 


After dividing through by A and solving for AN: 


d 
AN 2 Jb —YVb or N= Ya Ya 
Ya *—Ya log A 


Taking the log of both sides, recalling that log(4N) = N logA, then solving for N we obtain the 
second expression above. To obtain the form of Kremser's equation given as|(70), we substitute 


(75) |for A, 


Comment #1: In the xy-diagram shown next to |(70)) the equilibrium curve lies above the 
operating line. In other words, y,* > y, and yp* > yp. To obtain N>0, you must have y, > yp. If 
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you inadvertently (but consistently) interchange the subscripts a and b, you will obtain negative 
values for N. Kremser's equation can also be applied to situations in which the operating line lies 
above the equilibrium curve; then to obtain N70, the inequalities above must be reversed. 


Comment Z2: If the operating line and equilibrium curve should cross (which is physically 
impossible), when evaluating |(70)| you will find yourself trying evaluate the log of a negative 
number. 


Comment #3: The situation in which the operating line and equilibrium curve are parallel (i.e. 
A — 1) is a special case: becomes indeterminant. However, since the driving force y*-y is 
constant, the number of steps is just the total change in concentration divided by the driving 
force: 


for A=1: N2 a7 vb 
Va" -yp* 


Now we will return to discuss the motivation for Kremser's equation: approaching the 
intersection of OL and EC: 


Notice that, as x, approaches the point of intersection (see figure on p83], yp — yp and 


- 286 | oa 


log] | i log] | 


Thus you can never reach the point of intersection for any finite number of steps. In particular, a 
very large number of stages is required as we increase the purity of either of our two product 
streams. For example, as xp—1, the end-point (xp, xp) of our operating line approaches (1,1), 
which is always one point on the equilibrium curve: as xp1, N—»oo. 


CHOOSING THE REFLUX RATIO 


An important part of the design of any distillation column is the choice of reflux ratio. There 
are a few qualitative considerations which must be kept in mind. 
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I) Existence of R min 


The first thing you should know is that there is a minimum 
value that the reflux can assume and still achieve the desired 
separation. Recall that the slope of the ROL is 


slope of ROL = M 
R«1 


Notice that, as R 1s reduced, the intersection of the ROL and the 
Q-line moves toward the equilibrium curve. As it does so, the 


minimum number of stages required approaches infinity for the 
reasons mentioned above: you can never get to the intersection of 


int of pape ie 
p pol x the operating line and the equilibrium curve. 


Usually this minimum reflux ratio can be determined using the 
point of intersection between the q-line and the equilibrium curve. 
On some peculiar xy diagrams the R.O.L. or S.O.L. constructed in 
this way might cross the equilibrium curve at points other than 
where the Q-line crosses. Then the reflux ratio must be increased 
until neither operating line crosses equilibrium curve before 
reaching the q-line. Of course, the operating line can be tangent 
to the equilibrium curve. 


Xp Zr Xp 


This point of tangency is called a pinch point. In stepping off the 

stages, it will take an infinite number just to reach the pinch, but if 
the reflux ratio is increased every so slightly, then the separation can be achieved, although the 
number of stages might be large. 


minimum reflux ratio -- smallest value of R for which the desired separation can be achieved 
with a finite number of stages. 


The pinch point might also occur on the S.O.L. as shown at right. 
Regardless of where the pinch occurs, once you have drawn the q- 
line and operating lines, the minimum reflux ratio (thus determined 
graphically) can be calculated from (11.4-24) in Geankoplis: 


x € , 
kiya (16) 
y —X 


where (x^y^) is any point on the R.O.L. (x^v? can be chosen for 
convenience: for example, you might want to choose the y-intercept. 
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To determine the minimum reflux ratio using the McCabe- 
Thiele method, we might imagine drawing the ROL and SOL 
on the xy-diagram starting with a high reflux ratio (where ROL 
and SOL coincide with 45? line) and then decrease the reflux 
ratio until any one of three things happens: 


1. the ROL becomes tangent to the equilibrium curve 
2. the SOL becomes tangent to the equilibrium curve 


3. the intersection of the ROL and SOL occurs on the 
equilibrium curve 


An important thing to notice about the minimum reflux ratio is 
that, as we approach it, the number of stages approaches infinity: 


as R > Rmin N > © 


Clearly you don't want to operate near the minimum reflux ratio because it would cost too much 
to build a column with so many stages. Generally as you increase the reflux ratio, both operating 
lines move away from the equilibrium curve so that the steps become larger. Then fewer stages 
will be needed. 


as RÎ, NL 


2) Total Reflux: R—o gives Minimum Stages 


In the limit that A approaches infinity, the y-intercept of the R.O.L. becomes zero. The 
R.O.L. becomes the 45° line. Of course, the S.O.L. and R.O.L. always intersect along the q-line, 
which means that the S.O.L. also becomes the 45? line. This 1s the furthest the operating lines 
can ever be from the equilibrium curve, so this condition gives the minimum number of stages: 


as R> o: R.O.L.  S.O.L. > 45°-line 
and N >N min 


Although it might appear that this is the best value of the reflux ratio, it turns out that there is a 
"dark side" to this solution. The problem is what is happening to flowrates inside the column. 


Recall R =L/D 
Now to get infinite reflux ratio, one of two things has to happen: 


R> o: Loo or D0 
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In most designs, you are interested in processing a fixed amount of material, this implies that 
F,D, and B are fixed, which eliminates the second alternative above. A simple mass balance 
about the split of L and D tells us that the vapor flowrate must also get very large: 


L — o: V > oo, qe > ©, q, > ©, tower diameter — oo 


This means that the heat exchanger for both the condenser and the reboiler are getting very large 
as are the cost of utilities (steam and cooling 


water). 
total 


Thus the costs of building the tower becomes ¢ 
large either near the minumum reflux ratio 
(where the height of the tower becomes large) 
or near total reflux (where the diameter of the 
column becomes large). The costs of utilities 
increase about linearly with reflux ratio. 


heating 


year & cooling 


fixed 
costs 


The total costs will experience a minimum at 
some reflux ratio, which we call the 


optimum reflux ratio: minimizes total 
costs/year 


Typically the optimum occurs for: 
1.1 < Ropt'Rmin < 1.5 


However most plants actually operate at slightly larger R since the additional cost is not very 
much and it buys you flexibility in case you decide later you want to change xp or in case the 
feed composition changes, you have some room to maneuver: 


ED 250 
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CALCULATING FLOWRATES AND CONDENSER/REBOILER HEAT DUTIES 


To compute the heat duties of the condenser and reboiler, we will need to know the flowrates 
in various parts of the column. We will start at the top of the 
column and work our way down. First we use the definition 
of reflux ratio: 


L-RD and V=(R+1)D (77) 


Given F or D and the composition of the three streams (xp, 
Zr and xg), we can determine the two unknown flowrates 
just by doing a total mole balance and a component mole 
balance over the entire distillation column (the blue system 
in and figure at right). 


F-D-«B 
ZFF — xpD + xgB 


Solving these two equations for the two unknown flowrates, 
we obtain 


Pu 9H Xg8) 
Zp —Xpg 
and go 5D VE pa mno eg (79) 
Xp —-Xp Zzp —XBg 


Given F, xp, zr and xg, we can thus calculate D, V and L. In a total condenser, we are 
condensing all of the overhead vapor. Assuming the condenser output is a saturated liquid, the 
condenser heat duty is just 


qc — VÀ 


To calculate the reboiler heat duty, we need to evaluate 
V . A look at the streams entering the feed plate (see 
figure at right), we realize that we have at least two H, = (1-4). 
unknowns L and V . A total mole balance gives: 


V+F+L=V+L 


which represents only one relationship in two 
unknowns. The remaining relationship is obtained from an enthalpy balance about the feed tray: 


FH pf + LH; +VHy = LH, + VHy 
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Our analysis of countercurrent cascades is based on the equimolar overflow assumption, which 
contains some assumptions about enthalpy: see "Constant Molar Overflow" on pri] In 
particular, we neglect sensible heat and assume that all components have the same latent heat of 
vaporization, then we can evaluate the enthalpies of all saturated liquid or saturated vapor 
streams as: 


Hy - 0, Hy-À 


Recall the definition of q from pR1] This allows us to express the enthalpy of the feed in terms of 
q: 


|Hy-Hp k-Hy 


= or Hyg-(l1-g)A 
1= Hy H, X F-— (0-4) 


Substituting the values for enthalpy into the enthalpy balance about the feed plate gives: 


E y-V 
F(1l-q. -V X- Vx. or e cbed (80) 


Substituting this into the total mole balance and dividing by F: 


L-L 
or | -= =1 (81) 


Thus introducing the feed onto a plate causes a discontinuity in the flowrates of liquid and vapor 
in the column. Although the flowrates are the same for each stage in the cascade above the feed, 
and the same for each stage in the cascade below the feed, the flowrates are different in the two 
cascades. For example, for a saturated liquid feed, q=1, and the liquid leaving the feed tray has a 
higher flowrate than than entering by an amount equal to the feed flow rate: 


for q=1: L-L=qF=F and V-V -(1-g)F-0 


whereas the vapor flow is unchanged. On the hand, if the feed is saturated vapor, then the vapor 
flowrate leaving the stage will be higher than than entering by an amount equal to the feed rate. 


for q=0: L-L=0 and V-V =F 


The flows for other thermal conditions of the feed are summarized graphically below: 
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de 


q<0 q=0 0«q«1 q-1 gl 


Knowing all the flowrates, we can now calculated the reboiler heat duty. The reboiler has to boil 
all the vapor send into the bottom cascade. For a partial reboiler (where the bottoms product is 
withdrawn from the liquid remaining in the kettle), there can be no superheating of the vapor (it 
must be a saturated vapor). Then the reboiler heat duty is 


q, =V 


ROLE OF Q-LINE 


Recall (see “Btep 4] on phi that we said one point of the 
SOL is the intersection of the q-line (QL) and the ROL. In 
this section, we present the arguments supporting this claim. 


The equations of the operating lines are obtained by 
performing a mole balance on component A between some 
pair of streams in the middle of the column and the end of the 
column away from the feed (balances use the systems denoted 
by the yellow lines): 


R.O.L.: Vy = Lx + Dxp (82) 


S.O.L.: Vy-Lx-Bxg (83) 


We claim that the q-line represents the intersection of these 
two lines. Think of these as two equations in two unknowns 
(x,y). The point of intersection of these two lines is that point 
(x,y) which simultaneously satisfies both equations. 


Now we are not going to try to solve these two equations for (x,y). But we can find other linear 
relationships between x and y by combining these two equations. One linear combination is 
particularly convenient in that it yields an equation which we can plot without knowing any of 


the flowrates that appear in (82) and 
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(83)lminus (V -Vy =(L -Dx - (Bxg-Dxp) 


Notice that the last term in this equation represents the total rate of 
component A leaving the tower. This must equal the total rate at 


which A enters the tower. So the above becomes: 7 
(F Vy - (E -Dx - Fzp 
Dividing by F gives: 
V-V Der eee 
ae: m 
—— -—— 
q-1 q 
Substituting (80)|and (q-l)y = qx -Zf 
or y= q x+ ZF 
g-1 1-q 


which is the equation of our q-line. This equation represents a linear combination of the two 
equations which are satisfied at the intersection of the two operating lines. Thus this equation 
must also be satisfied at that point. It turns out that that QL represents the locus of intersections 
of the ROL and SOL for different reflux ratios. 


PLATE EFFICIENCY 


So far we have always assumed that the vapor and liquid 
streams leaving a given stage are at equilibrium. In 
stepping off the stages, we have always gone all the way 
from the operating line to the equilibrium curve. If the 
contacting of the two phases on the tray is good (that is, if 
the area of contact is large and the time of contact is long) 
we can approach this ideal state, but we will never quite 
reach it. For a nonequilibrium stage, our horizontal step 
will be a little shorter than for an ideal stage. Of course, we 
still end up at a point on the operating line, because this 
represents a mole balance which must still be satisfied. 


Because the real steps are shorter, it will take more of them to achieve the same separation. The 
ratio of the number of ideal to real stages required will be less than one and is called the: 


overall (column) efficiency -- Nigeqi/Nreal = Eoc < 1 
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While this efficiency is easy to use, the value of the efficiency is hard to 
determine. A more fundamental quantity is the 


Murphree (tray) efficiency -- 


HM = Ey (84) 
Xn T Yn+1 
where yg = actual y leaving stage n 


y, =y in equil with x, 


On an xy diagram, this represents the height of a real step 
divided by the height of an ideal step for the same x,. 
Alternatively, one could also define efficiency based on the 
change in liquid composition (rather than vapor composition): 


Ny aX 
Ey, =—=1—2 (85) 
X 


n-l "n 


This definition is based on the width of a step (involves 
changes in x) rather than based on the height of a step 


(changes in y) as in 


Some typical values of plate efficiencies are given in Table 18-4 (taken from Perry's 6" edition). 


06-361 page 97 Spring, 2001 


Table 18-4 Representative Plate Efficiencies 


Column Tray Static 
diameter, spacing, submergence, 
Disperser System ft in in Remarks Ref. 
Bubble cap Ethanol-water 1 
2 
Methanol-water 3.2 15.7 147 3 
Ethyl benzene-styrene 2.6 19.7 19 4 
Cyclohexane-n-heptane 4.0 24 147 5 
24 
50 
Cyclohexane-n-heptane 4.0 24 5 Tunnel caps 6 
24 
Benzene-toluene 15 15.7 14.7 n nhe 7 
Toluene-isooctane 5.0 24 14.7 || | ...... 8 
Ripple sieve Methanol-water 3.2 15.7 14.7 10.8% open 3 
Ethanol-water 25 14 14.7 10.4% 2 
Methanol-water 3.2 15.7 14.7 4.8% open 3 
Ethy! benzene-styrene 2.6 19.7 19 12.3% open 9 
Benzene-toluene 15 15.7 14.7 18% open 1 
Methyl alcohol-n-propy! alcohol- 
sec-butyl alcohol 6.0 18 18 10 
Mixed xylenes + Cg-Cig paraffins 
and naphthenes 13.0 2] 25 5 
Cyclohexane-n-heptane 4.0 24 5 14% open 13 
24 14% open 13 
4.0 24 5 8% open 12 
Isobutane-n-butane 4.0 24 165 14% open 13 
4.0 24 165 8% open 12 
40 24 300 896 open 12 
4.0 24 400 8% open 12 
Turbogrid valve Methanol-water 3.2 15.7 14.7 : 3 
z Ethanol-water 25 14 14.7 | ES eee 2 
Ethyl benzene-styrene 2.6 19.7 19 : ier 4 
Cyclohéane. n-heptane 40 2 20 50-96 Rect. valves " 
n-Butane-isobutene 40 24 165 104-121 Rect. valves 
Benzene-toluene 15 15.7 14.7 Satna y 2d 
References 


I. Kirschbaum, Z. Ver. Dtsch. Ing. Beih. Verfahrenstech., (5), 131 (1938); (3), 69 (1940). 

2. Kirschbaum, Distillier-Rektifiziertechnik, 4th ed., Springer- Verlag, Berlin and Heidelberg, 1969. 

3. Kastanek and Standart, Sep. Sci, 2, 439 (1967). 

4. Billet and Raichle, Chem. Ing. Tech., 38, 825 (1966); 40, 377 (1968). 

5. AIChE Research Committee, Tray Efficiency in Distillation Columns, final report, University of' Delaware, Newark, 1958. 

6. Raichle and Billet, Chem. Ing. Tech., 35, 831 (1963). 

7. Zuiderweg, Verburg, and Gilissen, Proc. Intn. Symp., Brighton, England, 1960. 

8. Manning, Marple, and Hinds, Ind. Eng. Chem., 49, 2051 (1957). 

9. Billet, Proc. Intn. Symp., Brighton, England, 1970. 
10. Mayfield, Church, Green, Lee, and Rasmussen, Ind. Eng. Chem., 44, 2238 (1952) 
11. Fractionation Research, Inc., "Report of Tests of Nutter Type B Float Valve Tray,” July 2. 1964, from Nutter Engineering Co., Tulsa. 
12. Sakata and Yanagi, Inst. Chem. Eng. Symp. Ser., no. 56, 3.2/21 (1979). 
13. Yanagi and Sakata, Ind. Eng. Chem. Process Des. Dev. 21, 712 (1982). 
“See Eq. (18-32). 
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NOTE: To convert feet to meters, multiply by 0.3048; to convert inches to centimeters, multiply by 2.54; and to convert pounds-force per square inch to kilopascals, multiply by 6.895. 


06-361 page 98 Spring, 2001 
DETERMINING NUMBER OF REAL STAGES 


Suppose we have a value for Emy or Emy, how do we use it to determine the number of real 
stages required? The procedure is similar to that for ideal stages: 


Steps 1-6) same as for ideal stages 


We locate the two operating lines and the q- 
line, just as before. But before we step off 
the stages we use the efficiency to determine 
a new curve between the operating lines and 
the equilibrium curve. 


Step 7) Plot psuedo-equilibrium curve using 


either | = or E i 


The blue curve in the xy-diagram at right 
was drawn using . At several x-values 
between xg and e draw a vertical line 
between the lowest operating line and the 
equilibrium curve. On each of these line 
segments, we locate a point which is a 
fraction of the distance up the line segment 
from the bottom, where that fraction equals the plate efficients. The locus of these points is the 
psuedo-equilibrium curve. 


Finally, we step off the stages as before, but using the psuedo-equilibrium curve in place of the 
true equilibrium curve. 


Step 8) Step off stages between operating lines and pseudo-equilibrium curve. 


Since we used definition (84)|based on the vapor mole fractions, we should start our stepping 
process at the bottom of the column, rather than the top. In other words, we start at (xp, xg) and 
step upwards until either the vapor or the liquid mole fraction exceeds xp. By stepping from the 
bottom up, the step from the operating line toward the equilibrium curve involves changes in y, 
whereas when we step from the top down, the step from the operating line toward the 
equilibrium curve involves changes in x instead of y. The latter is appropriate when the 
efficiency is defined in terms of x instead of y. 
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Chapter 12. Multi-Component 
Distillation 


The objective of distillation is the separation of the feed into 
streams of nearly pure products. In binary distillation, we have two 
components and two product streams. With more than two 
components in the feed, there is no way to achieve this with a single 
column (a single column might have a third product stream, but it's 
composition would never be nearly pure). 


A,B 


To separate a mixture of three components into A 
three nearly pure streams, we will need at least two 
separate distillation columns to get three streams. In A.B.C — B 
general, if there are N, components, we will need V.-1 — ' ' PM ae 


columns. Lx 
Nc components C 


N,-1 columns Alternative #1 


This assumes that each column can achieve a sharp 
separation between two components of adjacent 
volatility. 


Alternatively, the designer might try to first split off A 
component C in the first column — rather then A. The A,B 
choice between these two alternatives is a matter of > 
experience. One “rule of thumb” is to perform the B 
easiest separations first. In this course, we will focus A,B,C 
instead on the design of the individual columns after | C 
this decision has been made. 
Alternative #2 

The first step in the design of a single column is to Aine welatiie = A 

specify the product distribution. Let A,B,C,... denote (lowest b.p.) M 


components ranked in order of decreasing volatility. distillate 


to ; A t 
sharp split -- distillate contains only components A-M, while produ 


bottoms contains only components N-Z " 

In any real distillation, all components will be found in all N 

product streams in at least trace amounts. When we say “only” l . | bottoms 
components A-M are found in the distillate, we usually mean Gee B) ` [ product 
that the mole fractions of the remaining components (1.e. N-Z) ai LE 

can be neglected when summing the mole fractions to unity. 


Sharp separations are not always possible, but we will focus on 
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them in this introduction, for the sake of simplicity. In the description above, components M and 
N play a special role and are called 


key components -- M,N (for sharp split) 
The more volatile of the two is called the light key, while the other is the heavy key: 
light key -- M 


heavy key - N 


PRODUCT DISTRIBUTION FOR A SHARP SPLIT 


EXAMPLE 71: A mixture with 33% hexane, 37% heptane and 30% octane is to be distilled to 
give a distillate product of 1% heptane and a bottoms product of 1% hexane. Find D/F, B/F, and 
Xpj> Xp; for the other components. 


Solution: First we arrange the list of components in order of relative volatility. In this case, that 
is easy since we have a homologous series* of alkanes. Since vapor pressure general drops with 
molecular weight in any homologous series, we don't need to consult tables of boiling points or 
vapor pressures to order the components. 


Component Formula 
hexane C&Hj4 most volatile 
heptane C7H 16 
octane CgHig least volatile 


Let's take a basis of 100 moles of feed (F=100) and try to fill out the following table: 


Component X5: 

hexane 0.01 light key 
heptane -- heavy key 
octane ES 

totals — 100 1 D 1 B 1 


Comment: Generally the more volatile components end up in the distillate stream and the less 
volatile components in the bottoms. Being the most volatile, most of the hexane is expected to 
end up in the distillate; we are told that the bottoms actually contains a little hexane (xgg = 0.01). 
Heptane is of intermediate volatility and so might end up in either or both streams. From the 


"A “homologous series" is a set of very similar items, which differ in some property. In this case, the series is a set 
of components of similar molecular structure (the alkanes) which differ only in the number of CH» groups in the 
chain. 
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problem statement, we also know that the distillate contains a little heptane (xg7 = 0.01) and so 
presumably most of the heptane ends up in the bottoms. Being less volatile than heptane, we can 


reasonably assume that all of the octane ends up in the bottoms. 


In essence, we are then making a sharp split between hexane 
and heptane with the “1%” values specifying the sharpness of 
the split desired. In other words, hexane has been choosen as 
the light key and heptane as the heavy key. If we assume that 
there is virtually no octane in the distillate (octane is heavier 
than the heavy key) then we can finish specifying the 
composition of the distillate stream. 


Component 
hexane 
heptane 
octane 


totals = 100 1 D 1 


C, 


Co Ca C, 


C,, C, 


Knowing the concentration of hexane in both of the product streams allows us to determine the 


flowrates from mole balances: 


component C¢: 0.99D + 0.018 = 33 
total: D+B=100 
Solving simultaneously: D = 32.65, B = 67.35 
Component 
hexane 32.32 0.99 
heptane 0.01 
octane 
totals = 


Now using the fact that the balance of the distillate and bottoms streams must be heptane, we can 


fill out the table: 


06-361 page 102 Spring, 2001 


Component Bxp; Xp 
hexane 32.32 0.99 | 0.67 0.01 
heptane 0.33 0.001 | 36.67 0.5445 
octane 30 0.4454 
totals — 100 1 32.65 1 67.35 1 


Comment #1: In this example we have completed the specifications of the two product streams 
for a 3-component mixture in which the sharpness of split was specified by giving the mole 
fraction of the heavy key in the distillate stream and the mole fraction of the light key in the 
bottoms. Alternatively, we might have been given the recovery** of the two key components in 
either product stream. 


Comment #2: The only additional information that is needed if the number of components is 
larger than 3 is the mole fraction of those components in the feed. The initial guess for the 
distribution of additional components in the product streams is determined by the following 
rules: 


Rule 1) Components heavier than the heavy key 
are assumed to end up in entirely in the bottoms 


Rule 2) Components lighter than the light key are assumed to end up in entirely in the 
distillate 


MINIMUM NUMBER OF TRAYS 


Now let's see if we can get some idea of how many trays we will need. This is a difficult 
problem to do precisely. What can be done relatively easily is to estimate the number of trays 
needed at total reflux -- this is the minimum number of trays. Some gross assumptions: 


total reflux: D=0, L«V 
equal molal overflow: L=V=const for all trays 
constant relative volatilities: a,=const for all trays 


K. 
Recall cce. (86) 
Kj 


is the relative volatility of component i compared to j, where 


* For example, “5% of the heavy key is recovered in the distillate” means Dx 'D HK ^ 0.05xF ZF gy. 
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K, 224 (87) 


l 
i 
is the distribution coefficient for component i (or j). With 


these assumptions, we can start at the top of the column and 


work out the compositions of streams leaving every stage. D-0 
For a total condenser: *pi 


Vii = Xoi = Xpi (88) 


for i=1,...,.N.. The first subscript is the stage number (as in 
binary mixtures) and the second subscript is the component. 
With all the compositions of the vapor stream leaving stage 1 
known, we can determine the compositions of the exit liquid 
stream from the known values of the relative volatility: 


n 
^ Kj e 
Xlj 


" Xi loas y, i (89) 
l 
Xij Qj yj Yij 
Substituting E TH. O ji “Di (90) 
Xij x Dj 


"e Xli X Di 
for i=j: “Vig, PL gp 1=1 
B qu 
1 1 


which doesn’t tell us anything new; thus we have only N,-1 independent relations of this type 
instead of N.. Choosing one element as the reference (say j-N,.), we can write N-1 relations for 


the N,.-1 independent mole fractions. The missing equation is Xx; = 1. Thus the composition of 
the liquid leaving tray 1 has been completely determined. 


Next we determine y»; by mole balances. 
component i: Vy»; = Dxpj + Lxjj 


but for total reflux D~0 so that LV, leaving: 
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J2i ^?4i (91) 


for i-1,..,N.. Knowing the composition of the vapor leaving stage 2, we can calculate the 
composition of the liquid from the known volatilities. 


Rewriting|(89)| replacing the stage number 1 by 2: 


*2i Y2i 

X2j Y2j 

x9; Xij 

from konl PEE y B 

X2j Xlj 

: *2i 2 «pi 
Using sey! = Qj (92) 

X2j X pj 


From (90)]and we begin to see a trend. The generalization of this trend yields: 


Xi Xni 
ni _ a^ Di 
Xnj x Dj 


In particular, we are interested in the number of steps required to get from the distillate to the 
bottoms: 


X Bi = N min * Di 
= Ob jj —— 
x Bj Xx Dj 


The exponent represents the minimum number of steps on the xy diagram. One of these steps is 
the partial reboiler. Solving for N,,;,: 


In XDi/XBi 
XDj/¥Bj 


which is called Fenske's equation. Now in any real system, the relative volatility will not be 
constant. One way to patch things up is to use some kind of average volatility: 


Qij = 49 pij. Bij 


where aj; and aj; are evaluated at the bubble-points of the distillate or bottoms, respectively. 
In this case, an appropriate mean is the geometric mean of the volatilities evaluated at the 
conditions of the distillate and bottoms compositions. We use the bubble-point temperature 


Nmin = 
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rather than the dew-point temperature because the distillate and bottoms streams are usually 
saturated liquids. If they were saturated vapors, then we would use the dew-point instead. 


Example #2. Calculate the Nin for the separation in Example #1 if the column is operated at 
1.2 atm. 


Solution: 


Evaluating the relative volatilities 


We start by evaluating the relative volatilities in the feed, distillate and bottoms streams knowing 
their composition and thermal state. 


Time out: Vapor pressure data for pure components is available in the CRC handbook. C6, C7 
and C8 form a nearly ideal mixture and Raoult's law applies: 


pi = XP? 


where p; is the partial pressure and P;° is the vapor pressure. Then the distribution coefficient can 
be calculated from 


EJ. p. Xp 4E (93) 


A y 


and the relative volatility from Qj; = = = 


Vapor pressure depends very strongly on temperature: recall the Clausius-Clapeyron equation 
(S&VN, p182): 


daR? o hi din PP X 
dT — RT? d(/T) R 


where i; is the molar heat of vaporization of the pure component. Integrating this with respect to 
T (ignoring the weak temperature dependence of à;) yields 


Ek (95) 


js 
or pz! ex[- AL.) (96) 
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where A; is just some integration constant. Note the exponential sensitivity of vapor pressure to 
temperature. Vapor pressure data is often summarized in handbooks (e.g. Lange's Handbook of 
Chemistry, 13" edition, pages 10-29 to 10-54) using Antoine's equation: 


Bi (97) 
T+C; 


which is a semi-empirical generalization of (95)] that accounts for the variation of à; with T. 
Other handbooks (e.g. CRC Handbook of Chemistry & Physics, 57" edition, pages D-183 to D- 
215) give vapor pressures at several temperatures and leave you to develop your own 
interpolation scheme. 


Since vapor pressure depends strongly on temperature, so does the partition coefficient K; given 
by However, the relative volatility a is much less dependent on temperature. The reason 
becomes clear when we substitute (96) into 


NUES 
Qj = = , EXP 


J 


Note that when the two molar heats of vaporization are equal, the temperature-dependence 
vanishes. Even if they are only approximately equal, the temperature-dependence of the relative 
volatility will be much less than of the vapor pressures. For this reason, 


the assumption of constant relative volatility is frequently used. 0,990 
Time in: We took the values of the vapor pressure of each component onn 
given in CRC handbook and fit them to (97)|to find the constants 4;, B; 0.33 
and C; for each component. 0.37 
0.30 
Now the distillate is a saturated liquid containing nearly pure hexane, 0.010 
so its temperature is near the boiling point of hexane at 1.2 atm (the 0.545 
pressure throughout the column) or 75?C. 0.445 
C t P; (at o K; 
omponen Xp j? (atm) K = A ie Ve 
C6 0.990 1.200 1.000 2.574 
C7 0.010 0.466 0.388 1.000 
C8 0.000 0.187 0.156 0.401 


The bubble point of the bottoms is 115°C (see Geankoplis p681-2 for how to determine bubble 
points for multicomponent mixtures): 
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Component Xp Pj? (atm) p? m e 
K; = P Oj HK ~ Kuk 


Calculate Minimum Stages 


To calculate Nin, Fenske's equation is usually applied to the heavy key and the light keys. So 
we are interested in a y HK 


067 = 42.574 x 2236 = 2.399 


Fenske's equation is 


1 d | 099/0.01 ) 
X pj [Xj 0.01/0.544 
Naw YN, / - 982 


In(o.;;) ... 1n(2399) 


So at an absolute minimum, we need 10 ideal steps, or 9 ideal stages plus the reboiler. Knowing 
the number of trays, we can also use Fenske's equation to check if indeed the octane fraction in 
the distillate is negligible as we assumed in Example 71. 


&g7 = V0.401 x 0469 = 0.434 
| Ms ie] 


w -1o (001/0544 
mm ^^ 1n(0.434) 
Solving for xpg, we obtain Xpg = 1.93x10-6 


which indeed is negligible compared to the other mole fractions in the distillate. 
MINIMUM REFLUX RATIO 
Method 1) Psuedo-Binary Method 


Assuming a saturated liquid feed (q=1) and constant relative volatility, the minimum reflux 
ratio turns out to be 
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1- 
=F zd (98) 
where xp and zp are the molefraction of the light key in the psuedo-binary mixture and 1-xp and 
l-zp are the molefraction of the heavy key in the psuedo-binary mixture: 


ZF LK ZF HK 


ZF = and 1 ZF = 
ZF LK tZF,HK ZF LK tZF,HK 


where zp;y and Zp yg are the mole fractions in the multi- 
component mixture. 


Proof: We can obtain the minimum reflux ratio from the point of |, "4 
intersection of the q-line and the equilibrium curve. Recdll (16) 
from page 89! 
xp—y' 
Race (99) 
y'-x 


which was derived using two points on the ROL to determine the 
slope and knowing that the slope is R/(R+1). For a saturated s Eg Xp 
liquid feed, the q-line is vertical so that 


x'2zp (100) 


The y-coordinate of this point is determined knowing that ( x', y' ) lies on the equilibrium curve. 


For constant relative volatility, we can express the equilibrium curve as a simple equation using 
the definition of relative volatility: 


Ky _yalxa _ Y 


U AR = 
Kg ypgixg x[(l-x) 
Solving for y: y= = 
14 (a. — Dx 
Substituting x = x' = z p , we can calculate 


1-4 (a 4p - zp 


* To derive the relations above, take a basis of 1 mole of feed. then zp z g and zp yg represent the number of moles 
of LK and HK in the feed. 
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Substituting (101)|and (100) [into [99)]leads to 
Method 2) Underwood's Method 


This method makes similar assumptions to the psuedo-binary method except that all 
components are considered, not just the two key components. In particular, the assumptions 
made by Underwood's method are: 


1) equimolal overflow 
2) constant relative volatility 
3) asharp split with only one “pinch” 


First, we will just state the method and then we will take a quick look at its proof. 


Recipe: 


Step 1) Find the root (denoted $4) in the interval aye «6 «ark which satisfies the 
polynomial (order of polynomial = number of components) 


y SEP ela (102) 
i Qj -9 


Step 2) Once we have $4, we can compute the minimum reflux ratio from 


Rmin +1= yee (103) 


where the sum is just over those components present in the distillate (1.e. we ignore components 
heavier than the heavy key). 
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Proof: A mole balance on component i about the top n stages of 
the rectification section yields: 


Veni = Lyn i + Dxp,j (104) 
Recall that the distribution coefficient is defined as: 


Yn,i 


^ 


K;-yix; or K,;— 
Xn i 


We put the tray subscript on K because its value depends on the 
temperature, which varies from tray to tray. Thug (104) can be 
rewritten as: 


Yn,i 
Vna+1Yn+Li =L, *Dxpj 


n,i 


For constant molal overflow, the flowrates will be the same for 
each tray in the rectifying cascade, just as with binary mixtures. 
Thus we can drop the subscripts on flowrate. 


Yy 5 
Vy = LE * DxD (105) 


n,i 


When we approach the minimum reflux ratio, we will encounter a "pinch" somewhere in the 
column. At the pinch, we take an infinite number of steps without changing the mole fractions: 


at a pinch: Yt i 7 Yn (106) 


(106)|into and solving for y;: 


D * i 
— 107 
Vin Kooi 


We have put a subscript “oo” on each flow rate to emphasize that this corresponds to the 
minimum reflux, which requires an infinite number of trays. If we require to be satisfied 
for each component and sum over all components: 


D X pj 
3-197 FG (108) 
i dum — 


Now we can write K; in terms of the relative volatility: 
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Q; = 
Kref 


While the flowrate ratios can be related to the reflux ratio: 


Lo Rmin g 22! 
Vo Rami Vo Rmint! 
aes becomes: paid = Rmin +1 aes) 
a; — 


where 


If Kep were known, then Rmin could be determined as the root of this equation. Unfortunately, to 
calculate Kep we need to know the temperature at the pinch, which is not known. 


Generally there is also a pinch in the stripping section of the tower. For sharp splits, the 
temperature at this second pinch will be the same. We can then derive a second equation like 
103)|by performing our balances about the bottom n stages of the column to obtain a second 
equation. 


where the overbars denotes flowrates in the stripping section. Then we will have two equations 
and two unknowns: K,ef and Rmin- Their solution is the value of $ which satisfies the following 
equation: 


where zp; is the composition of the feed and q is its thermal condition. There are several 9's 
which satisfy this equation. It turns out that we want: 


QHK <>< ALK 
Once 6 is determined, we can calculate R,,,;, from This is called Underwood's method. 


Example #3. Calculate the Rmin for the separation in Examples #1 and #2 if the feed is 60% 
vapor (q = 0.4). 
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Solution: Recall the feed and product compositions (see sketch at 


right). Using Underwood's method, we first determine the root of 0.990 
102) 0.010 


0.000 
"P 0.33 
> = 1-4 0.37 
7 trý 0.30 
0.010 
Substituting into known values: One 


2399x033 1x037 0.434x 030 
2399-6 1-6  0434-¢ 
"nc 


=1-04 


which is a cubic equation. The left-hand side 
Ko) is plotted at right. This equation has 
three singular points at 


= a4, o and ag mE I- 


These are shown as the 3 vertical lines in the 
plot. Notice that there are also 3 different $'s 
for which mS 


fe) = 1-q - 0.6 | | | 


We are interested in the root in the interval: 


OG gy «6 «ark or 1<o<2.399 


which turns to be $ = 1.739. Once we have the value of $, we calculate the reflux ratio from 
103) 


where the sum only includes those components present in significant quantities in the distillate. 
In this case, octane (which is heavier than the heavy key) was neglected in the distillate for 
Example #1. It should also be neglected in this sum, although in this case, it’s inclusion does not 
change the result. 


_ 2399 x 0.99 " 1 x 0.09 


ntl= = 359 
mun 299-1739 1-1.739 


ANSWER: Rar 29 
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NUMBER OF IDEAL PLATES AT OPERATING REFLUX 


Ref: MSH p588-609 


Method 1) Underwood's method 


Underwood's method can be extended to calculate the number of plates analytically when the 
equal molal overflow assumption and constant relativity are appropriate. But even then, the 
calculations are sufficiently involved that a 
computer would be helpful. If we are going to 
use a computer, we might as well do the rigorous 
tray-by-tray calculations (see Method 3); thus we Lac 
will not extend Underwood's method further in 
these Notes. 


Method 2) Gilliland correlation 


An empirical but much easier to use method is 
the Gilliland correlation. This is an empirical 
relation between the number of ideal stages and 
the operating reflux ratio. You also need to 
know the minimum number of ideal stages (at 
total reflux) and the minimum reflux ratio. 
Because of its simplicity, this correlation is 
widely used for preliminary estimates. 


DH 002 00d C10 Oe 
Fale 
1 
Re 


A second type of empirical correlation is Erbar dL clc 


& Maddox (see Fig. 11.7-3 in Geankoplis): 


Nmin R Rmin 
Rel Rmin +1 


Method 3) Rigorous Solution of MESH Equations 


Reference: Henley & Seader, Equilibrium-Stage Separation Operations in Chemical 
Engineering, Chpt. 15. 


The final design of multistage equipment to multicomponent distillation usually requires 
rigorous determination of temperatures, flow rates, and compositions at each stage. This must be 
done numerically by computer. We will use a computer program called CHEMSEP to look at a 
couple of typical columns. Although the program is designed to be used like a "black-box", I'd 
like to give you a little peek at what goes on inside. 
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A rigorous mathematical description of a multistage cascade involves material and energy 
balances and some representation of the vapor-liquid equilibrium. The resulting set of algebraic 
equations is called the MESH Equations, for reasons which will soon become apparent. 


Consider a general stage in the cascade, which we 
will label stage n. It has the usual liquid and vapor 
streams entering and leaving. It might also be a feed 
stage, in which case F,, z 0. If it is not a feed stage, 
we just set F, — 0. There might also be liquid or 
vapor product streams U, or JW, as well as a heat 
removal Q,. The MESH equations result from 
balances done about this arbitrary stage. 


M Equations are Mole balances for each component i 
(N, equations for each stage): 


Ln-1Xi,n-1 + Vn+1Yi,n+1 + FnZi.n = (Ly *U, om + (V, +W,, Win 
Note that the first subscript is the component i and the second subscript for the stage n. 


E Equations describe the liquid-vapor Equilibrium relations among the two outlet streams, which 
are assumed to be at equilibrium (N, equations for each stage): 


Xin = Ki nXin 
where K; , is the distribution coefficients. 


S Equations are the mole fraction Summations (2 equations for each stage): 
ECT -] and 3 dux =1 
i i 


H Equations is the energy balance for each stage (1 equation for each stage): 


L 


nl psp Vn Pe Pin = (y PU Hn Vat ea to 

If there are a total of N, stages in the column and N, components, then the number of equations is 
number of equations = N,(2N, +3) 

Now we consider the following to be unknowns: 

Xim Yin Lm Vm Tn for i= La, and nN, 


where T, is the temperature of the liquid and vapor streams leaving stage n. 


number of unknowns = N (2N; +3) 
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Of course, there are more variables involved in the equations. We assume that the feeds and any 
side drawoff flows are specified: 


given: ym Zin HF n Un, Wr Qn 


Also, we assume that the thermodynamic relations exists so that we can compute the enthalpy of 
saturated liquid and saturated vapor, given the composition, pressure and temperature: 


known functions: Hin =HL,n (T, P, xis a 


Hy n = Hy „(Tp Pp Vin YN. n) 
Kin = KDB ya) 
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Chapter 13. Computer Simulation of Distillation 
We will be using a computer program called *ChemSep" to simulate distillation. Today I'd like 


to introduce this program. Compared with industrial software, it is user-friendly and written for 
student use. Below is the welcoming screen. 


Arno Haket and Harry Kooijman are from the University of Technology in Delft (The 
Netherlands), while Ross Taylor is a Professor of Chemical Engineering at Clarkson University 
(New York state). 


After a copyright screen, we reach the main menu: 
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Results 


nits F Options F InoF 


SOME PRELIMINARY STEPS 


Although a mouse can be used with the program, the default is to use the cursor keys to move 
around the screen. When the desired menu item is highlighted, press ENTER to select it. You 
can also select menu items by typing the first letter of the desired item. Below is the result of 
choosing the File menu. 


Results 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 
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Notice that a sub-menu has appeared. You can probably guess what most of these commands do. 


Before going any further, you should make sure that any files we create will be placed on the 
Floppy Disk drive A:\. This can be accomplished by changing the default directory to A:\. 
Select the Directory command. You will then be presented with the current directory name 
which can be edited to read: 


Options Help 
Chemsep v3.50 


Rename 
Erase 
DOS shell 


:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


You should also save your configuration file on your floppy disk. This can be accomplished by 
going to the options menu and selecting the Save options command. You should edit the 
file name to read “A: NChemSep.cn£": 


Results Options Help 
Chemsep v3.50 


Solve options 
Units 

Macros 
Interface 
Directories 
Video device 
Output device 
Save options 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info Fl0:Main menu 
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Now any changes you make in the configuration can be saved in this file and re-loaded when you 
return to ChemSep at some later time. 


LOADING A FILE (PROBLEM) 


We are going to create a new file which simulates the solution to Prob. 3 on Hwk. #6. We start 
by selecting the New command from the File menu. You will then prompted for a file name: 
type “A:\HWK6P3.SEP”. The "A: V" will force this file to be stored in the root directory of the 
floppy drive, which is the recommended location for your personal files. 


File Input Results Options Help 
A:NHWK6P3.SEP (New problem) Chemsep v3.50 


Components 
Operation 
Properties 
Specifications 
Solve 

ChemProp 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info Fl0:Main menu 


The first step in inputing the problem description is to select the components. This particular 
homework problem is based on MS&HS Prob. 18-3 in which we are to separate a mixture of 
acetone and methanol. Under the input menu, select Components, followed by several 
“Enter’s” until you are prompted for a name (see below): 
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File Input Options Help 
A:NHWK6P3.SEP (New problem) Chemsep v3.50 
--« Components »- 


| None 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


Enter the first component name (acetone) and hit “enter”. This will bring up another short menu 
containing “acetone” at the bottom: Select “acetone” and hit “enter”. 


File Input Results Options Help 
A:NHWK6P3.SEP (New problem) Chemsep v3.50 
--« Components »- 


-» Restart search 
-» Refine search 
Acetone 
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This will add acetone to the list of components. Next select ->Restart search and type the 
next component name (methanol). After selecting methanol and hitting “enter”, we have two 
components in our list: 
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File Input 


Options Help 
A:NHWK6P3.SEP (New problem) 


Chemsep v3.50 
--« Components »- 


| Acetone 


| Methanol 


-» Restart search 
-» Refine search 
Methanol 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


Hit “Esc” twice brings us back to the Input menu. This time, select the operation command, 
then the Column subcommand. The cursor is now on the “*” 


on the line with operation. 
Hitting the “enter” key leads to the following screen: 
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File Input Options Help 
A:\HWK6P3.SEP (New problem) Chemsep v3.50 


| Components 
| Operation 


Condenser Simple Distillation 
Reboiler Extractive Distillation 
Stages Azeotropic Distillation 
Feed stages Simple Absorber/Stripper 
Sidestream stages; Reboiled Absorber/Stripper 
Pumparound stages; Refluxed Absorber/Stripper 

Single Column Stage 

Simple Extractor 

Single Extraction Stage 

Complex Column 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info Fl0:Main menu 


Move the cursor to simple distillation and hit “enter”. This moves the cursor to the ^" on 
the line with Condenser. Hitting “enter” brings up a menu of options for the operation of the 
condenser: 


File Input Results Options Help 
A:NHWK6P3.SEP (New problem) ChemSep v3.50 


Components 
Operation 


Condenser 

Reboiler 

Stages Total (Liquid product) 
Feed stages Total (Subcooled product) 
Sidestream stages; Partial (Vapour product) 
Pumparound stages; None 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


Of course, in the homework problem we want to subcool the condensate, so we choose Total 
(Subcooled product). After making a few other choices in the remaining items, then selecting 
Show flowsheet, we have: 
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File Input Help 
A:NHWK6P3.SEP (New problem) Chemsep v3.50 


Components 


Operation 


Condenser Total (Subcooled prj 
Reboiler Partial (Liquid proi 
Stages 23 

Feed stages 17 

Sidestream stages 

Pumparound stages 


Show flowsheet 


23+----+ 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


We next select Return and the cursor moves to the Properties submenu of Input. Hit “enter” 
twice to enter the Thermodynamic models dialog box. Since this is not a course in 
thermodynamics, I will always tell you which thermo models to us. After making a few choices, 
the dialog box appears as follows: 
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File Input Options Help 
A:\HWK6P3.SEP (New problem) Chemsep v3.50 


Components 
Operation 
Properties 


Activity coefficient 

Equation of State Peng-Robinson 
Vapour pressure 

Enthalpy Excess 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info Fl0:Main menu 


Selecting Return takes us to the Load data submenu. Some thermo models do not require any 
further input from the user, but the cubic EOS does require a value for the interaction parameter. 
After selecting Load data, then Cubic Eos leads us to the following dialog: 


File Input Results Options Help 
A:NHWK6P3.SEP (New problem) ChemSep v3.50 


Components 
Operation 
Properties 


+ 
Thermodynamic models 
Load data 
--+ — ee ee ee ee ee 
Load data for: 
4--|---------------- 
| Cubic EOS 


Component i competent j 


Acetone Methanol 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


Move the cursor to Correlation and hit "enter". A value for kij has been entered from the 
library, then select Return. Hitting “enter” three more times brings us back to the 


Input/Specifications submenu. Hit “enter” and choose Pressures. Make the selections 
shown as follows: 
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File Input Options Help 
A:\HWK6P3.SEP (New problem) Chemsep v3.50 


Components 
Operation 
Properties 
Specifications 


Analysis 
Pressures 


Condenser pressure 

Column pressure Constant pressure 
Top pressure 1.000 

Pressure Drop (atm/stage) 

Bottom pressure 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info Fl0:Main menu 


The units for pressure might not show up as atm. If not, you can change the units for entering 
any quantities hitting the function key r5. After setting the pressure, choose Return which moves 
the cursor to Heaters/Coolers. No heat is to be added or removed except for the condenser 
and reboiler, so we move the cursor on down to Efficiencies. Here we enter the Murphree 
tray efficiency. We enter “1” as the default to obtain ideal trays throughout the column. 


File Input Options Help 
A:NHWK6P3.SEP (New problem) ChemSep v3.50 


Components 
Operation 
Properties 
Specifications 


Analysis 
Pressures 
Heaters/Coolers 
Efficiencies 


Efficiency (-) 
Default 1.000 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


Using the Insert command, we could specify different efficiencies for each tray. Choosing 
Return moves the cursor to the Input/Specifications/Feeds command. This is where we 
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specify the thermal condition, composition and flowrate of the feed. The thermal condition can 
be specified either as temperature and pressure or as vapor fraction and pressure. Since our feed 
is 70% vapor, we use the latter option. After entering the values appropriate to our problem, we 
have: 


File Input Results Options Help 
A:NHWK6P3.SEP (New problem) Chemsep v3.50 


Components 
Operation 
Properties 
Specifications 


Feed stage 
State: 

Pressure (atm) 
Vapour fraction (-) 
Temperature (°C) 


Flowrates (mol/s) 
Acetone 
Methanol 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


A second feed could be added using the Insert command. Selecting Return, brings us to the 
Input /Specifications/Operation command. This is where we enter specifications regarding 
the two remaining degrees of freedom in the column design. Hitting "enter" twice brings up a 
range of choices for specifying the first degree of freedom. 
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File Input Options Help 
A:NHWK6P3.SEP (New problem) Chemsep v3.50 


Components 
Operation 
Properties 
Specifications 


Analysis 
Pressures 
Heaters/Coolers 
Efficiencies 


Reflux ratio 
Heat duty of condenser 


* 


Distillate flow rate 


[wo] 


Reflux flow rate 


* 


Component flow 


Mole fraction of a component 


Component recovery 


Fraction of combined feeds recovered 


+ 
+ 
l 
I 
I 
I 
I 
1 
I 
I 
I 
I 
I 
i} 
1 
1 
I 
1 
I 
i 
+ 
+ 


Split between two components 
Flexible 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


We can choose only one of these. If we specify the mole fraction of acetone to be 78% (as in our 
problem statement), then the program adjusts the reflux ratio in an effort to obtain this 
composition. If we specify the reflux ratio, then the program determines the composition of the 
distillate product. Here we will choose the reflux ratio (to be 8.62) and see what happens. Since 
we specified a subcooled reflux, we next get to choose its temperature or the number of degrees 
of subcooling. We specify the temperature as 25?C. Hitting "enter" when the cursor is over the 
only remaining “*” brings up the following list: 
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File Input Options Help 
A:NHWK6P3.SEP (New problem) Chemsep v3.50 


| Components 

| Operation 

| Properties 

| Specifications 


Analysis 
Pressures 
Heaters/Coolers 


Boilup ratio . 
Heat duty of reboiler 
Temperature of reboiler 


ES 


Bottom product flow rate 


- 


Reboiled vapour flow 


* 


Component flow 


Mole fraction of a component 


Component recovery 


Fraction of combined feeds recovered 


+ 
+ 
l 
I 
I 
L] 
I 
1 
I 
l 
I 
I 
I 
i} 
1 
1 
I 
I 
I 
i 
+ 
+ 


Split between two components 
Flexible 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


Now we get to specify one property of the bottom stream. We choose to specify the mole 
fraction of the bottoms product. Now the Operation specification is complete. 
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Input 


A:NHWK6P3.SEP (New problem) 


+ +--------------------+ + 


Components 
Operation 
Properties 
Specifications 


Analysis 
Pressures 
Heaters/Coolers 
Efficiencies 
Feeds 

Operation 


Reflux ratio (-) 
Temperature of reflux (^C) 


Mole fraction of a component (-) 


Options 


Acetone 0.010000 


Spring, 2001 


Help 
Chemsep v3.50 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info Fl0:Main menu 


Choosing Return twice places the cursor on the Input/Solve command. Hitting "enter", the 


program saves the description and attempts to solve the MESH equations. If all the data has been 
entered correctly, the next screen appears: 


File 


Input Results 


A:NHWK6P3.SEP 2 seconds 5 iterations 


Options 


Help 
Chemsep v3.50 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info Fl0:Main menu 
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Now we are ready to display the results. Choosing Results/Graphs gives us a menu of choices: 


File Input Results Options Help 
A:NHWK6P3.SEP 2 seconds 5 iterations Chemsep v3.50 


Liquid phase composition profiles 
Vapour phase composition profiles 
K-value profiles 

Temperature profile 

Pressure profile 

Flow profiles 

Mass transfer rates 

Stripping factors 

Key ratio profiles 

Relative volatilities 
McCabe-Thiele diagram 

Ponchon Savarit 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


PLOTTING THE MCCABE-THIELE DIAGRAM 


Lets first check out the McCabe-Thiele 
diagram: 


cCabe-Thiele diagram for Acetone — Methanol 


The diagram at right seems to resemble that 1-7 | | 
obtained in Hwkkey #6, except that the 
azeotrope occurs at a mole fraction of about 0.76 
instead of 0.80. This difference might occur 
because the thermodynamic model we chose is 
not the best one. In any case, it means that the 
distillate composition will never reach the 0.78 
required in Prob. 3. 


a 


im! 


a 


Let's take a closer look at the profile. To get a 
better look at what's happening at the top of the 
column, we might want to expand that section of 
the McCabe-Thiele diagram. The scale as well 
as the format of any graph in ChemSep can be 
customized by hitting the SPACE key when the 
graph is displayed. The following dialog box 
appears: 


a 


Y Acetone (Acetone+Methanol > 


oo# l l l l 
D.o 0.2 0.4 0.6 0.8 1.0 


X Acetone” (Acetone+Methanol} 
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File Input 


lst Key component 
2nd Key component 
Flip 
Minimum/Maximum X 
X Tic interval 

X Label 
Minimum/Maximum Y 
Y Tic interval 

Y Label 

Grid 

Axis color 

Stages lines 
Equilibrium Line 
Operating lines 

Q lines 

Printer Setup 
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Results 
A:\HWK6P3.SEP 2 seconds 5 iterations 


Options 


Help 
Chemsep v3.50 


McCabe-Thiele diagram 
Acetone 

Methanol 

NO 

0 1 

0.2 0 

X Acetone/(Acetone+Me 
0 1 

0.2 0 

Y Acetone/(Acetone+Me 
Off 

lGray 

1-Blue 

1-Red 

Yel low 


+ + 
Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 
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I have shaded the two lines which can be edited to expand the scale. Changing the range on both 
x and y from the interval (0, 1) to the interval (0.6, 0.8) and then hitting the ESC key yields the 
second McCabe-Thiele diagram below at the left. A second diagram appears below at right in 
which we have expanded the lower end of the diagram. 


cCabe-Thiele diagram for Acetone - Methanol 


cCabe-Thiele diagram for Acetone - Methanol 


0.20 i 


a 


Y Acetone tAcetone+hethanol t? 
a 


a 


Y Acetone/ (Acetone+Methanal } 
a a 


0.60 | l | 
0.60 0.65 0.70 


X &cetone/taácetone«Methanol? 


0.75 Be 


n.nn | l l l 
d o.02 o.04 0.06 


x Acetone tAcetone+Methanol) 


0.08 ü. 


10 
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PRINTING THE GRAPHS 


Once you get the desired graph on the screen, you probably want to make a hard copy. To 
print, hit the SPACE key to bring up the graph customization dialog box shown above. Notice 
the "Print" and "Printer Setup" items on the bottom line of the dialog. To print, you move the 
cursor over the "Print" item and hit the ENTER key. ChemSep will send the graph to the printer 
identified by the "Printer Setup" menu. 


To get the graphs into this document, I instructed ChemSep to write a Windows bitmap into a 
file which I later imported into MS Word. I edited the "Printer Setup" (by moving the cursor 
over this item and hitting ENTER) to read as follows: 


File Input Results Options Help 
A:\HWK6P3.SEP 2 seconds 5 iterations Chemsep v3.50 


McCabe-Thiele diagram 
lst Key component Acetone 
2nd Key component Methanol 
Flip NO 
Minimum/Maximum X O 0.1 
X Tic interval 2.0E-02 0 
X Label X Acetone/(Acetone+Me 
Minimum/Maximum Y O :2 
Y Tic interval 5.0E-02 0 
Y Label Y Acetone/(Acetone+Me 


Device 

Port A:NTEST.BMP 
Mode Half 

Work path 


+ + 
Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


After you have Output Device Setup reading like that shown above, hit ESC to return to the 
graph customization dialog, select Print and hit ENTER. You will be asked to supply a 
filename. Use a different name for each graph. 


When you "Quit" the ChemSep program, you can run Word. Bitmap files can be imported 
into Word using the Insert/Picture command. The picture can be cropped (eliminate unwanted 
portions of the picture around edges) and resized using the Format/Picture command. You can 
also use the mouse to do this. 


e In Word6 (Office 95), "Cropping" is done with the mouse by dragging one of the six 
"handles" while holding down the SHIFT key. "Scaling" is done with the mouse by dragging 
one of the handles without using any keys. If you use one of the four corner handles for 
scaling, the ratio of height to width will be preserved. 


+ 
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e In Word7 (Office 97), "Cropping" is done with the mouse using the “Crop Tool" which is 
shown by the icon at right. The Crop Tool can usually be found on the Picture toolbar, which 
can be displayed with the View/Toolbar command. 


Finally, you can position the graph anywhere on the page: 


e In Word6, you can enclose the picture in a frame and the move the frame. Make sure the 
picture is selected (a single click of the mouse) and the six handles appear around it. Use the 
"Insert Frame" command from the "Insert" menu. The program might ask if you want to 
switch to "Page View". Answer "yes". Once in pageview, you can drag the picture to the 
desired position with the mouse. You can also use the "Format Frame" command to position 
the picture and to control "text wrapping" etc. 


e Word7 uses Textboxes instead of frames. When inserting the picture from file, check the 
“Float over text" option of the Insert Picture dialog box. This puts the picture inside a 
Textbox which can be dragged anywhere on the page and text will wrap around it. If you 
have already imported the picture and did not check this box, you can still do so by selecting 
the picture with a single click, then use the Format/Object command. In the Format Object 
dialog, click on the Position tab and check the “Float over text" option. 


You can add text, equations, drawings, and more bitmaps. Once you have the document looking 
the way you want, print it using the "File/Print" command. 


TESTING THE MCCABE-THIELE METHOD 


Now that we have some idea how to run ChemSep, let's use it to test the validity of the 
assumptions of the McCabe-Thiele method. The main assumption is "equimolal overflow." We 
showed that if 1) you neglect sensible heat changes compared to latent heat and 2) you assume 
that the latent heat of vaporization is independent of composition, then the molar flowrates of 
liquid and vapor streams do not vary from stage to stage. 


Profiles of liquid and vapor flowrates can be plotted using the Results/Graphs/Flow 
Profiles command. For the description above, we obtain the graph below on the left: 
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Flows (mol/s} 


Flows (mol/s} 


OLiquid flow Vapour flow 


OLiquid flow Vapour flow 


The discontinuity of flowrates on the 17” stage occurs because this is the feed tray. Because the 
feed is 30% liquid in the bottom cascade is higher than in the top; because the feed is 70% vapor 
flowrate in the top cascade is higher than in the bottom. There is also a drop in both liquid and 
vapor flows at the top of the column. This drop occurs because of the subcooled reflux. 


Beyond the expected changes in flowrates at the feed tray and at the top, there is a continuous 
drop in flowrate of both streams as we move down either the rectifying cascade or the stripping 


cascade. These changes in flowrates within a cascade are not consistent with the equimolal 


overflow assumption. 


The second graph on the right above is what would be predicted by making the equimolal 
overflow assumption. We forced Chemsep to make the equimolal overflow assumption by 
altering slightly the thermodynamic model. Navigating back to the 


Input/Properties/Thermo-dynamic models menu, we change the Enthalpy setting to None: 
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File Input Results 


A:\HWK6P3.SEP 5 iterations 


Components 
Operation 
Properties 


Activity coefficient 
Equation of State 
Vapour pressure 
Enthalpy 
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Options Help 
ChemSep v3.50 


Fl:Help F2:Last help F3:Load F4:Save F5:Units F6:Options F9:Info F10:Main menu 


Selecting Return and re-Solveing the model leads to the flowrate profiles shown on the right 


side. 


These variations in flowrate within a cascade 
stem from the fact that the enthalpy of 
saturated liquid is not independent of 
temperature or acetone content as is required 
for "equimolal overflow." The plot at right 
shows the liquid and vapor enthalpies, and the 
mole fraction of acetone in the liquid, versus 
stage number.* 


For equimolal overflow, the enthalpy profiles 
should be two vertical lines. Clearly, they are 
not vertical lines. The departure arises because 
the heat of vaporization of acetone and 
methanol are not quite equal. The following 
values were read from the ChemSep library: 


acetone = 29.5 kJ/mol 


Amethanol = 35.3 kJ/mol 


Enthalpies (kJ/mol? 
-4ü -—3ü -zü -10 ü 10 
m ] T 
di J 
sL pE 
] pu 
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E dH 
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mi j J; 
[a] 
m BE 
nis f 
i 4 
20 
Oo pa 
ch | | | l 
0.0 ü.z 0.4 0.6 0.8 1.ü0 


Liquid composition (-) 


O Acetone OLiquid enthalpyt Vapour enthalpy 


* To obtain these plots, I customized the "liquid phase composition profile" graph -- deleting the methanol 
composition, and adding the enthalpies using the “top axis" for enthalpy. 
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How important are these enthalpy differences? On the McCabe-Thiele diagram, these 
differences show up as curvature in the operating lines. 


The plot at right shows points taken from — os; 
the operating lines in the rigorous 
ChemSep separation and compares them 97+ D 
with the ROL and SOL calculated using -——— 
the McCabe-Thiele method.* These  ?5| [—Eauilib 
calculated operating lines are the straight, 


solid lines on the plot. 


o 
a 


o 
p. 


There is very little difference. A closer 
examination reveals that the maximum 
difference in mole fraction is 0.00473, 
which occurs on stage 19. Put another 
way, the main prediction for this problem 
is what the distillate composition and 
flowrate will be 


Acetone in Vapor 


o 
w 


02 + 


0.1 + 


Enthalpy Model Xp D/F 0 i is å ED ne ae " di 


Excess 0.7161 0.3399 Acetone in Liquid 


None 0.7142 | 0.3408 


The difference between the two is 
insignificant. This insignificant difference is typically the case when the molar heats of 
vaporization differ by less than 20%. 


* This graph was produced by EXCEL after importing the WK1 file produced by "Results Spreadsheet" menu of 
ChemSep. 
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THE TEMPERATURE PROFILE 


Another interesting profile which can be 
plotted by ChemSep is temperature. The plot 
at right is obtained for the current problem. 
Notice that the hotter temperature occur at the 
bottom of the column. This is because the 
bottoms contain the less volatile component 
which has the higher boiling point. 


Stage number 


zür 


| | | l | 
56 58 5üu 62 64 


DEGREES OF FREEDOM Liquid temperature É C) 


OLiquid temperature 


ChemSep lets you specify a number of 
attributes of the distillate and bottoms streams. We saw this list when we specified the condenser 
and the reboiler. In particular, there were 10 different kinds of specifications which can be made 
regarding the condenser, product stream or reflux (see page 127]: 


1. Reflux ratio 

2. Heat duty of condenser 

3. Temperature of condensate 

4. Distillate flow rate 

5. Reflux flow rate 

6. Component flow 

7. Mole fraction of a component 

8. Component recovery (-) 

9. Fraction of combined feeds recovered (-) 
10. Split between two components 


There are also 10 different specifications which can also be made about the reboiler (see page 
128). 


Boilup ratio 

Heat duty of reboiler 
Temperature of reboiler 
Bottom product flow rate 
Reboiled vapour flow 
Component flow 

Mole fraction of a component 


TTO AR p 


06-361 page 138 Spring, 2001 


8. Component recovery 
9. Fraction of combined feeds recovered 
10. Split between two components 


How many of these are independent? In ChemSep, you must specify the number of stages in 
each section, the pressures everywhere, and a complete description of all feeds. Then it turns out 
that there are really only two degrees of freedom left: there are only two additional aspects of the 
solution which you can specify. 


We conclude by comparing the list of inputs and outputs for the two techniques we have 
examined. The comparison assumes we have ideal stages, constant pressure, a total consender 
producing a saturated liquid product, and a partial reboiler. 


Solution by McCabe-Thiele Graphical Method 


Given: 

l. x D 

2. desired xp 

3. R>Rmnin 

4. Nç or use of optimum feed 
Find: 

1. Ng 

2. (actual xp) 


Solution by Chemsep 


Given: 
1. XD 
2. Xp 
3. Ng 
4. Ng 
Find: 
1. R 
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Chapter 14. Gas Absorption/Stripping gas A.B 


gas absorption -- 2 components of a gas are separated by contact with liq C 


y gas abs 


a liquid (in which one component is preferentially soluble) 


stripping -- 2 components of a liquid are separated by contact with a 


gas 


An example of gas absorption is the removal of 
ammonia from air by contact with liquid water. 
Ammonia is very soluble in water whereas air is only 
slightly soluble. 


Both gas absorption and stripping involve at least 
three components. Usually only one of these 
components crosses the phase boundary. In the 
example of ammonia and air, ammonia is the 
component whose molar flowrate changes by the 
largest percentage of the inlet value. Although some 
air will also dissolve in water, and some water will 
evaporate into the air, the molar flowrates of air and 
water change by negligible fractions: their flows can 
usually be considered constant. 


These are the two main differences between them and 
distillation: 1) at least 3 components, and 2) often 
only "one transferrable component". By constrast, in 
distillation, all of the components are present in both 
phases. 


EQUIPMENT FOR ABSORPTION/STRIPPING 


Although liquid and gas streams for absorption or 
stripping could be contacted using a tray column (like 
that used in distillation), tray columns are seldom 
used. The reason is that tray efficiencies are generally 
much lower for absorption and stripping than for 
distillation (perhaps only 5% instead of 50%). 


Because of this very low efficiency, very large 
numbers of trays are required -- perhaps 100's or 
1000's. Fabrication costs just become prohibitively 
expensive. Fortunately, a viable alternative exists: 
the packed tower. 


gas C 


stripping 
liq A,B 


= out 


_~ Liquid distributor 


— Packing restrainer 


Shell 


Random packing 


Liquid 
re-distributor 


Packing support 


*— — Gas in 


— Liquid out 


06-361 page 140 Spring, 2001 


A packed tower is simply a tube or pipe, which is filled with some sort of "packing." The 
packing typically consists of particles around an inch in diameter. In commercial packed towers, 
the usual choice are particles with one of three different shapes: 


BIBIBI- E 


ia) Raschig (4) Lessing {c} Partition (e) Berl saddle (el Intalox saddle 
ring ring ring (Courtesy, Maurice (Courtesy, United States 
A. Knight) Stoneware Co.) 


(£1 Teilerette ig) Pall ring 
(Courtesy, Ceilcote (Courtesy, United States 
Company, Inc.) Stoneware Co.} 


e raschig ring (which is just a piece of pipe which has been cut into segments, whose length 
and diameter are about the same) 


2Dnal do 
Le Das to 15 inches 


e Berl saddle 
e Pall ring 


Although, in a pinch, almost anything you have laying around would do -- ping-pong balls, golf 
balls, etc. The purpose of the packing is to promote good contact between the liquid and vapor 
streams which are being brought together to permit interfacial mass transfer. 


The liquid stream is usually fed into the top of the tower while the vapor is fed into the bottom. 
Thus we have countercurrent flow of the two streams, which has the same advantages for mass 
transfer as it did for heat transfer. 


The packing promotes good contact between the phases by dividing the two feed streams into 
many parallel interconnected paths. Ideally, you would like the liquid to flow downward as a 
thin film over the surface of the packing. This would give the maximum surface area of contact 
between the gas and liquid. 


If you just pour the liquid from the end of the pipe onto the top of the packing in tower having 
much larger diameter than the pipe, most of the packing will not even be wet. Only some of the 
channels will be carrying flow. This is called: 
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channeling -- maldistribution of liquid flow 


Liquid in 
4 


T 
892 
Ow 
E 
bi 
oO 
eae 


o 
S] 


LLY 


wy Irrigated 
packing 


fatal 


Fig. 6.23 Liquid distribution and packing irrigatio: Fig. 6.14 Weir-trough liquid distributor. (The United States Stoneware 
(a) Inadequate ; (b) adequate. Company.) 3 


So some sort of device to distribute the flow over the entire cross section of the tower is needed. 
This device is called a distributor. 


Even if the flow is evenly distributed at the top of tower, channeling might still develop as the 
fluid trickles down. When two thin films converge they tend for form a thick film and a dry 
patch, which results in a reduction in contact area. So redistributors are placed every 10-15 feet 
along the length of the tower. 


A TYPICAL ABSORBER DESIGN PROBLEM 


To motivate the next few lectures, let’s pose a typical design problem. For this example, I'm 
going to take Prob. 22-1 from McCabe, 


Smith & Harriott. 
water with 


Problem: treat 500 SCFM of air containing 0.02% acetone 


14 mol% acetone to remove 95% of the L=11L 
acetone by absorption in liquid water in a 
packed bed operating at 80°F and 1 atm, 
with l-inch raschig rings. The feed water 
contains 0.02% acetone and the flowrate is 
1.1 times the minimum. 


80° F 
] atm 


500 SCFM 
air containing 
14 mol% acetone 


recover 95% 
of acetone 


The partial pressure of acetone over an 
aqueous solution at 80?F can be calculated 
from 


Pa=P4y 4X4 where ]|ny4- 195(1- x)? 


where P4 — 0.33 atm is the vapor pressure of acetone at 80°F. As the designer, you must select 
the following: 
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e flowrate of water 
e diameter of tower 
e height of packing 


Solution Overview: 


Choice of liquid flowrate: L=1.1 to 1.5 times Lmin 


Just as with distillation, there is a minimum L/V required to achieve the desired separation. In 
typical operation of an absorber, the liquid rate is chosen to be just above the minimum. 


inch H 
Tower diameter: Ap = 0.25 to 0.5 E UM 
Zr ft of packing 


The diameter of the tower is usually chosen on the basis of pressure drop. Generally, for a 
particular L and V, the smaller the column diameter, the larger the mass velocities will be, and 
the larger the pressure drop. Generally, large mass velocities are desirable because they give high 
mass transfer coefficients, but too large mass velocities cause “flooding” which severely 
decreases mass tranfer rates. The pressure drop above is about optimum. 


V/S 

Tower height: Zr= / J dy 
Kya y-y* 
—— a 
Ho, Nop 


The height of the tower is determined by mass transfer rates. Basically, the gas and liquid phases 
need to be in contact for a certain time for the acetone to diffuse from the gas phase into the 
liquid. The equation above is called the "design equation". 


vis 


y= Kya 


is called the height of one transfer unit, where Kya is the overall mass transfer coefficient times 
the interfacial area per unit volumn of packing, and 


NO J y-y* 
Ja 


is called the number of transfer units," where y is the mole fraction of the transferable 
component and y-y* is the local driving force for mass transfer. 


* This particular formula is for dilute solutions, which is probably not really applicable to the current example, but it 
is applicable in the homework problems. 
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Detailed Solution: 
First of all, what is “SCFM”? This stands for “standard cubic feet per minute” 
SCFM = standard (0°C, 1 atm) cubic feet per minute 


where “standard” means that the volume is evaluated at standard conditons. We can convert 
SCFM into moles. Starting with the ideal gas law, 


PV =nRT 
R(273°K i 3 
2: VOLUI. ( Nw liters _ ft 
n P (1 atm) gmol Ibmol 


depending of which set of units are used to express R. Recall that a gmol is the quantity of 
material whose mass equals the molecular weight in grams. Similarly, a Ibmol is the quantity of 
material whose mass equals the molecular weight in pounds. For example, air has a molecular 
weight of 29: 


MW of air = 29-8_ = 29 _© 
gmol Ibmol 


Thus we can convert 500 SCFM into a molar flowrate: 


y 


S00ft?/min | „g omol _ 
359 ft /Ibmol min 


Now let’s move on to determine the liquid flowrate required. Recall that in distillation the slope 
of the ROL is given by 


slope of ROL = 


L 
V R«1 
You should also recall that there exists a minimum value of the reflux ratio for any given 
combination of product and feed specifications. This minimum value of R translates into a 
minimum allowable L/V. There is a similar minimum L/V for gas absorption. To find the 
minimum reflux ratio in distillation, we needed to plot an operating line and an equilibrium 
curve. 


EQUILIBRIUM CURVE 


In Prob. 22-1, McCabe gives us equilibrium data in the form of partial pressure of acetone in 
the gas phase, pa, as a function of the mole fraction of acetone in the liquid, x: 
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pa = PAY 4x (109) 
P^ — vapor pressure of acetone at 80?F — 0.33 atm 


y4 = activity coefficient for liquid mixture 


Activity coefficient is a measure of nonideality of the liquid phase. For ideal solutions the 
activity coefficient is unity: 


y4 = 1 (ideal solution) 


For ideal solutions, reduces to Raoult’s law. Our solution of acetone in water is not ideal, 
but McCabe kindly gives us a model 


Iny 4 = 195(1- x? 


Note that y 4 1 (i.e. Iny4 > 0) asx > 1. This isa 
general rule: pure components also behave ideally 
since nonideal behavior results from interactions 
between different molecules; there are no different 
molecules when the fluid is pure. We can obtain an 
expression for the mole fraction of acetone in the 
gas just by dividing the partial pressure by the total 
pressure: 


yD 


where we use P = 1 atm. Repeating this for 
different x’s to obtain y’s up to 0.14 (the feed 
concentration), we obtain the curve at right, where 


x = mole fraction of acetone in the liquid (2 component is water) 


= mole fraction of acetone in the gas (2"* component is air 
y g p 


OPERATING LINE 


Suppose we are trying to absorb acetone from a mixture of air and acetone by contacting the 
air mixture with water. Let 
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L,V — total molar flowrates 


x,y — mol.frac. of transf. comp. lean gas weak liquor 


where the transferrable component is acetone in this 
example. Performing a component balance on the 
transferrable component about the top section of the 
tower yields: 


Laxa + Vy = Lx + Vyyg (jlo) 7% 
rich gas strong liquor 
Solving for y: 
L Va Jac LaXa 


=—x+ 
a V 


This relationship y(x) between the mole fractions in the liquid 
and gas streams is again known as the operating line, since the 
relationship is imposed by a component mole balance. However 
as acetone is transferred from the gas to liquid, L and V change. 
In particular, their ratio changes. Thus the operating line is not 
straight. This makes it difficult to determine the minimum 
water flowrate 


There are two limiting cases where we can easily predict the 
change in L and V. 


Approx. #1: Only One Transferable Component 


It is often possible to assume that only one component is undergoing tranfer between the 
liquid and gas streams. For example, in our problem acetone is being transferred from the air to 
the water: 


transferable: acetone 
non-transferable: air, water 


Although some of the water will evaporate when it contacts the air and some of the air will 
dissolve in the water, the molar rates of transfer of these components can often (but not always) 
be neglected compared to rate of acetone transfer. When you can neglect the molar transfer rate 
of all but one component, then significant simplification can be made. If we can neglect 
evaporation of water and dissolution of air, then the moles of water (the nontransferable 
component) in the liquid stream must be the same at all elevations: 


L' 


(l-x)L =L'=const so L= 
1-x 


(111) 
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Similarly, the moles of air in the gas stream must be the same at all elevations: 


V' 


(1-y)V s Y' = const so V= 


A L' È, L' V' 
krolanaletalinto Ere May Vy _ Lx 4 Ja 


log, 13" l-x deu 


(112) 


where L',V' are constants. If we now express concentrations in terms of mole ratios instead of 
mole fraction: 


yz- * moles of A in liquid 
= l=x moles of non-A in liquid 
y= moles of A in gas 


| ]-y . moles of non-A in gas 
then this equation becomes very simple: 


L'X, +VY=L'X+ VY, 


or y= Ex +(Q- 2X0 | (113) 


Since L'/V' = const, this is the equation of a straight line. 
Thus for the special case of one-transferable component, 
the operating line is straight on mole ratio coordinates. 
Note that mole ratios need not be smaller than unity: 


O0<x<l 


but 0<X<a 


although they are in this example. 


[1 


Lets apply this one transferable water 
component" approximation to Prob. 22-1. xg — 0.0002 
One point on the operating line is (X,,Y,). 

We are told that the inlet water contains a 

small amount of acetone: 80° F 


x,=0.0002 > X,= Te = 0.0002 BEIM: 


air containing recover 9596 
14 mol™ acetone of acetone 
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The concentration of acetone in the gas phase is determined from the specification that we want 
to remove 95% of the acetone from the feed: 


YaVa = 0.05 yp Vp 


[i [i 


pe ndr. 
1-y 1— yp 


a 


Since V’ (the molar air flowrate) is the same at either end of the column, we can cancel it out, 
leaving 


Y, =0.05 Ý%  -000814 


14 2: 
» =0.1628 
1-0.14 


The second point on the operating line must lie somewhere along the line 
Y= Y, = 0.1628 


When you plot this up, you find that the 
equilibrium curve is below the operating tr y,70.4 or Y, =0.1628_ _ 
line. This is generally true for gas I 
absorption and makes sense when you [ 
realize that the gas must be richer in the O12 + ws 
ammonia than at equilibrium — otherwise I \ 
the acetone would not spontaneously [ 
absorb into the liquid. 


X= 0.0894 
X= 0.081 


Note that there is a minimum slope which 0.04 il _ Y, = 0.0081 (0.05 Y,) 
the operating line can have before a < f 
“pinch” developes with the equilibrium 0 ! | | i 
curve. This minimum slope represents the 0 002 004 006 008 Ol 
minimum value of 
X 
(=) _ 01628- 0.00814 _ 59, 
V') min 00081-00002 ` 
Ibmol 
conio =i P aseo ERS 
min min 


Multiplying this by the molecular weight of water (18 Ib/Ibmol) and dividing by the density of 
water (8.33 Ib/gal) yields a minimum water flowrate of 5.75 gal/min. 
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Approx. #2: Dilute Solution 


In other problems, the solutions might be very dilute. If the solutions are sufficiently dilute, then 
mole ratios are virtually equal to mole fraction: 


if x<<1 

then l-x2=1 

and X= x eX 
1-x 


When this is also true for the gas stream, then (113) can be approximated as 


L' L' 
P= ht | Sa 


In fact, we can also drop the prime, since Z' ~ L: 


L L 
dni d Ya y'a 


Thus the operating line will also be linear on mole fraction coordinates if the solutions are 
sufficiently dilute. This linear operating line on mole fraction coordinates will also apply even if 
we have more than one transferrable component, provided both phases are sufficiently dilute in 
all the transferrable components. 


INTERFACIAL MASS TRANSFER: REVIEW 


An important design parameter is the depth of packing required. In distillation, the tower 
height was determined by the number of plates required times the plate separation (which is 
usually 1-2 feet). The number of plates required for a given separation is determined by the 
operating and equilibrium lines. For packed towers, the height also depends on the operating and 
equilibrium lines, but in addition it is inversely proportional to the mass transfer coefficient, 
which in distillation plays only a minor role in determining plate efficiency. 


Now that we have established the importance of interfacial mass transfer in packed towers, let's 
talk about modelling mass transfer across a phase boundary. 
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Let's start by recalling the driving force for heat transfer 
across an interface. Suppose I contact a hot gas with a 
cold liquid. The temperature profile near the interface 
will look something like that show at right. There are 
two characteristics of this sketch which are important: 


1. heat flows from high to low temperature 
2. temperature is continuous across the interface 


Interfacial mass transfer is similar to interfacial heat 
transfer, but it is also different. gas 
1. mass transfer occurs in the direction from 

high to low chemical potential (not 

necessarily from high to low Y 

concentration”) 


2. chemical potential is continuous across 
interface (concentration is generally not 
continuous) 


The main difference is evident in the sketch of the 
concentration profile near the interface. Note the discontinuity in the mole fraction at the 
interface.* 


Ji **i 
Tyj = Ty = Tj 


The reason for this discontinuity in concentration across the interface has to do with 
thermodynamic criteria for phase equilibrium. Recall: 


phase equilibrium: iy = u^ for j=1,...,N- 
thermal equilibrium: TV — TL 
y 
u 


where u; is called the chemical potential which plays 
the role of temperature in mass transfer. Unfortunately, 
there exists no "thermometer" for measuring chemical gas 
potential. Instead, we are forced to measure chemical 
concentration. While chemical potential usually 


* Within a single phase, transport is usually from high to low concentration. 


TED 
l 


* Instead of component i, the subscript will be used to denote quantities evaluated at the interface between two 
phases. In place of i, we will use j to denote components. 


06-361 page 150 Spring, 2001 


increases with concentration within any given phase (this is why diffusion of a solute occurs 
from high to low concentration), when comparing chemical potentials between two phases, there 
is no general correlation between chemical potential and concentration. 


To illustrate this, recall the simplest case of VLE: an ideal gas mixture in equilibrium with an 
ideal solution. This leads to Raoult's law: 


yP-x,Pj 
. P? 
Note that: a =} #4] 
Xj P 
Thus yj*xj for VLE 


even in the simplest case of vapor-liquid equilibrium, the mole fractions are not equal, except in 
the trivial case when you have only one component and the total pressure is the vapor pressure. 
Definitions of Transfer Coefficients 


Recall from the first part of this course, the local heat flux through the interface can be 
related to the local temperatures using any one of three types of local heat transfer coefficients: 


rate of heat transfer _ h, (Tj, - T) = hy (7; - T.) = U(Tp - Te) (114) 


interfacial area 


where hyh 


,, = one-phase local heat transfer coefficients 


U — overall local coefficient. 


Basically these equations say that the rate of heat transfer is proportional to the driving force 
(which is the departure from equilibrium) and the proportionality constant is the heat transfer 
coefficient. 


Similar definitions of transport coefficients can be made for mass transfer: 


molar rate of transfer 


: - =k, (xj -x)= ky (y - yi) = Ky(x*-x) = Ky(y—y*) (115) 
interfacial area 


where kky = one-phase local mass transfer coefficients 
KK, = overall local mass transfer coefficients 


Comparing the mass transport expressions in (115)| with their heat-transfer analogues in (114) 
there is a good deal of similarity — especially in the one-phase relations. For a single phase, the 
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driving force is the difference between the concentration or temperature in the bulk and in the 
concentration or temperature at the interface. 


But the driving force for the overall coefficients look a little different. The overall driving force 
for heat transfer is just the difference between the temperature of the hot and cold fluid 


Th -Te 


By simple analogy, you might expect the overall driving force for mass transfer to be the 
difference in concentration of the two phases: 


y - x + overall driving force 
Instead y-y*andx*-x 


appear in| (115) as the overall driving force. The *'s are 
defined as shown on the figure at right. Basically y - y* and 
x* - x are two different measures of the distance between the 
operating line and the equilibrium curve at that elevation in 
the packed bed where the liquid and gas concentrations are 


(xy). 


x-y does not represent the overall driving force for mass transfer, because the transport rate does 
not go to zero when x-y = 0. The tranport rate is zero only at equilibrium, and x-y + 0 at 
equilibrium. 


Determining the Interfacial Concentrations: (xj,yj) 


Although overall mass-transfer coefficients are the most easy to use in design, correlations 
are sometimes more available in the form of single-phase coefficients. We may have a need to 
calculate overall coefficients from the single-phase coefficients. For example, in Prob. 22-1 we 
are given the values of H, and H,. If we had instead the value of Ho, we could use the value of 
Noy we computed last time to calculate the height of packing. In what follows, we derive a 
relationship among Hoy, Hy and H,,. 


Consider once again, the rate of interfacial transport at some arbitrary 
elevation in the absorber, denoted by the dotted line in the figure at right. 
The profile of concentration of the 
tranferable component near the 
interface at this location in the 
absorbed looks like the sketch below 
at left. 


At steady state, the flux of ammonia 
through the gas film must equal the 
flux of ammonia through the liquid 


mass 
transport 
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film: 


N4 = kx) = kyO-yi) (116) 


If the fluxes were not equal, we would have ammonia building up at the interface. We will 
denote this interfacial flux by N4. 


Suppose the bulk compositions — x,y — are known, together with the two single-phase mass 
transfer coefficients. To determine the flux, we must also determine the interfacial compositions 
— Xjyj Think of these two interfacial concentrations as two unknowns. We need two 
equations. One equation is provided by the requirement that (x; y;) must lie on the equilibrium 
curve. The second relation is which can be re-written as a linear relationship between y; 
and x;: 


k k 
yj 2-—-xjtyt—-x (117) 
ky ky 
slope intercept 


The intersection gf (117) and the equilibrium 
curve gives the interfacial concentrations. 


Example: In Prob. 22-1, we are given values of 
the heights of both single-phase coefficients (Hy 
and H,) and need to evaluate the height of the overall gas-phase coefficient (Hoy). 


Solution: Consider the following problem: 


Given: ky and k, 
Find: K, 


Using the definitions of the k's and of Ky, 


N4 = kx) = k(v-vi) = Kyr-y*) (118) 
The relationship among the concentrations is —» slope = H 
shown in the figure at right. Adding and ( ky 


subtracting y; from the overall driving force y- 
y*: 


y-y*-y-yjtyj-y* (119) 
—— —— 

N4 N4 

Ky ky 


Using we can assign meaning to two of 
the three differences appear above, leaving: 
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N 4 
- +y; —y* 120 
K 5 yi-y (120) 


(x; -x)m 


The remaining difference y;-y* can be related to x;-x and the local slope of the equilibrium curve 
m: 


yj - y* - (x; - x)n 2 —m (121) 


In the second equality above, x;-x is expressed in term of the remaining single-phase mass 


transfer coefficient using Substituting into and dividing by N4: 


n (122) 
K, k, k 

ERN x 

ny — avg. slope of E.C. 


where m= 


Of course, if the equilibrium curve is straight (as it will be in dilute solutions), then m is its slope. 
This is similar to the expression for overall heat transfer coefficients for a double-pipe heat 
exchanger: 


E + E (123) 
U hy, hy 

We said that 1/U represented the total resistance to heat transfer through the two phases, which is 
just the sum of the resistances of each phase. The main difference between (122)|and (123)]is the 
appearance of m in For heat transfer, the slope of the equilibrium line is unity (m — 1) 
because at thermal equilibrium T, = m T, = T,. 


In short, (122) |states that the total resistance to mass transfer equals the sum of the resistances of 
the gas and liquid phases. We could also have showed, in a similar fashion, that 


= + (124) 
K, mk, ky 
where now: m= 5 
x* =x 


L 


If the operating line is straight, these two slopes would be the same and we could relate the two 
overall mass transfer coefficient: 
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straight E.C.: ——- 


Equimolar Counter-Diffusion vs. Diffusion through Stagnant Fluid 


In the form most analogous to Fourier's law of heat conduction, Fick's law of diffusion of a 
binary mixture of components A and B, the flux (moles/area/time) of component A in the +z 
direction is given by 


where c4 is the molar concentration of component A and D 4p is the diffusion coefficient. This 
gives the flux of A relative to a reference frame which moves with the mole-average velocity (i.e. 
the average of the species velocity of each component) of the mixture (A+B). More useful in 
calculations is the molar flux of A relative to a stationary reference frame is 


dc N4, +N 
Naz — —D Bp "n = CA (125) 


—————— 
v=mol average velocity 


where c = c 4*cg is the total molar concentration. 


Case I: If we use a capillary tube to connect two gas reservoirs having the same total pressure 
but different amounts of gases A and B (see Fig. 6.2-1), we will obtain equimolar counter- 
diffusion of A and B (i.e. N4; = -Np,). In this way, the mole-average velocity v = 0 and the total 
pressure in both reservoirs remains constant. For Ny, = -Np., [125Yfives 


N 4z =—-D4p— =J 4z (126) 


Case II: If instead, we have diffusion of benzene vapor above its liquid through air (which is 
virtually insoluble in the liquid), the air must remain stagnant (i.e. Ng; = 0) since it cannot enter 
the liquid (see Fig. 6.2-2a). This is diffusion of A through stagnant B. For Np, = 0, (125) |gives 


dc C dc 
Ng, =-D4gp—4+Ny,4 or (l-y4)N4; 2-Da4g —4 
dz c dz 
ya i 
a 
Finally we have Ny = (127) 
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For the same driving force (concentration gradient), J T is the same for both cases, but Ny, are 
different. Since l-y4 is always less than one, we see that equimolar counter-diffusion [126)]is 
slower than diffusion through a stagnant fluid [a27]. This can be qualitatively understood as 
follows. Suppose that to get to class, you need to walk down a corridor that's crowded with other 
students. If everyone else was standing almost still (stagnant fluid), it would be easier to walk 
around them than if everyone is walking toward you (counter diffusion). 


All of the flux expressions above apply locally at every point in the fluid. For 1-D stead 
equimolar counter-diffusion (from a reservoir having c4, to a second reservoir having c 45), 
integrates to 


D D , 
Naz = J4, = " = (ca 7€42)7 E 2S (Xa 742) * (vai 7 42) (128) 


where z5-z, is the length of the capillary through which diffusion is occuring and &', is the 
single-phase (gas) mass transfer coefficient to be used with mole fractions and equimolar-counter 
diffusion. The analogous expression for diffusion through stagnant film of B is (see Section 
6.2C): 


Npe CP 4B =k E (129) 
dace av (Xa 742) ^j (vai - Y42) 
(L4). -(L- 4); 
h 1- Z 
where ( yaly (1-74), 
In 
(1- y4), 


is the log-mean of yp evaluated at either end of the diffusion path. Comparing (128) |and 
we see that 
ky 


es (1=74)i 


y 


Similar relations exist for diffusion in the liquid phase: 
ky 


Nee (basa) 


x 


and between the overall mass-transfer coefficients for equimolar counter-diffusion and diffusion 
through a stagnant film: 
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M MEE 


where (I- y4),, = 


EE) (1- xar.) 


and where y4g is the mole fraction in the bulk of the gas, y4 is the mole fraction which would 


be in equilibrium with the bulk liquid having a mole fraction of x4z, and xy is mole fraction in 
the bulk of the liquid which would be in equilibrium with the bulk gas having a mole fraction of 
YAG: 


HEIGHT OF A PACKED TOWER 


The analysis which follows has as its goal the determination of the height of packing 
required. The approach is similar to that used in the design of double-pipe heat exchangers in 
which the goal is the determination of the area of heat-exchange surface required. The main 
complication arises from the fact that the mass flux is not proportional to y-x; instead, the mass 
flux is proportional to y-y*. 

We have to chop up the tower into pieces which are small E, L, 

enough so that the driving force is virtually uniform throughout Xa 
each piece. Since the compositions change only with z, we chop 
up the tower in such a way that we produce pieces which have 
z=const, which is a thin horizontal slice. 


Now let's take a closer look at what happens inside this slice of 
the tower. The slice contains solid packing as well as liquid and 
gas streams. In what follows, we will ignore the solid and treat Y, X; 
the liquid and gas phases as if they were completely separated, 

rather than interspersed in each other. 


double-pipe hxer analogy: 


Let's do a mass balance on the acetone in the gas phase_only 
contained within our slice of column. Besides the liquid streams 
entering and leaving the slice, we have acetone crossing the 
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interface between the gas and liquid phases. ^ The rate of 
absorption per unit area can be expressed in terms of the local 
overall mass transfer coefficient: 


molar rate of transfer Hu 


K — y* 130 
interfacial area y o») mo 


We could also have used one of the other expressions in (115) 
this one proves later to be more convenient than some of the 
others. 


Just like the overall heat transfer coefficient U, depends on the 
flowrates of both fluids being contacted in the heat exchanger, 
the overall mass transfer coefficient K, depends on both 
flowrates 


Rd LV) 


y 
To obtain the rate of transfer across the interface in our slice of 
column, we need to multiply by the interfacial area in that 
slice. 


In a heat exchanger, the heat transfer surface is fixed by the 
geometry of the equipment selected: it is just the area of the pipe 
wall or the tubes. In particular, the heat transfer area does not 
depend on the flowrates of the hot and cold streams. On the 
other hand, the boundary between liquid and gas in a packed bed 
is very complex and very hard to measure directly. Most 
importantly, the area also depends on the flowrates of the gas 
and liquid streams (L, V). 


The interfacial area is usually expressed as a, the interfacial area 
per unit volume of packing: 


interfacial area 


—=a=a(L,V) (131) 
volume of packing 


One empirical correlation relating area to flowrates is the 
Schulman Correlation (Table 6.3 of Treybal): 


a= mG? Gh 


where G, and G, are the mass flowrates of the liquid or gas 
stream divided by the cross-sectional area of the tower 


Spring, 2001 


Uo - Te) 


nD,dz _ 4 


n(D,/2))dz Po 
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mass/ time 
nD? 
4 


and where m, o, B are constants which depend on the type and 


size of packing used. Multiplying by 


rate of transfer interf area 
x 


G= = mass velocity 


= K,(y-y*)xa 
interf area vol. of packing y=") 
molar rate of transfer 4U 
= - = Kya(y - y*) (132) (T, — T.) 
volume of packing D, 
D; 
Me » 


faz 


If we now multiply by the volume of this slice of the packed 

column, we will obtain the rate at which acetone crosses the 

interface in this slice: 

dq =U, dA, (7j - T.) 
—— 


rate of transfer = r S Az = dN, xD, dz 
2 
: TD 
where S — empty tower cross-section — T 


S Az = volume of slice 


At steady state, the rate of transport of acetone into the vapor 
must equal the rate out: 


in = out 
| | T MPC py T(z)= 
Vy) -= (V + rSAz 
(Vy) z 7 (Vy) a4 (133) ip, Ts + Az) 4 dq 
Dividing by -Az and letting Az > 0, (133)]|becomes: 
mpc em = ud 
d(Vy) = -rS (134) P dz dz 


dz = nU, D, (T, - Te) 
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Now the total molar flowrate (V) will change as the gas phase loses acetone, but the flowrate of 
acetone-free air (V’ ) does not change with Z: 


y = FEE pa 
l-y l-y 
(9) - P2) reg a J- =" Way ay 
Ie ty (1- y) 1-y 
» | V dy 
becomes: = -rS = -K -y*)S 135 
i l- y dz , yal y y") d 


The second equation is obtained by substituting our expression for the rate of absorption r from 
(132)] We can now solve this equation for the height dZ of the slice: 
—Vdy 


di (I= »» - » *(K,a)s 


So if we know the local mole fractions y and y* and the change in the mole fraction dy of the gas 
phase which occurred between the top and bottom of this slice, we can calculate the height of the 
slice. The total height of packing is the sum of the height of each slice: 


e jp (1- y TERIS a)S 


Factoring out the S, which is constant along the entire height of the packing, and changing the 
order of integration (since y,>y,): 


NS S Vdy 
7s J ya 07 xv 7» *)(Kya) 


Now K, depends on G, and G,, which in turn depend on gas composition. For the case of a 
dilute gas stream, we know that 


y<<l l-y#l V «V' = const Ka x const 


After making these approximations, we have 


Yb 
Ties Í dy (136) 
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which serves as the basic design equation for packed bed absorbers. 


For very dilute solutions, this integral can be evaluated analytically (1.e. numerical integration 
is not required). If both x««1 and y<<1 hold, then the molar flowrates of the two streams won't 
change much as the solute is transferred across the phase boundary. In particular, the ratio of the 
flowrates will be constant: 


then L/V ~L'/V'= const 
This means that the operating line will be straight, even on mole fraction coordinates: 
y= (Vy + [ya - (L/V)xa] 
Also, for sufficiently dilute solutions, Henry's law always applies so that 
y*=mx 


will also be a straight line. When both operating and equilibrium curves are straight, then the 
integration can be done analytically instead of numerically. The result is: 


straight OL & EC: Nj e db d. (137) 


where ( y-y *) z 18 the log-mean of yp-yp* and Ya-Ya” which represent the driving force at the 


bottom and top of the tower, respectively. 


TRANSFER UNIT 


Meaning can be given to the integral. Even if the integrand varies significantly over the 
domain of integration, we can still pull it outside the integral, provided we replace it with some 
appropriate value, which can be thought of as the mean value for this range. Recall the Mean 
Value Theorem of Calculus: 


Joy cu n I eg VOR ue 


Yb 1 Yb z 


Of course, it is not usually known what type of average to use, but for the present purpose, the 
precise value is not important. The integral can be thought of as the total change in gas mole 
fraction divided by the average driving force. 


(Gas Phase) Transfer Unit -- a slice (not necessarily thin) of an absorber in which the gas 


undergoes a change in y equal to (y — y gee which represents the average driving force 


over the entire absorber. 
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No. of TU's -- the number of these it takes to make the entire absorber (not necessarily an 
integer) 


Height of TU -- height of a slice of an absorber corresponding to one TU 


From we see that the integral clearly represents the number of transfer units: 


Yb 
Noy | dy : (139) 
Mg 4 — V 


This is called the number of overall gas phase transfer units since it is based on the overall 
driving force expressed as the gas-phase mole fraction. Since the product in (136)|represents the 
total height of packing, the coefficient of this integral must represent the height of a transfer 


unit: 
V/S 
Hp, =| —— 
id ej 
avg 


So that our design equation can be rewritten as: 


Zr = HoyNoy 


which is called the HTU method for sizing an absorber. Throughout the derivation above we 
made use of one type of mass transfer coefficient: 


r= K,a(y-y*) 


Of course, other types of mass transfer coefficient are commonly encountered and can also be 
used to compute the height of packing required: 


r= Kya(x*-x) = k,(x-x) = ky(y-yj) 


to mention a few. Of course, there are still more definitions based on driving forces expressed as 
differences in molar concentration or partial pressure. Any one of these mass transfer 
coefficients can be used in the HTU method: 


ZT = HoyNoy E Ho Nox = H,N,— ALN, 


where 
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Xa dx L/S 
Nox -Í 7 Hox E ) 
x, x *—x ad sn 
Xa L/S 
N, = dx H, -(4| 
x, Xj—X xs 
Yb V/S 
N= VY H= VIS 
Va X Vi kya avg 


Note that the N's are not equal to each other; thus the number of transfer units depends on which 
driving force you use. 


Ny ZN, Oy 
, = water 
EVALUATION OF NTU'S Ya = 0.00807 x, = 0.0002 
Let's see how this works by (95% recovered c: 


continuing the example we started a in liquid) 

couple of lectures back (see page 

141). The problem is summarized in 

the figure at right. We want to mE 

evaluate Noy fro | (139). The y in ice Ee 

this expression is the mole fraction of ,, Je 
V' — 1.39 Ibmol/min 

acetone in the absorber at some 

particular elevation. If the 

corresponding mole fraction of 

acetone in the water is x, then two are related by the operating line: 


y): 7j. RE UNES. SINN RT ANNE (140) 
lex qx l-y 1—yg 


To obtain x, = 0.07 (specified in problem statement) with x, = 0.0002, y, = 0.00807 (5% of 
acetone remains with gas) and y; — 0.14, we find that 


n = 2.06 


With this ratio known, given a value of y in the interval y, < y < yp, we can use the operating line 
to solve for x(y). On the other hand, y *(x) is the mole fraction of acetone in the gas, which would 
be in equilibrium with a liquid having mole fraction x. 


x, — 0.07 


y*(x) = Pace Hey x = 033exp|1.95(1 - x)? |: (141) 
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A plot of y(x) and y*(x) is shown below at left. Neither curve is a straight line on these 


coordinates, although the operating line is close. 


To evaluate the integral, we need to compute the integrand at various points in the domain: 


1 


y xfrom(140)| y* from (141) 


y-y* 

Ya = 8.07E-03 2.00E-04 4.64E-04 131.401 
0.022 7.18E-03 0.016 168.635 
0.036 0.014 0.031 204.747 
0.05 0.021 0.045 225.474 
0.063 0.028 0.059 219.902 
0.076 0.035 0.000 | 192.77 
0.089 0.042 0.083 158.147 
0.102 0.049 0.094 126.193 
0.115 0.056 0.105 100.38 
0.128 0.063 0.115 80.583 

yp = 0.14 0.07 0.125 65.615 


A plot of the integrand might look something like that shown above at right. The area under this 
curve represents Noy. 


No 


Yb dy 
y xil = 210 


y,y-y* 


To evaluate this integral in Mathcad, it convenient to pose y and y* as functions of x: y(x) is the 
operating line given by (140) and y*(x) is the equilibrium curve given by (141)] Then 
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d 
dy 2 dx=y'(x)dx 


-{” dy dx _ [m ano 
X 


N = 
n Ya Jet a Yæ) -y *(x) 


yY 


The maximum value of the integrand occurred at the point in the column where the operating line 
comes closest to the equilibrium curve. The denominator y-y* is a measure of this distance. 
Because the equilibrium curve turns downward in this case, the "pinch" point occurs in the 
middle of the range. In other cases (e.g. NH3 in water), the equilibrium curve turns upward 
instead. 


Special Case: y4««1 with concave downward EC. 


The figure at right shows a similar plot you 
will obtain for Prob. 22.2 on the Homework. 
In this case, the denominator of the integrand 
becomes very small near the lower limit; thus 
the integral has much of its area under a 
"spike" located at the lower limit. 


To integrate such a function accurately, you 
will need a very small grid spacing — which 
means a lot of points — or at least very small 
spacing under the spike — which means a 
nonuniform grid. One solution is to transform y 
the integral by multiplying and dividing the 

integrand by y: 


y y y 
No zi zi =|” 7 2f” Y dny 
Y yk ES 

Ya yTy Ya YTY y 


Thus we plot y/(y-y*) versus lny. Since both 
the numerator and the denominator become 
small at the lower limit, the quotient tends to 
be more near constant with y, allowing the 
integration to be done with fewer points and a 
uniformly spaced grid. 


y - yeq 


ANALOGY WITH DOUBLE-PIPE HXER 


Comparing (137) with (138), we see that 


for very dilute solutions, the average driving 
force is the log-mean. This result is quite In(y) 
analogous to our design equation for double- 
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pipe heat exchangers. To see this analogy, let's compute the total rate of ammonia transfer over 
the entire packed tower: 


qT = V(ys-ya) (142) 


V/S Yp-Ya 
Kya (yv-y*), 


Recall the design equation: Zr = Ho, No, = (143) 


This contains the product that appears in (142)| Solving|(143) |for Vivp-y,) 


V(yp -y»4)* Kia SZr (y-y*), 
total volume 
of packing 


where SZrT is the total volume of packing in the absorber and a is the interfacial area per unit 
volume of packing. Thus 


aSZry = total interfacial area in entire absorber = AT 


and substituting in (142) yields: 
qr- KyAr(y - » *) 
which is analogous to: qr = UATAT], 


Thus, for dilute solutions, the design equation for an absorber is identical in form to that for heat 
exchangers: with the overall mass transfer coefficient K,, analogous to the overall heat transfer 
coefficient U and the log-mean of y-y* at the two ends of the absorber analogous to the log-mean 
AT. 


RELATIONSHIPS AMONG Hox, Hoy, Hx AND Hy 


In on page we obtained a relationship between the overall mass-transfer 
coefficient and the single-phase mass-transfer coefficients. This relationship can be easily 
modified to relate the corresponding heights of transfer units. 


Dividing both sides of (122) by a (the interfacial area per volume of packing), we obtain: 
1 1 m 


= + 
Kya kya kya 


(144) 


where m = (y-y*)/(x;-x) = avg. slope of E.C. 
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Of course, if the equilibrium curve is straight (as it will be in dilute solutions), then m is its slope. 
This is similar to the expression for overall heat transfer coefficients for a double-pipe heat 
exchanger: 


LS + = (145) 
U h, h, 

We said that 1/U represented the total resistance to heat transfer through the two phases, which is 
just the sum of the resistances of each phase. The main difference between (144)|and is the 


appearance of m in For heat transfer, the slope of the equilibrium line is unity (m — 1) 
because at thermal equilibrium 7), = m T, = Ty. 


In short, states that the total resistance to mass transfer equals the sum of the resistances of 
the gas and liquid phases. Using we could also have showed, in a similar fashion, that 


1 1 1 


= + 
K,a  mkya ka 


where now: m = (y-y;)/(x*-x;) 


If the operating line is straight, these two slopes would be the same and we could relate the two 
overall mass transfer coefficient: 


straight E.C.: 1/K, = m/K, 


If we now multiply by V/S we have the height on an overall gas transfer unit on the left- 
hand side: 


V/S _V/S mV/S VÍS V L/S 


Ka k ka k Lk 
ya ya xa yd xa 


Hoy H, Hy 


Recalling the definitions for various heights of a transfer unit we end up with 


4 
Hoy = Hy +m- H; 
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1077 —74 71 Zi 4L 


i — VE dM eee 
“LZ Z 


Z I 
/ f | Column dia. = 30 in. 
A A VFI Packing height = 10 ft 


AP, INCHES WATER/FT PACKING 


Liquid rate, lb/ft? - h, 
as parameter 


100 200 300 400 600 1,000 2,000 5 000 
AIR MASS VELOCITY, Gy, Ib/ft? - h 


Fig. 8. Pressure drop across a packed bed of 1-inch ceramic Intalox saddles as a function of 
air and water flowrates. Taken from Fig. 22-4 of MS&H5. 


PRESSURE DROP IN PACK BEDS 


Usually, flow in a packed-bed absorber is countercurrent. 
Gravity causes the liquid to flow down through the packing, 
whereas a small pressure drop drives the gas flow upward. This 
pressure drop plays an important role in packed beds. Without 
liquid present, the pressure-drop across dry packing increases 
approximately with the square of the gas flowrate (see Fig. 8]. 


packing 


Ib / hr of gas E My,V 
cross - section of tower S 


G, = mass velocity of gas = 
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where M, is the average molecular weight of the floodin 
gas. Once liquid is flowing down the column as log(Ap) 
gas is flowing up, some of the space between the 


ME Lu loading 
packing particles is taken up by liquid, leaving less x " 

* & 1.8 to 2.0 
space for the gas to flow. Forcing the same gas constant / (turbulent) 
flow through a smaller opening will increase the holdup 
frictional contribution to the pressure drop. Thus r$ 
the curves for increasing liquid flow are above that wet (G,>0) dry (G,=0) 


for dry packing although the curves tend to be 
parallel. log(G,) 


Note that the curves with liquid flow curl up at 
high gas flow rates. This can 
be explained as follows: 


The upward flow of the gas 
exerts a shear force on the 
liquid, retarding its downward 
motion. As the gas flowrate 
increases so does this shear 
force. When the shear force 
becomes comparable to gravity, 
the liquid flow might slow 
down and liquid begins to 
accumulate in the tower. 


packing 


stagnant packing 
particle i 


flowing 
air particle i 


air 


The amount of liquid residing or accumulating in the tower is called: 
(liquid) holdup -- fraction of interstitial volume occupied by liquid 


intersticial volume -- "empty" space between and inside packing particles; that volume 
occupied by liquid or gas. 


loading -- an increase in liquid holdup caused by an increase in gas flowrate 


For a given liquid flowrate, there is a maximum flowrate of gas which can be forced through the 
column. If you exceed this maximum the shear forces on the liquid are so high that they exceed 
the weight of the liquid. Then the net force on the liquid is upward and you blow the liquid back 
out the top of the column. This is called: 


flooding -- liquid downflow is essentially stopped by high gas upflow 


Since you can't get the liquid through the column, flooding is clearly a condition to be avoided. 
Generally, to get a high area of contact between liquid and gas, you want to operate well below 
the flooding velocity: 


at flooding: ApríL « 2 to 3 inch H5O/ft of packing 
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loading: Apr/L ~ 0.5 inch H5O/ft 
normal oper.: Apr/L ~ 0.25 to 0.5 inch H20/ft 


Normal operating pressure drops are just below those for which loading begins. 


Notice (from Fig. 8) hat when we increase — *"? 


the liquid flowrate, flooding occurs at 
lower gas flows. If we increase the size of 
the packing particle, we can tolerate higher 2000 
gas flows. While higher gas flows can be 
obtained with larger packing particles, the 
amount of interfacial area is generally less 


(ie. a is decreased). These effects are 

summarized by Fig. 9at right. " E 
Fig. 8]and Fig. 9 lare for a particular 

packing (1" Intalox saddles) and for a 

particular temperature and pressure (20 °C 200 

and 1 atm). A more general correlation of EE i e" 

pressure drop is provided by Fig. 10. |The EN lim 

. . . 100 
size and type of packing is accounted for 1000 2000 5000 10,000 ^ — 20,000 50,000 
in this correlation through the parameter RELIER 


Flooding Velocity G,, Ib/ft?-h 


Fig. 9. Air flowrate at which flooding 
occurs for ceramic Intalox saddles at 
different water flowrates. Taken from Fig. 


F, = packing factor 


Values of the packing factor for various 22.5 in MS&HS. 
! types of packing are given in Fig. 11. The y-coordinate of this 
Symbol Units graph in Fig. 10 |s not ens ea i need to use the units 
Gy, Gy Ib/ft*-s summarized in the table at right. 
Hx cP 
Px Py Ib/ft? 


g, | 32.2 Ibe fi/Ib-s2 


TOWER DIAMETER 
Once the total molar liquid and gas flowrates (L and V) are known, we can choose the 
diameter of tower we need. The diameter of the tower is usually chosen such that the pressure 


drop is some prescribed value below flooding: 


choose Dy such that: E = 2 & 0.25 to 0.5 inch H,O/ft of packing 
Z T 


The pressure drop depends on the mass velocities: 


A 
AP =AG,, G,) 
Zr 
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Parameter of curves (s pressure 
drop in inches of wafer/foof 
of packed height 


es 
Ba 


H 


N 
i 
1 


0.01 0.020.040.0601 O2 0.40.6 
Gx Py 
Gy Py-Py 


Fig. 10. General correlation of pressure drop across packed beds. Taken from Fig. 22.6 of 
MS&HS. 


oe M 
hel 62 l mass flowrate of liquid | | M,L M,L 


cross - sectional area of column S Di 
4 


and where M, is the average molecular weight of the liquid. 


MV 
NES 


Similarly, 
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Packing 
factors 
Nominal Bulk Total Porosity, 
Type Material size, in. density,+ lb/f ^ area,t fÜ/fÉ £ F f 
Raschig rings Ceramic j 55 112 0.64 580 1.528 
1 42 58 0.74 155 1.368 
IE 43 37 0.73 95 10 
2 41 28 0.74 65 0.928 
Pall rings Metal 1 30 63 0.94 56 1.54 
14 24 39 0.95 40 1.36 
2 22 31 0.96 27 1.09 
Plastic 1 5.5 63 0.90 55 1.36 
14 4.8 39 0.91 40 118 
Berl saddles Ceramic 4 54 142 0.62 240 1.588 
1 45 76 0.68 110 1.368 
14 40 46 0.71 65 1.07§ 
Intalox saddles Ceramic 1 46 190 0.71 200 227 
1 42 78 0.73 92 1.54 
14 39 59 0.76 52 118 
2 38 36 0.76 40 10 
3 36 28 0.79 22 0.64 
Super Intalox Ceramic 1 = = — 60 1.54 
saddles 2 — — — 30 10 
IMTP* Metal 1 — — 0.97 4] 174 
14 — — 0.98 24 137 
2 — — 0.98 18 119 
Hy-Pak Metal 1 19 54 0.96 45 1.54 
14 — — — 29 1.36 
2 14 29 0.97 26 1.09 


t Bulk density and total area are given per unit volume of column. 
ł Factor F, is a pressure-drop factor and f, a relative mass-transfer coefficient. 
§ Based on NH,-H,O data; other factors based on CO,-NaOH data. 


Fig. 11. Characteristics of common packings. Taken from Table 22-1 of MS&HS. 


So to get a certain pressure drop for a certain set of flowrates, we must choose a particular value 
for the tower diameter. Although it looks like you have to guess the diameter Dy in order to 
calculate G, and G, to get Ap, a trial-and-error can be avoided by noting that the abscissa 
(x-coordinate) does not depend on the diameter: 


M,L 
Ge le MIL 


85 H MyV 
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Note that the unknown cross-sectional area S conveniently cancels out. What remains is known. 
So the procedure is as follows: 


Given: L, V 
Find: Dr 
Solution: 
1. calculate the abscissa 


2. locate the point on the curve of Fig. 22-6 which has the desired Ap and this value 
of the abscissa 


3. read the corresponding ordinate of this point from Fig. 22-6 
4. calculate G, which gives this value for the ordinate 
5. calculate the desired tower diameter from 


D? _ 9 My 
4 G 


y 
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Chapter 15. Membrane Separation 


In distillation, gas absorption, and extraction we have 


two immiscible phases which we bring into intimate A 
contact. à 
Operation Phases contacted A 
distillation liquid and its vapor 
gas absorption liquid and another gas 
solid extraction solid and a liquid 
liquid extraction two immiscible liquids 


Another unit operation for mass-transfer is the membrane 

separator. Here the two phases are not immiscible — usually they are both gases or both are 
aqueous solutions. They are not intimately contacted; rather, a membrane is placed between 
them. 


This membrane serves two functions: gas gas 

A,B A,B 
1. keeps phases from mixing high P low P 
This is also the primary function of the pipe wall in HXer; A 


however, the membrane has a second, even more important, 
function. Whereas the pipe wall plays a passive role in heat 


transfer, the membrane plays an active role: 
sfer, the membrane plays ctive role membris 


2. controls selectivity 
Because of differences in molecular size or solubility of the gas in the membrane material, 


different species diffuse through the membrane at different rates. This difference is the basis of 
separation. 


EXAMPLES 


With the proper choice of membrane material, a membrane separator can achieve significant 
separations in a single stage. Indeed, as we shall see, it usually turns out that multistage 
operations with membranes are not economical (compared to cryogenic distillation). Thus most 
industrial applications of membranes do not require very high purity products. 


Gases 


e separation of 238UF, from 235UF,. This was a very important part of the atom bomb 
project in World War II. 
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separation of O, from N; in air. A number of industrial processes require a chemically 
inert gas (like nitrogen). Membranes can produce 95-99% nitrogen from air (79% 
nitrogen, 21% oxygen) in a single stage. Sometimes its the oxygen-rich stream that useful: 
40-5096 oxygen can be obtained from air, which is useful for medical patients having 
respiratory conditions. 


recover H, from purge streams in ammonia, methanol, and hydrogenation plants. 
separation of He from CH, (natural gas). Because of its very small molecular size, half 


of the helium in a feed stream containing <1% He can be recovered in a single stage, 
resulting in a permeate stream contain 30 times greater concentration of He. 


Liquids 


dialysis (or liquid permeation). Separation of urea and other small molecules from protein 
and cells in whole blood is accomplished in artifical kidneys. Recovery of caustic (NaOH) 
from hemicellulose produced in the manufacture of rayon. 


electrodialysis. Brackish water (containing less salt than sea water) can be made potable 
by applying electric field across membrane. 


reverse osmosis. Brackish water can also be made potable by applying pressure to force 
water through a membrane which is impermeable to ions. 


In the 1950's, President Eisenhower was convinced that drinking water was the key to world 
peace in many parts of the world. Eisenhower's successor as President was John Kennedy. In a 
famous speech, President Kennedy expressed his intention to “land a man on the moon and make 
the deserts bloom." The second part of this quote has been obscured by the first part. 
Nonetheless, Kennedy did follow through by creating the Office of Saline Water, which provided 
financial support for much of the early research on membranes. 


EQUIPMENT 
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Permeate 
pe pipe containing . 
collection holes Antitelescoping device 
Feed solution N 
Concentrate 


wa — Permeate out 
A > Concentrate 
2 


iN 


Permeate ^" | 
Y JJ 


S 


Feed flow 
across feed 
channel spacer 


p" 


Feed solution 


Membrane » 
Permeate flow (after passage 


through membrane into permeate 


Permeate ollection material : 
collection material) 


Membrane Covering 
Feed channel spacer 


Fig. 12. A spiral-wound membrane 
Fig. 13. Magnification of a hollow-fiber 


member. Outer diameter between 50- 
200 microns. 


i Hollow, thin-walled 


LA plastic tube 


Retentate out 


Fig. 14. A bundle of hollow fibers 
resembles a shell-and-tube heat 


exchanger. 


Spring, 2001 


06-361 page 176 
Epoxy 
Open ends tube sheet 
Concentrate FI soap of fibers Porous 
Snap ring outlet ad m back-up disc 
Oring seal Snap ring 


Permeate 


End plate Oring seal 
End plate 


‘Fiber Shell Porous feed 
distributor tube 


Fig. 15. A commercial cartridge containing hollow fiber membranes. This is a Permasep™ 
cartridge made by DuPont. 


Fig. 16. Flow patterns inside bundle 


"i 
o| S s => (side views). (a) shows hollow fibers 
a E — with left ends sealed and right ends 
i = i open to exit manifold. (b) shows hollow 
| fibers with both ends open to exit 
A | 
(a) manifolds. 


Fig. 17. Flow patterns inside bundle (end view). 
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SEPARATION OF GASES 


Porous Membranes 


Perhaps the simplest basis for separation of two gases is molecular size. Gases of different 
molecular weights have different diffusion coefficients and will diffuse at different rates through 
the membrane. For membranes having pores whose diameter is comparable or smaller than the 
mean-free path of gas molecules, the diffusion coefficient is inversely proportional to the square- 
root of molecular weight: 


T 
A 


for Knudsen diffusion: D4 =9700r 


(146) 


where r = pore radius (cm) 
T = absolute temperature (K) 
M , = molecular weight of diffusing gas A 
D, = diffusion coefficient (cm2/s) 


Notice that larger molecules (i.e. larger M4) have smaller 
diffusion coefficients. 


For membranes having fine pores through which the gas cai 


can diffuse, the fluxes of two gases, A and B, through a porous 
membrane are given by C ((z) 


C4(0) - C4(/) 
P te 


Jah (147) 


Cp (0) - Cg(/) 
Jp; = Dee B B 
te 
where / is the thickness of the membrane and « (which is greater than unity) is a dimensionless 
factor called the tortuosity which accounts for the fact that the pore length is generally larger than 
the thickness of the membrane. 


Even if there are equal driving forces for the two components to diffuse across the membrane, 
they will have unequal concentrations one the downstream side because of their different 
diffusion coefficient. This is the basis for separation in a porous membrane. 
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The best known example of gas separation by a porous membrane is separation of 238UF¢ from 
235UF&. Most gases have much smaller molecules than these and polymer membranes are used 
instead of porous membranes 


Nonporous Polymer Membranes polymer 


gas 
Many polymer membranes do not have pores. Instead, 


the gas must first dissolve in the polymer. Then the 
dissolved gas diffuses through the membrane, which can be  C4(Z) 
thought of as a liquid-like film. As a result, the 
concentration profile of the diffusing component A is not 
continuous (see figure at right). 


Generally the mole fraction of the dissolved species is very 
small so that the equilibrium solubility of the gas in the 
solid can be described by Henry's law: 


C4-S4P, and Cp =SpPp (148) 


where C4 = molar concentration of A in membrane 
P, — partial pressure of A in gas 
S , — solubility constant* 


Again, at steady state the concentration profile inside the membrane will be linear and Fick’s law 
yields: 
dC 4 


CT dex 


(149) 


,£49-040. ps 


P4(0) - Pa 
4 A 


For the purpose of comparison, the flux through a porous membrane [given by|(147)]] can also be 
written in terms of the partial pressures. For an ideal gas, the molar concentration C4 can be 
related to the partial pressure of the gas P4 by 


Pyv = n4RT 
P 
C,- PA ME d 
V RT 
D,P -P 
L475 becomes: J Az = 4 PAQ)- Pa porous membrane 
RT £ 
P = P4 (2 
"m J47; - D4S 4 40 7 4(0 polymer membrane 


* reciprocal of Henry's law constant 
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Although these two expressions have the same form, the coefficients are different. In particular, 
the solubility constant $4 varies much more between various gas solutes than does the diffusion 
coefficient D4. This implies that polymer membranes tend to be much more selective in 
separation various gas species than porous membranes. Unfortunately, diffusion coefficients in 
solids and liquids are much smaller than for gases so this increase in selectivity is often traded off 
with lower permeation rates. 


The product D 4S 4 is known as the 
permeability coefficient: q4- DAS , [7] Barrer 


which has units of flux per unit pressure gradient. Since these complex units occur frequently in 
such problems, a new unit is defined, which is known as the Barrer: 


3 3 
(pea io Oo gg S oe: 
S cm? cmHg s- cm - cmHg 
R — —— 
rate l/area dr 
flux 


Often the thickness of the commercial membranes is not known. Then it 1s easier to measure the 
flux per unit pressure, rather than per unit pressure gradient. The flux per unit pressure drop is 
known as the 


_ D484 


permeability: Q4 7 


EXAMPLE: PRODUCTION OF ENRICHED AIR 


It turns out that most polymers are more f 

. eed concentrate 
permeable to O, than to N2. Let's see how this, —519, O, 
fact can be exploited to enrich the oxygen (79% N.) 
content of air. Consider the apparatus at right, 
which has one feed and two product streams. permeate 
One of the product streams has been pulled membrane 
through a polymer membrane with a vacuum — >) y 
pump. BERGE 


The feed is normal atmospheric air which 

contains 21% oxygen and 79% nitrogen. We 

will denote these as components A and B. From our analysis of steady diffusion through the 
membrane, we know the fluxes of the two species must satisfy: 


J4 5 Q4(P41 - P42) 2 QA| xB — y B & xQ4P (150) 
<<P, 
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Jp RI (1- x)Og P, 


Because we are using a vacuum pump to pull the gas through the membrane, the pressure on the 
permeate side of the membrane will be much below atmospheric, thus: 


for vacuum (on side “‘2”): Py << P| 


This leads to the approximate expressions for the fluxes (above). Since the entire permeate 
stream must come through the membrane, the mole fraction of A in the permeate stream is 
determined by the relative amounts of the two species which come through. In time ¢ the number 
of moles of species A which comes through a membrane of area A is J,At. Similarly, the number 
of moles of species B which comes through is JpAt. The mole fraction of A is the number of 
moles of A divided by the total number of moles. The area and time cancel leaving 


o adt 20 d 
d J 4At+ Jp At JA4tJpg 


(151) 


Substituting (150)|into 


e xQ4P _ xO4 
xQ4P, «(1- x)OgP,  xQ4 *(1- x)Op 


Notice that the total pressure cancels out. Finally, we divide the numerator and denominator by 
the permeability of gas B: 


OX Ox 
Yee) (Geet et 
where "EAD 
Op 


is called the separation factor. It's role is similar to relative volatility in distillation so we give it 
the same symbol. If we did not make the assumption of a vacuum on the downstream side, the 
analysis becomes a little more complicated, but still tractable. More generally, becomes 


J4 = Q4RP (x - Ry) 
Jg = Og A |[(1- x) - &(1- y)] 


where R = —- 


and becomes the root of the following quadratic equation: 


R(a -1)5? —{(a -1x +R)+1I}y+ax=0 (153) 
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Note that 1f R — 0, reduces to|(152). At right is a plot of y versus a for several different 
R’s. A couple of trends are apparent. First, 


yox as a1 


a=] means there is no difference in 
permeability of the membrane to oxygen 
or nitrogen: then the permeate has the 
same concentration as the feed. 


Second, note that the permeate tends to 
pure oxygen if the permeability ratio is 


high enough: 


for R<x: yo l as ao 


A closer look reveals that pure oxygen is 
only obtained if the pressure ratio is 
small enough. In particular, the ratio 


Permeate Mole Fraction 


must be smaller than the feed 
molefraction. 


Permeate concentration y in the limit 
as a — oo as a function of pressure 
ratio R. 


but P= RP, 


ktssylintolesa] 


Permeability Ratio, a 


The plot at left shows how lim (y) depends on the 


0, — 00 
pressure ratio R. when R<x, then y-1. When R=1 (i.e. 
no driving force for flow through membrane), the 
permeate has the same concentration as the feed. 


R=1 and ao: y=x 


This behavior can be rationalized as follows: the partial 
pressure of oxygen in the permeate can never exceed 
that of the feed: 


os ah ase 
Po, in permeate Po, in feed 


(155) 


* 
R 
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If this inequality is not met, we would have diffusion uphill. The upper bound in oxygen partial 
pressure is reached when a — o: 


l x 
lim y= mini =) 
100 R 


Stage Cut 


Now let’s look at a couple of complications. The fraction of the feed which permeates 
through the membrane is known as the stage cut: 


. total molar permeate rate Vout feed residue 
~ total molar feed rate Lg Lin +y i Wy Y Lout 
permeate 
When the stage cut is very small, then very little of Vout 


the feed passes through the membrane 
Co-current flow through separator 
Lout 8 Lin 


This also means the composition of the feed doesn't change much as a result of passing though 
the separator, so that 


stage cut «€ 1:  Xout™ Xin feed residue 
Lig Lout 
Then every point on the membrane surface sees permeate Y TE Y n 


the same feed concentration x. This is what we 
have implicitly assumed in all of the above 
analysis. Under these conditions, it doesn't 
matter whether we operate the separator 
cocurrently or countercurrently; we get the same y for both flow patterns. But this is not a very 
economical way to operate: very little of the feed is converted into product. 


Vout 


Countercurrent flow through separator 


If the stage cut is not small, we will partially deplete the feed in the more permeable component: 
oxygen and x varies locally along the surface of the membrane. 


for significant stage cut: our Xin 


This complicates the analysis, because now there will be a continuous decrease in x and L as the 
feed gas moves along the surface of the membrane toward the residue. 
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permeate 


permeate 


LLL LLL LLL 


LLL LLL LLL LL 


c E 
feed reject feed 


y 


FiGURE 13.3-4. Ideal flow patterns in a membrane separator for gases: (a) 
complete mixing, (b) cross-flow, (c) countercurrent flow, (d) 
cocurrent flow. 


Permeate 
purity, Yp 


Stage cut, 0 


FIGURE 13.8-2. Effect of stage cut and flow pattern on permeate purity. 
Operating conditions for air are as follows: xp = 0.209, a* = 
10, palpi = 380 cm Hgl76 cm Hg = 5. P4 = 500 x 107 
cm"(STP): cmls - cm? «cm Hg. (1) countercurrent flow, (2) 
cross-flow, (3) cocurrent flow, (4) complete mixing (W5). 
[Reprinted from W. P. Walawender and S. A. Stern, Sep. Sci., 
7, 553 (1972). By courtesy of Marcel Dekker, Inc.] 


Fig. 13.8-2 shows the effect of stage cut on the permeate concentration for each of the 4 flow 
configurations illustrated in Fig. 13.3-4. These are calculated by carving the membrane into a 
large number of differential slices and then integrating the resulting equations. 


Note that in general permeate concentration decreases with stage cut. This occurs because the 
local concentration on the feed side of the membrane is lower than at the inlet of the feed. 
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Chapter 16. Liquid-Liquid Extraction 


liquid-liquid extraction — separation of two components of a liquid (the *feed") by contact with 
a second immiscible liquid (the solvent") 


Extraction is used primarily when direct distillation is not economical (since the solvent usually 
has to be recovered by distillation). This might be the case with mixtures have very similiar 
boiling points, or substances that cannot withstand the temperatures of distillation, even under 
vacuum. Examples include: 


e extraction of penicillin from fermentation broth by contact with amyl or butyl acetate 


e recovery of acetic acid from dilute aqueous solutions by contact with ethyl acetate 
(fingernail polish solvent) or ethyl ether 


e separation of aromatics (rings) from aliphatics (straight chains) by contact with 
triethylene glycol 


e separation of high-MW fatty acids from vegetable oil by contact with liquid propane 
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Contacting of the two liquids might be performed 
in a countercurrent cascade of stages. Sometimes this 
can be accomplished in a sieve-tray tower like that 
used in distillation (see figure at right), where the 
difference in density drives the heavy fluid down while 
the lighter one floats upward. 


interface 


extract 


feed 


raffinate 


More often a mechanical agitator is required to form 
the emulsion of the two immiscible liquids at each 
stage. One common design is the mixer-settler 
(shown above). The two liquid phases are fed to the 
tank where the mixer imparts mechanical energy and 
thus disperses one of the phases as fine droplets in the 
second phases (usually continuous). This state of 
mixtures is called an emulsion. 


This emulsion is then transferred into a settler tank 
where the droplets begin to coalesces after the solute 
has redistributed across the liquid-liquid interface. 
Once the two phases have separated, the lighter one is 
drawn off the top (shown as "extract" in figure above) 
of the settler while the heavier one is drawn off the 
bottom (“raffinate”). 


Spring, 2001 


poor liquid out 


Principal 
interface 


e B 2,9 I 


Coalesced 
dispersed 


Perforated 
plate 


Downspout 
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Stage Stage 
2 3 


g -< 
Final 
extract 


Final 
roffinate 


Solvent 


Feed 


Fig. 10.5 Flow sheet of three-stage countercurrent mixer-settler extraction cascade. 


Mixers-settlers can also be arranged as a countercurrent cascade of stages, as shown in the figure 
above. 


TRIANGLE DIAGRAMS 


Many extraction systems involve three components which form two immiscible phases. For 
example, acetone and water are miscible in all proportions. When an organic solvent, methyl 
isobutyl ketone (MIK) is added to a water-rich acetone solution, two phases are formed: a water 
rich phase and an MIK-rich phase, both phases containing all three components. The 
composition of the two phases can be graphically represented by means of tie lines on a "triangle 
diagram." 


An equilateral triangle has the property that the sum of the A 
lengths of the three normals drawn from any interior point 

to each side equals the height of the triangle, regardless of 

which interior point is chosen. 


If we normalize the lengths of these three normals by the 
height, the normalized length can represent a mass fraction 
or a mole fraction of one of the three components. The 
geometry of the equilateral triangle assures that 


X4 t*tXptxc-l 


Thus each point inside the triangle represents the 

composition of a possible three-component mixture. A point at one of the three corners 
represents a pure component, while the normal distance between an interior point and the side 
opposite the corner for pure A represents the mole fraction or mass fraction of component A in 
the mixture. 
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The phase diagram for a mixture of acetone, 
methyl isobutyl ketone (MIK) and water at 25?C 
is shown at right. The shaded region (below the 


100 Acetone 


a" : AASA 
curve) represents compositions for which two KANAAT 
d pud So SOA d 
immiscible liquid phases form, while the eS IS 3 
E > á 
unshaded region (above the curve) are $ yy, Raffinote phases S M 
ege. . N XX 05 
compositions for which only one phase forms. e 5 EEO % 


: ; : : OEY 

The nearly horizontal lines with the shaded region COO 

represent tie lies connecting the compositions of x ene 
(ee COD. 


two liquid phases that co-exist at equilibrium. OE XS S 
XX ese NO B 


In this system, acetone is completely miscible in; ^ 4 os os or 06 os 04 03 o2 ol E 2 
both water and MIK, but MIK is only slightly REPRE 

soluble in water and water is only slightly soluble 

in MIK. This is an example of a Type I triangle 

diagram. 


100 Aniline 
A very different looking triangle diagram results P 
when one of the liquids is not completely miscible 
with either of other two. The diagram at right is 
for the system aniline-heptane- 
methylcyclohexane. This is an example of a Type 
IT triangle diagram. 


NOMENCLATURE 


Notice in the Type I triangle diagram that as ARN ANATAYA ATATA Ta 
the fraction of acetone gets higher, the tie lines get AIIIN IN inate phase SK KKK 

shorter. Eventually, the tie line becomes Pee ed ee E 
infinitesimally short, such that the composition of — /?/ BAPE RECON HERTANE (a) 


the two immiscible phases become equal at the 


plait point — point on the triangle 
diagram at which the composition of the two immiscible phases are 


equal. 
Under normal operation we contact the "feed" stream 
; rs s : Sa wah ae Solvent Extract 
with a “solvent” stream, producing two immiscible liquid MIK 
phases which are separated to form two streams, which 
are called the “extract” and the “raffinate.” 
Raffinate Feed 


extract — product stream which is richer in the solvent 
(the remains of the solvent stream after the 
solute has joined it) 


acetone + water 
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raffinate — product stream which 1s leaner in the solvent (what remains of the feed after the 
solute has been extracted) 


C (eth: 
10 a er) 


RECTANGULAR DIAGRAMS 


tie line 


2 extract layer, yc vs. Ya 
While triangle diagrams are commonly used to 5 
d Oo two-phase region 
present the phase diagram for three-component g 0.5 
mixtures, the geometry of an equilateral triangle E one-phase region 
complicates construction. Three-component E | 
s $ = i à 
phase behavior can also be presented using a aiina ei le YS tA 
: : : B H4 i A (acetic acid) 
rectangular diagram, an example of which is (water) 0 05 1.0 
shown at right for a ternary system of Type I. Mass fraction A (x4, 74) 
3 
As with an (equilateral) triangle, the point at each ERN 
apex of this (right) triangle represents one of the E 
three components in the pure state. Each of the PE 
& 8 


three sides of the triangle represents mixtures in 
which one of the three components is completely 

absent. For example, the base of the right triangle — "9" 253. eme 20 (ver r ether (C) lid uid phase 
at right represents mixtures of acetic acid and 

water which do not contain any ether. 


Water layer (raffinate) composition, x 4 


The hypotenuse of the right triangle corresponds to 
x4 recall or xg=0 
where the subscripts denote components. Lines parallel to the hypotenus correspond to 
X4+xXc=const<1 or xp>O 
if the line lies below the hypotenus and 
X4+xc=const>1 or xpg«0 
if the line lies above the hypotenus. Clearly, the latter situation is impossible so that points 


outside the right triangle are extraneous. 


FLASH MIXING OF Two LIQUIDS 


Suppose we take two immiscible liquids (formed of 
the same three components: A, B and C) and mix them. 


js L, X4, X¢ 
What will be the composition of the mixture? Let 
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L, V — total mass flowrate of the two liquid input streams 
M - total mass flowrate of the output stream (ignoring the phase distribution) 
X4, xc = mass fraction of components A,C in the stream having flowrate L 
Y4, yc = mass fraction of components A,C in the stream having flowrate V 
X4M> XCM = mass fraction of components A,C in the stream having flowrate M 
EXAMPLE #1: Given: L, V, x4, xc, y4 and yc 
Find: M, x44 and xcy 


Solution: Here is a receipe for a graphical solution using a 


rectangular diagram (see figure at right). = 
o 
Step 1: Locate points L and V on the rectangular diagram E m 
[coordinates are (x4, xc) and (v4, yc)] "i 
= yc 
Step 2: Locate point M on a straight line connecting L and V 
using the inverse lever-arm rule: JA XAM XA 


Mass fraction A 
L 
M = 1 
LM y ey 


where LM or VM denote the length of the lines connecting points L and M or V and M, 
respectively. 


Step 3: Read coordinates of point M as (x4), xc). 
Proof: A total mass balance is 
L+V=M (157) 


and mass balances on components A and C yield: 


Lx 4+ Vy4 = Mx4y (158) 
Lxç+ Vyc = Mxcm (159) 
Substituting (157)|into Lx4 + Vy4 = (L-V)x4M 
Collecting terms, we obtain L(x4—xa4y) = Vam- ya) 
or L YA 4M (160) 


V XAM —%XA 
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Similarly combining and we obtain 


L Yc-XcM (161) 
V XCM —Xc 


Equating L66ylana tress) YA-*XAM _YC7~*CM 


XAM~*A %*CM~%*C 


X = XOX 

" CM JG a *Q- XM aua 
XAM YA X4 —XAM 
1 r 


slope of line VM slope of line LM 


This shows that line segments LM or VM have the same slope and therefore point M must lie 


on line LV . Moreover, an be rewritten as 
XAM TXY 7^4 TAF 
= dno. Vi 
V XA—XAM LM 


sin 8 


which proves 


SINGLE-STAGE EXTRACTOR 


If point M in the previous example lies in the 1-phase region of the phase diagram, then we 
are done because we know the composition (x 4,4, xcm) of the single phase leaving the mixer. If 
the point M lies in the two-phase region and the two phases are allowed to reach equilibrium and 
separate in the settler, we will have (not one, but) two product streams. 


Let's re-label the two input and two output streams as in the figure y, Va 
at right. In a typical extraction problem, stream Lp is an aqueous 
solution containing a solute (component A) we are trying to extract by 
contact with an organic solvent (component C). So let's assume that the Ło Li 
second feed stream (labelled V5) is pure solvent. Let's try to determine 
the composition of the two product streams. 


EXAMPLE #2: Given: Lo; V5, X 40> XCO 7 0, YA2 = 0 and yc =] 


Find: Lj, Vi, x 441, xcj, y41 and yc} 
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Solution: Here is a receipe for a graphical solution using a 
rectangular diagram (see figure at right). 


Step 1: Locate points Lọ and V3 on the rectangular diagram 
[coordinates are (x 49, xco=0) and (y 4570, yc2=1)] 


Step 2: Locate the mixing point M as in Example #1 


Step 3: Find the tie line which passes through M 


Step 4: The endpoints of this tie line represent L, and V 
[coordinates are (x41, xc) and (v41, yc1)] 


Step 5: The flowrates L4 and V, can be determined using the reverse of the procedure for 
determining the mixing point M in Example 71 


Essentially the mixing point gives the composition of a hypothetical single-phase mixture formed 
from the two feed streams. According to a mass balance, this is the same mixture which would 
be formed by mixing the two product streams. In other words, M is the mixing point for either 
Lo and V5 or L; and Vj. 


MULTI-STAGE EXTRACTOR 


If we are satisfied with the amount of solute we have been able to extract from the solution in 
the single-stage extractor, then we are done. If we desire to extract more of the solute (i.e. we 
wish to achieve a smaller value of x41) then multiple stages might be used. Let's specify an 
arbitrarily small value for x 4y, the concentration of solute remaining after N stages. 


extract solvent 
Vne1 Vy VNy+1 


Ly Ly. I LN 
raffinate 


EXAMPLE #3: Given: Lo, VN+1> X40» XCO = 0, YAN] 7 0, YCN+1 7 1 and XAN 


Find: Ly, Vi, xc Vat and yc} 


06-361 page 192 Spring, 2001 
Solution: 
Step 1: Locate points Lp and Vy. 


Step 2: Locate their mixing point M. Read coordinates (x4, 
xcy) and calculate M = Lo + Vy. 


This will also be the mixing point for the two product streams Ly 
and V4, although these streams are not coming from the same 
stage (ie. they contain the same material as in the two feed 
streams). Assuming all the stages in the cascade are equilibrium 
stages, we also know that the points Ly and V, lie somewhere on 
the equilibrium curve, but not necessarily on opposite ends of the 
same tie line (since these streams are not coming from the same stage, they don't have to be in 
equilibrium with each other). 


Step 3: Using the known value of x4w, locate the point Ly on the lower portion of the 
equilibrium curve (i.e. on the raffinate portion of the curve). Read xcy. 


Step 4: Draw a straight line connecting Ly and M; then extend this line until it crosses the 
equilibrium curve again. 


Step 5: The intersection of this line with the equilibrium curve is point V4. Read the coordinates: 
O41, yc). 


Step 6: Calculate the flowrates Ly and V, as the simultaneous solution of the following two 


equations: 
total mass balance: Lyt =M 
component C balance: Lyxcy + Vivo, = Mxcm 


Everything in these equations is known except for the two unknowns: Ly and V. 

EXAMPLE #4: How many equilibrium stages are required to perform the separation described 
in Example #3? In other words, what is the value of 
N? 


Solution: 


Step 1: Perform steps 1-6 in example #3 to locate point V1. tie line 


Step 2: Locate point L; at the other end of the tie line passing through 
Vi 


After we remove extraneous lines and markings (for clarity), the L, 
rectangular diagram looks as shown at right. 
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Step 3: Knowing the flowrates and compositions of the two output streams for stage #1, we can 
calculate the mixing point M, using the method of Example #1. This is also the mixing point 


for the two feed streams: Lo and V5. 


Step 4: Locate the point V by connecting Lo and M, with a straight line 
and extending it until it crosses the equilibrium curve. 


es 


Step 5: Locate L, at the other end of the tie line passing through V». tie line 

Step 6: Repeat steps 3-5 (after incrementing all the subscripts by +1) 
until point Ly has an x-coordinate x4, which equals or is smaller 
than that specified as x 4y for the entire cascade. 


DIFFERENCE-POINT METHOD 


Although the procedure in Example #4 works well (provided solvent flowrate Vj, is above 
the minimum required), the coordinates of the mixing points Mj, M3 etc. have to be calculated 
manually. This tedious step can be avoided by observing that all the “operating lines” (i.e. L9 V», 
L, V3 etc.) have a common point called the difference point. 


Consider a total mass balance around some arbitrary stage in the y, Fati 
cascade: 


Lii + Va =La + Vn 


L n-1 Ln 
which can be rearranged to obtain 


Ln-Vn=Ly- Viti =A (162) 


This difference between the flowrates of two “adjacent” streams in the cascade has exactly the 
same value for every stage in the cascade (i.e. n = 1,2 ... N). We denote this difference in 
flowrates by A. If instead of a total mass balance, we perform a mass balance on component A: 


Ln-1XAn-1 + Vn+1YAn+1 = LnXAn t VnYAn 


which can be rearranged to obtain 


Ln-1X4n-1 7 Vn An = Ln¥An 7 Vint Ant = AXAA (163) 


Once again the difference the mass flowrates of component A in any two adjacent streams is 
exactly the same, regardless of which stage we pick. We denote this difference as Ax44: the 
product of A and a ficticious mass fraction which we call x44. The value of this ficticious mass 
fraction can be calculated by dividing|(163) by (162) for n=1 (choosen because we know the 
flowrates and compositions of the two adjacent streams at the end of the cascade from Example 
#2): 
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-Liox4 -hya Qu re Loxco - ya (164) 


X 
xis Lo- Lo- 


An equation similar to this can be obtained for component C, which leads to the second equation 
above. The point (x44, xca) is called the difference point and is denoted by A. All the operating 
lines pass through this common point. 


EXAMPLE #5: Repeat Example #4 using the difference- 
point method A«— operating point 


Solution: VN 


Step 1: Locate points Lọ and A. The coordinates of the 
latter are calculated from 


XC; Jc 


Step 2: Connect Lo and A with a straight line. 


Step 3: Locate V, as the point where this line crosses the 
(upper part of) the equilibrium curve 


Step 4: Determine the tie line passing through V4. Locate Lu 
L4 as the other end of this tie line. 


XA» JA 


Step 5: Connect L, and A with a straight line. 


Step 6: Locate V» as the point where this line crosses the 
(upper part of) the equilibrium curve 


Step 7: Determine the tie line passing through V». Locate L, as the other end of this tie line. 
Step 8: Repeat steps 5-7 until the desired composition is obtained for Ly. 


Proof (of difference point method): 
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